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Abstract 

The  control  of  aerial  gymnastic  maneuvers  is  challenging  because  these  maneuvers 
frequently  involve  complex  rotational  motion  and  because  the  performer  has  limited 
control  of  the  maneuver  during  flight.  A  performer  can  influence  a  manuever  using  a 
sequence  of  limb  movements  during  flight.  However,  the  same  sequence  may  not 
produce  reliable  performances  in  the  presence  of  off-nominal  conditions.  How  do 
people  compensate  for  variations  in  performance  to  reliably  produce  aerial 
maneuvers?  In  this  report  I  explore  the  role  that  passive  dynamic  stability  may  play 
in  making  the  performance  of  aerial  maneuvers  simple  and  reliable. 

I  present  a  control  strategy  comprised  of  active  and  passive  components  for 
performing  robot  front  somersaults  in  the  laboratory.  I  show  that  passive  dynamics 
can  neutrally  stabilize  the  layout  somersault  which  involves  an  “inherently  unstable” 
rotation  about  the  intermediate  principal  axis.  And  I  show  that  a  strategy  that  uses 
open  loop  joint  torques  plus  passive  dynamics  leads  to  more  reliable  11/2  twisting 
front  somersaults  in  simulation  than  a  strategy  that  uses  prescribed  limb  motion. 

Results  are  presented  from  laboratory  experiments  on  gymnastic  robots,  from 
dynamic  simulation  of  humans  and  robots,  and  from  linear  stability  analyses  of 

these  systems. 
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Abstract 

The  control  of  aerial  gymnastic  maneuvers  is  challenging  because  these  maneuvers 
frequently  involve  complex  rotational  motion  and  because  the  performer  has  limited 
control  of  the  maneuver  during  flight.  While  a  performer  can  execute  a  sequence  of 
limb  movements  during  flight  to  produce  a  desired  maneuver,  the  same  sequence  may 
not  produce  reliable  performances  in  the  presence  of  off-nominal  conditions.  How 
do  people  compensate  for  variations  in  performance  conditions  to  reliably  produce 
aerial  maneuvers?  It  is  possible  that  people  sense  errors  in  performance  and  actively 
compute  appropriate  responses  to  them  to  produce  reliable  maneuvers.  However, 
it  is  also  possible  that  the  maneuver  is  inherently  stable,  that  the  body  naturally 
compensates  for  variations,  and  that  the  athlete  does  little  active  computation. 

In  this  thesis  I  explore  the  role  that  passive  dynamic  stability  may  play  in  making 
the  performance  of  aerial  maneuvers  simple  and  reliable.  I  consider  the  control  of 
the  tucked  somersault,  the  layout  somersault,  and  the  11/2  twisting  front  somer¬ 
sault.  I  present  a  control  strategy  comprised  of  active  and  passive  components  for 
performing  robot  front  somersaults  in  the  laboratory.  I  show  that  passive  dynamics 
can  neutrally  stabilize  the  layout  somersault  which  involves  an  “inherently  unstable” 
rotation  about  the  intermediate  principal  axis.  And  I  show  that  a  strategy  that  uses 
open  loop  joint  torques  plus  passive  dynamics  leads  to  more  reliable  11/2  twisting 
front  somersaults  in  simulation  than  a  strategy  that  uses  prescribed  limb  motion.  Re¬ 
sults  are  presented  from  laboratory  experiments  on  gymnastic  robots,  from  dynamic 
simulation  of  humans  and  robots,  and  from  linear  stability  analyses  of  these  systems. 
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Chapter  1 


Introduction 


In  both  the  1992  and  1994  winter  Olympics,  the  hnal  jump  performed  by  the  men’s 
champion  freestyle  skier  was  a  quadruple  twisting,  triple  somersault.  In  the  per¬ 
formance  of  these  ’jumps’  the  athletes  skied  off  of  a  ramp  from  which  they  soared 
approximately  45  feet  into  the  air,  remained  aloft  for  approximately  3.0  sec,  rotated 
three  times  about  a  horizontal  axis  and  four  times  about  their  body  vertical  axis,  and 
hnished  by  landing  squarely  on  their  feet  so  that  they  could  continue  skiing  down  the 
hill  in  a  controlled  fashion.  It  is  incredible  that  these  athletes  can  perform  a  maneuver 
like  this  with  such  accuracy.  How  do  they  do  it? 

Two  related  issues  that  make  these  performances  challenging  are  the  controllabil¬ 
ity  and  stability  of  aerial  maneuvers.  Controllability  refers  to  the  ability  of  an  athlete 
to  influence  the  outcome  of  an  aerial  maneuver  once  it  has  been  initiated.  An  airborne 
performer  can  not  apply  any  external  forces  or  torques  to  the  body.  So  how  can  a 
performer  control  his  body  orientation?  Previous  research  of  this  topic  has  revealed 
movement  techniques  that  performers  can  use  to  produce  complicated  aerial  maneu¬ 
vers.  These  techniques  involve  moving  the  limbs  to  reconhgure  the  body  during  flight 
(Figure  1-1).  However,  while  a  prescribed  set  of  body  conhgurations  can  produce  a 
desired  aerial  maneuver,  this  sequence  may  not  lead  to  reliable  performances. 

The  stability  of  aerial  maneuvers  concerns  their  reliable  performance  in  the  pres- 


12 


CHAPTER  1.  INTRODUCTION 


13 


•  Change  inertia 
about  an  axis 


From  Frohlich  1980. 


Figure  1-1:  A  limited  ability  to  influence  the  outcome  of  a  ballistic  maneuver  arises 
from  the  relative  movement  of  the  limbs  and  torso  during  flight.  One  way  to  influence 
the  outcome  is  to  change  the  moment  of  inertia  about  an  axis  in  order  to  change  the 
rotation  rate  about  that  axis.  The  standard  being  the  ice-skater’s  spin.  A  second  way 
to  influence  the  outcome  is  to  use  momentum- free  rotations  [Smith  67,  Frohlich  79]. 
This  technique  allows  a  structure  to  be  reoriented  while  maintaining  zero  angular  mo¬ 
mentum  by  performing  a  sequence  of  limb  movements.  A  third  way  is  to  reconhgure 
the  system  so  that  the  principle  axes  of  inertia  are  reoriented  relative  to  the  inertially 
hxed  angular  momentum  vector.  This  allows  sharing  of  momentum  between  principal 
axes.  This  procedure  is  frequently  used  to  introduce  or  remove  twist  in  somersaults 
[Batterman  68,  Frohlich  79,  Yeadon  84].  The  existence  of  these  mechanisms  makes 
it  possible  to  actively  adjust  the  outcome  of  an  aerial  maneuver  once  it  has  been 
initiated. 
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ence  of  off-nominal  conditions.  Stability  is  important  because  the  accumulation  of 
even  small  errors  during  the  relatively  long  flight  times  of  an  aerial  maneuver  could 
result  in  a  disastrous  outcome.  Making  a  performance  reliable  requires  that  the  ath¬ 
lete  compensate  for  inaccuracies  in  movement,  variations  in  equipment,  and  external 
disturbances.  It  is  possible  that  people  produce  reliable  maneuvers  by  sensing  these 
variations  and  actively  computing  appropriate  responses  to  compensate  for  them. 
However,  the  complexity  of  this  feedback  control  approach  would  appear  to  place 
great  demands  on  the  athlete’s  perceptive  and  motor  control  abilities.  Is  there  a 
more  simple  approach  to  producing  reliable  maneuvers  that  does  not  require  active 
compensation  by  the  athlete? 

The  focus  of  this  thesis  is  on  the  use  of  passive  dynamic  stability  as  an  alternative 
or  a  complement  to  active  control  for  producing  reliable  aerial  maneuvers.  Passive 
dynamic  stability  means  that  a  maneuver  is  inherently  stable  by  virtue  of  the  natural 
dynamic  interaction  of  the  limbs  and  body.  The  precise  limb  movement  that  a  per¬ 
former  makes  during  a  maneuver  will  depend  upon  not  only  his  motor  activity  but 
also  on  environmental  forces.  Is  it  possible  that  the  passive  forces  that  arise  due  to  off- 
nominal  conditions  could  provide  a  built-in  correction?  If  so,  passive  behavior  could 
automatically  compensate  for  errors  in  initial  conditions,  in  control  movements,  or 
from  external  disturbances.  If  this  were  the  case  then  the  demands  for  active  control 
by  the  athlete  could  be  dramatically  reduced.  A  passive  control  strategy  is  appealing 
because  it  could  relieve  the  human  performer  of  sensing  small  variations  in  movement 
and  computing  control  responses  fast  enough  to  produce  accurate,  stable  maneuvers. 
Instead,  the  athlete  may  need  only  to  “play  back”  a  pre-recorded  set  of  motor  actions. 
This  feed  forward  command  combined  with  the  passive  dynamic  response  of  his  body 
may  allow  the  maneuver  to  “unfold”  on  its  own. 

To  see  if  passive  dynamic  stability  could  play  a  role  in  aerial  performances,  I 
consider  the  control  of  three  gymnastic  maneuvers,  the  tucked  front  somersault,  the 
back  layout  somersault,  and  the  front  somersault  with  one  and  a  half  twists.  I  show 
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that  passive  dynamics  can  play  a  significant  role  in  making  these  maneuvers  reliable. 
Neutrally  stable  passive  dynamics  make  the  control  of  the  tucked  somersault  simple 
since  an  active  controller  need  only  consider  rotation  about  the  major  principal  axis. 
The  layout  or  straight  body  somersault  is  considered  to  be  inherently  unstable  because 
it  requires  rotation  about  the  middle  principal  axis  of  inertia,  an  unstable  equilibrium 
for  a  rigid  body.  I  show  that  the  layout  somersault  can  passively  be  made  neutrally 
stable  by  tuning  the  compliance  of  the  shoulders  of  the  performer.  Finally,  I  present 
results  that  suggest  that  a  passive,  compliant  model  of  the  human  body  that  uses 
open  loop  torque  control  can  produce  more  reliable  11/2  twisting  somersaults  than 
a  model  that  uses  prescribed  limb  motion.  I  obtained  these  results  using  simple 
analytic  models,  non-linear  dynamic  simulation,  and  laboratory  robots.  Next,  I  briefly 
describe  the  results  of  experiments  on  each  of  these  maneuvers. 

The  somersault  is  a  maneuver  in  which  a  performer  jumps  into  the  air  and  rotates 
once  about  a  side-to-side  axis  before  landing  on  the  ground.  The  main  requirement  for 
the  somersault  is  to  land  in  a  balanced  manner.  This  in  turn  requires  the  performer  to 
land  with  a  precise  body  attitude.  The  tucked  somersault  exhibits  passive  directional 
stability  in  rotation  since  it  involves  rotation  about  the  major  principal  axis  of  inertia. 
This  means  that  imprecise  initiation  of  a  somersault  will  not  dramatically  affect  the 
orientation  of  the  spin  axis  during  flight.  However,  avoiding  over-rotation  or  under¬ 
rotation  of  a  somersault  about  the  spin  axis  may  require  compensation  of  rotation 
rate  during  flight.  Active  control  of  the  inertia  can  be  used  to  correct  errors  in 
somersault  rotation  rate.  I  present  results  from  somersault  experiments  using  a  60 
lb,  one  meter  tall  laboratory  robot  that  runs  on  two  springy  legs  (Figure  1-2).  The 
robot  was  programmed  to  initiate  the  somersaults  using  a  pre-programmed  pattern 
of  action.  To  avoid  large  tilt  and  twist  angles  of  the  3D  Biped  during  take-off  and 
landing  we  use  a  wide  double  stance  of  the  robot  and  insure  that  the  feet  touch  down 
simultaneously.  During  flight  the  robot  actively  controls  rotation  rate  by  “tucking” 
or  “untucking”  its  legs  to  manipulate  the  robots  inertia.  The  robot  actively  adjusts 
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Figure  1-2:  Photograph  of  the  3D  Biped  Robot. 

the  position  of  its  feet  relative  to  its  center  of  mass  just  prior  to  landing.  This  element 
was  necessary  to  compensate  for  errors  in  the  estimated  landing  time  of  the  robot. 
On  its  best  day  the  robot  did  successful  somersaults  and  continued  running  on  7  out 
of  10  attempts. 

The  layout  or  straight  body  somersault  is  considered  to  be  inherently  unstable 
because  it  requires  rotation  about  the  middle  principal  axis  of  inertia,  an  unstable 
equilibrium  for  a  rigid  body.  A  rigid  body  that  is  somersaulting  about  the  middle 
principal  axis  will  always  exhibit  a  sequence  of  half  twists  about  the  minor  principal 
axis.  Despite  this  fact,  athletes  regularly  perform  this  maneuver  with  apparent  ease. 
Previously,  biomechanics  researchers  have  assumed  that  the  athlete  senses  the  insta¬ 
bility  of  the  maneuver  and  actively  compensates  for  it  with  movements  of  the  arms 
and  body  [Nigg  74,  Hinrichs  78,  Yeadon  90]. 

I  show  that  the  layout  somersault  can  be  a  passive,  neutrally  stable  maneuver. 
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Figure  1-3:  Photograph  of  the  mechanical  doll  used  for  experiments  with  layout 
somersaults.  The  doll  could  consistently  perform  stable,  triple  layout  somersaults. 

Passive  stabilization  of  the  layout  somersault  results  from  the  natural  dynamic  inter¬ 
action  of  the  limbs  and  body  during  movement.  Stabilization  arises  from  the  inherent 
tendency  of  the  arms  to  tilt  in  response  to  twisting  movement  of  the  body.  The  arm 
tilt  forces  the  principal  axes  of  the  system  to  move  in  a  direction  that  compensates 
for  tilt  and  twist  errors.  This  built-in  correction  eliminates  the  divergent  tendency  of 
the  system  as  long  as  the  compliance  of  the  shoulders  cancels  the  unstable  centrifu¬ 
gal  forces  on  the  arms.  I  verify  this  result  with  linear  stability  analysis,  non-linear 
dynamic  simulation,  and  experiments  on  a  human-like  doll  built  and  tested  in  the 
laboratory  (Figure  1-3.)  The  doll  has  spring-driven  arms  but  has  no  other  control 
system,  sensors,  or  actuators.  The  doll  routinely  exhibits  triple  somersaults  about  its 
middle  principal  axis  without  twisting. 

While  twisting  is  to  be  avoided  in  the  layout  somersault,  it  is  a  feature  in  other 
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Figure  1-4:  Images  arranged  in  right-to-left,  top-to-bottom  order  from  a  dynamic 
simulation  of  a  1  1/2  twisting  front  somersault.  We  found  that  open  loop  torque 
control  of  a  compliant  model  led  to  relatively  reliable  performances. 

tricks.  I  present  results  from  dynamic  simulation  of  a  human  performing  a  front 
somersault  with  one  and  a  half  twists.  The  twisting  maneuver  is  started  from  a  front 
somersault  by  reconhguring  the  body  mid-maneuver.  While  a  sequence  of  prescribed 
limb  movements  can  be  found  to  produce  a  twisting  maneuver,  simulation  results 
show  that  the  performance  is  sensitive  to  small  variations  in  initial  conditions.  If, 
instead,  open  loop  torque  control  is  used  with  a  compliant,  passive  dynamic  model  of 
the  human  the  maneuver  reliability  can  be  signihcantly  improved. 

I  also  present  results  from  dynamic  simulation  of  a  front  somersault  with  one  half 
twist  performed  by  the  3D  Biped.  We  found  that  to  make  the  simulated  3D  Biped 
perform  a  satisfactory  front  somersault  with  one  half  twist  we  had  to  add  weight  to 
the  robot  to  make  its  inertia  more  like  that  of  a  human. 

The  ultimate  goal  of  this  research  is  to  develop  a  theory  of  passive  dynamic. 
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aerial  maneuvers.  By  showing  that  passive  dynamics  can  play  a  signihcant  role  in 
the  performance  of  reliable  aerial  maneuvers  we  argue  that  they  reduce  the  need  for 
active  control  from  the  athlete.  The  performance  of  a  maneuver  may  be  simplihed 
by  relegating  some  responsibility  for  control  to  the  natural,  mechanical  behavior  of 
the  body.  The  appropriate  use  of  passive  dynamics  may  have  the  added  beneht 
of  producing  natural  looking,  coordinated  movement.  Perhaps  athletes  and  other 
people  use  a  performance  strategy  that  seeks  to  maximize  passive  dynamic  behavior. 
It  is  difficult  to  know  what  control  strategies  people  may  or  may  not  use,  but  in 
the  laboratory  we  can  examine  the  feasibility  of  a  strategy  by  testing  it  in  a  real  or 
simulated  system. 

1.1  Background 

Relevant  background  material  for  the  study  of  gymnastic  maneuvers  comes  from 
holds  such  as  biomechanics,  biology,  robotics,  aeronautics,  and  astronautics.  Several 
researchers  have  explicitly  studied  the  performance  of  gymnastic  maneuvers.  These 
studies  have  revealed  the  salient  features  of  known  gymnastic  techniques  for  produc¬ 
ing  maneuvers.  Some  robotics  researchers  have  studied  gymnastic  maneuvers  using 
dynamic  simulation  and/or  laboratory  robots  in  order  to  develop  strategies  for  con¬ 
trol.  The  study  of  passive  dynamic  stability  has  roots  in  the  helds  of  aeronautics 
and  astronautics.  It  is  common  for  airplanes  and  spacecraft  to  be  designed  to  exhibit 
passive  dynamic  stability.  The  study  of  human  locomotion  provides  an  inspirational 
example  of  passive  dynamic  stability  and  the  rich  behavior  that  it  can  produce.  In 
the  following  sections,  I  briehy  discuss  background  material  from  each  of  these  helds. 

1.1.1  Gymnastic  Maneuvers 

Most  researchers  investigating  the  control  of  aerial  maneuvers  have  been  concerned 
with  explaining  the  physics  of  the  maneuvers.  During  the  aerial  phase  of  a  maneuver 
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the  performer  can  not  produce  any  external  forces  on  the  system  so  the  trajectory  of 
the  center  of  mass  and  the  angular  momentum  are  hxed  from  the  moment  of  take¬ 
off  to  the  moment  of  landing.  Knowing  this,  early  researchers  were  compelled  to 
hnd  explanations  of  two  aerial  maneuvers  that  at  hrst  appeared  to  violate  conser¬ 
vation  of  angular  momentum:  1)  a  cat  when  dropped  with  no  net  angular  momen¬ 
tum  will  right  itself  before  landing  on  the  ground  [Marey,  Kane  69]  and  2)  spring 
board  divers  who  leave  the  board  with  rotation  only  about  a  side-to-side  (somer¬ 
sault)  body  axis  can  subsequently  initiate  rotation  about  their  head-to-toe  (twist) 
body  axis  [Batterman  68,  Frohlich  79,  Yeadon  84],  Where  did  the  extra  angular  mo¬ 
mentum  come  from?  Resolving  these  apparent  discrepancies  led  researchers  to  hnd 
movement  techniques  that  could  be  used  to  perform  these  and  other  interesting  ma¬ 
neuvers.  These  techniques  involved  reconhguring  the  body  in  hight  (Figure  f-f). 

Takashima  [Takashima  90]  studied  high  bar  maneuvers.  He  used  the  control  of 
rotation  rate  about  the  somersault  axis  to  produce  a  balanced  landing  of  a  simulated 
human  dismounting  from  the  horizontal  bar.  His  algorithm  could  produce  accurate 
tucked,  multiple-somersault  dismounts.  In  hight,  he  used  the  feedback  control  of 
posture  to  control  rotation  rate.  He  used  a  combination  of  feed  forward  and  feed 
back  control  to  execute  the  landing. 

While  they  appear  to  be  similar  maneuvers,  the  tucked  somersault  and  the  layout 
somersault  are  dynamically  very  different.  The  tucked  somersault  involves  rotation 
about  the  major  principal  axis  of  inertia,  a  stable  mode  of  rigid  body  rotation.  Rota¬ 
tion  about  the  major  and  minor  principal  axes  of  inertia  exhibits  a  limited  form  of  sta¬ 
bility  called  directional  stability.  Directional  stability  refers  to  the  fact  that  the  spin 
axis  will  maintain  roughly  the  same  inertial  orientation  when  deviated  slightly  from 
that  orientation.  The  layout  somersault  involves  rotation  about  the  intermediate  prin¬ 
cipal  axis  of  inertia,  an  unstable  mode  of  rigid  body  rotation  [Crandall  68,  Hughes  86]. 
When  a  rigid  body  is  spun  about  its  intermediate  axis,  that  body  axis  will  exhibit 
large  excursions  from  its  initial  orientation.  These  cyclic  excursions  are  a  series  of 
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nearly  180°  rotations  away  from  and  returning  to  the  initial  axis  orientation.  The 
different  modes  of  rotation  can  be  easily  demonstrated  to  the  reader  by  spinning  a 
video  tape  or  a  book  (with  a  rubber  band  around  it)  about  each  of  its  axes  of  sym¬ 
metry.  You  will  see  that  when  the  body  is  spun  about  the  major  and  minor  principal 
axes  the  spin  axis  roughly  maintains  its  inertial  orientation.  However,  when  the  body 
is  spun  about  its  middle  principal  axes  it  will  exhibit  a  sequence  of  twists  about  the 
long  body  axis.  This  rigid  body  instability  has  lead  researchers  to  conclude  that  the 
layout  somersault  is  inherently  unstable. 

Nigg  recognized  that  the  layout  somersault  may  be  unstable  since  it  involved 
rotation  about  the  middle  principal  axis  [Nigg  74].  Using  cinematographic  techniques 
Hinrichs  [Hinrichs  78]  measured  the  body  conhguration  of  an  athlete  performing  a 
layout  somersault.  He  conhrmed  that  the  layout  somersault  involved  rotation  about 
the  middle  principal  axis  of  inertia.  He  hypothesized  that  the  athlete  made  small 
corrective  movements  of  the  arms  and  torso  in  flight  to  stabilize  the  somersault. 

Yeadon  [Yeadon  84]  used  a  combination  of  cinematography  and  dynamic  simula¬ 
tion  to  study  the  control  of  aerial  maneuvers.  In  his  research  he  developed  a  mass 
properties  model  of  the  human  form  that  provided  an  estimate  of  inertial  parameters 
from  anthropometric  measurements.  He  also  hlmed  highly  skilled  athletes  performing 
complex  aerial  maneuvers  and  digitized  this  data  to  determine  the  body  attitude,  con¬ 
hguration,  and  angular  momentum  during  hight.  Then  he  used  the  inertia  parameters 
and  digitized  conhguration  and  momentum  data  as  input  to  a  dynamic  simulation  of 
the  human  body  during  hight.  This  dynamic  model  computed  the  gross  body  atti¬ 
tude  during  hight  as  a  function  of  the  measured  internal  conhguration  and  angular 
momentum.  Using  this  system,  Yeadon  could  numerically  study  the  effect  of  changes 
in  body  conhguration  on  the  performance  of  complex  aerial  maneuvers. 

Yeadon  found  that  by  piking  or  arching  (bending  at  the  waist  in  the  sagittal 
plane)  during  a  layout  somersault  an  athlete  can  change  his  or  her  inertia  enough 
to  make  the  somersault  axis  an  axis  of  maximum  inertia.  This  in  turn  implies  that 
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arched/piked  somersaults  are  passively  stable.  Because  many  athletes  perform  layout 
somersaults  with  arch,  it  is  possible  that  this  is  done  to  stabilize  the  maneuver. 

However,  in  Yeadon’s  analysis  of  a  double  layout  somersault  performed  by  Carl 
Furrer,  the  1982  World  Trampoline  Champion,  he  found  that  Furrer  was  in  fact  ro¬ 
tating  about  the  middle  principal  axis  of  inertia  during  nearly  all  of  the  maneuver. 
Furthermore,  dynamic  simulation  of  this  layout  somersault  exhibited  the  character¬ 
istic  twist  instability  of  rotation  about  the  middle  principal  axis  while  the  actual 
human  performance  exhibited  no  such  instability.  This  difference  implied  that  the 
small  digitization  errors  in  translating  him  conhguration  data  to  the  simulation  were 
responsible  for  the  change  in  performance,  a  fact  that  would  point  to  an  inherently 
unstable  system.  These  results  suggest  that  the  athlete  uses  some  form  of  stabilization 
to  perform  the  layout  somersault. 

Yeadon  proposed  a  specihc  technique  for  stabilizing  the  layout  somersault.  He 
realized  that  asymmetric  movement  of  the  arms  in  the  frontal  plane  could  be  used  to 
stabilize  rotation  about  the  middle  principal  axis  of  inertia.  He  designed  a  stabilizing 
feedback  controller  that  used  the  sensed  twist  angle  as  the  feedback  signal  to  drive 
the  arm  abduction/adduction  angular  rates.  Using  a  linear  model,  Yeadon  found 
the  athlete  would  have  to  respond  to  a  growing  twist  instability  within  0.28  of  a 
somersault  (sa  200ms)  [Yeadon]  in  order  to  maintain  stability. 

The  ballistic  nature  of  many  gymnastic  performances  makes  an  open-loop  (feed 
forward)  component  a  likely  part  of  any  control  strategy.  An  open  loop  strategy  is 
one  in  which  the  performer’s  control  motions  are  not  derived  from  the  current  state 
of  motion  but  are  instead  produced  from  a  pre-programmed  pattern  of  action.  Since 
important  parameters  of  ballistic  motion  are  hxed  from  takeoff,  an  open-loop  strategy 
is  required  to  anticipate  the  maneuver  in  order  to  set  up  appropriate  initial  conditions. 
Raibert  and  Hodgins  programmed  a  planar  biped  robot  to  perform  front  somersaults 
in  the  laboratory  using  a  combination  of  open  loop  and  feedback  control  strategies. 
The  biped  robot  used  an  open  loop  strategy  for  initiation  of  the  somersault  and  the 
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majority  of  flight.  A  feedback  controller  was  used  to  accurately  place  the  feet  just 
prior  to  touch  down.  Using  this  combination  of  feed  forward  and  feedback  control,  the 
robot  could  successfully  perform  the  front  somersault  followed  by  balanced  running 
on  90%  of  the  trials. 

1.1.2  Aeronautics,  Astronautics  and  Celestial  Mechanics 

Some  modern  aircraft  are  designed  to  require  active  control  for  stabilization.  These 
aircraft  use  digital  computers  to  manipulate  the  aircraft  control  surfaces  to  render  the 
planes  flyable.  This  inherent  instability  is  tolerated  to  provide  highly  maneuverable 
aircraft.  However,  passive  dynamic  stability  is  commonly  built  into  general  aviation 
aircraft  and  spacecraft.  More  precisely,  they  have  an  equilibrium  condition  like  flying 
straight  and  level  for  an  aircraft,  or  spinning  about  the  major  principal  axis  for  a 
satellite  that  is  stable  so  that  the  craft  can  tolerate  external  disturbances  without 
diverging  from  the  stable  equilibrium  condition.  The  origin  of  this  topic  has  its  roots 
in  celestial  mechanics. 

The  moon  always  presents  the  same  face  to  the  Earth.  However,  it  does  not  do 
so  exactly.  Librational  stability  refers  to  the  stability  of  the  oscillations  of  the  moon 
about  its  center  of  mass  as  it  circles  the  Earth.  Galileo  was  the  hrst  to  notice  these 
oscillations.  Newton  conjectured  that  a  reason  for  this  behavior  would  be  that  the 
moon  was  elongated  towards  the  Earth.  However,  it  took  Louis  Lagrange  to  develop  a 
mathematical  theory  describing  this  phenomenon  [Lagrange].  Using  a  linear  analysis, 
Lagrange  derived  a  set  of  four  inequalities  involving  the  inertia  of  the  moon  that  must 
be  satished  for  the  moon  to  exhibit  stable  librational  motion. 

In  1885,  Henri  Poincare  [Poincare]  realized  that  a  linearized  analysis  could  not 
be  conclusive  in  determining  librational  stability.  This  result  gave  rise  to  the  use  of 
the  Hamiltonian  as  a  Lyapunov  function  candidate  in  determining  Lyapunov  stability. 
This  approach  established  the  stability  of  a  satellite  conhguration  called  the  Lagrange 
satellite  but  it  could  not  establish  the  stability  of  the  Delp  satellite  which  was  also 
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stable  according  to  linear  analysis. 

The  only  other  example  that  the  author  has  found  of  stable  rotation  about  the 
intermediate  axis  is  from  a  study  of  large  cellestial  structures.  Duncan  and  Levison 
[Duncan  89]  simulated  the  behavior  of  a  self-gravitating  system  of  2048  bodies  in 
order  to  determine  if  it  was  stable.  They  found  an  example  of  a  simple  spherical 
system  initially  in  dynamic  equilibrium  that  experienced  an  instability  producing  a 
hnal  equilibrium  state  of  stable  rotation  around  the  intermediate  principal  axis  of 
inertia  of  the  system  of  particles.  This  result  was  considered  noteworthy  because  it 
conflicted  with  the  rigid  body  analogy  of  unstable  rotation  about  the  middle  principal 
axis.  No  explanation  of  the  source  of  stability  was  offered. 

1.1.3  Ballistic  Walking 

While  the  signihcance  of  passive  dynamic  stability  is  recognized  in  studies  of  celestial 
structures  and  in  the  design  of  aircraft  and  spacecraft  its  relevance  to  the  control  of 
movement  in  biology  and  robotics  is  only  beginning  to  be  explored.  One  possible 
reason  is  that  these  former  examples  typically  involve  the  stability  of  an  equilibrium 
conhguration,  i.e.  no  accelerations.  Animals,  people,  and  robots  frequently  move 
with  signihcant  accelerations.  The  stability  analysis  of  non-equilibrium  movement  is 
a  much  more  difficult  process.  Some  progress  has  been  made  with  the  analysis  of 
walking. 

Mochon  and  McMahon  and  later  McGeer  showed  that  passive  dynamic  stability 
may  be  important  to  human  locomotion.  First  proposed  by  Mochon  and  McMa¬ 
hon  [Mochon  80],  a  ballistic  walker  uses  only  gravity  and  the  dynamic  interaction 
of  the  swing  and  stance  legs  to  produce  a  repetitive  walking  pattern.  The  passive 
pattern  accounted  for  the  folding  and  unfolding  of  the  legs  and  the  positioning  of 
the  foot  forward.  McGeer  [McGeer  89]  showed  the  viability  of  ballistic  walking  by 
building  passive,  planar,  anthropomorphic  linkages  with  no  sensor  or  actuators  that 
demonstrated  stable  walking  down  an  incline. 
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1.2  Organization  of  Thesis 

The  remainder  of  this  thesis  is  organized  as  follows. 

Ch.  2  This  chapter  provides  a  review  of  rigid  body  rotation.  Rigid  body  rotation 
provides  an  important  simplihed  model  for  the  analysis  of  multi-body  rotational 
systems.  Analysis  and  tools  introduced  here  will  be  used  throughout  the  thesis. 

Ch.  3  In  this  chapter  I  discuss  somersault  experiments  with  a  3D  biped  robot.  The 
robot  somersault  axis  is  coincident  with  the  major  principal  axis  of  inertia  so  the 
maneuver  exhibits  some  passive  stability  properties.  The  robot  actively  controls 
landing  attitude  by  retracting  or  extending  its  legs  during  the  maneuver  to 
change  somersault  rate.  The  robot  has  successfully  performed  front  somersaults 
in  the  laboratory. 

Ch.  4  In  this  chapter  I  present  an  analysis  of  the  layout  somersault.  I  show  that  the 
layout  somersault,  involving  rotation  about  the  middle  principal  axis  of  inertia, 
can  be  passively  stabilized  by  tuning  parameters  of  a  passive  dynamic  model  of 
the  human  body.  Using  the  simplest  possible  model  of  the  layout  somersault,  I 
explain  the  fundamental  dynamics  of  passive  stabilization. 

Ch.  5  I  describe  layout  somersault  experiments  with  a  human-like  doll.  These  exper¬ 
iments  verify  that  the  layout  somersault  can  be  consistently  stabilized  for  at 
least  three  and  one  half  somersaults. 

Ch.  6  I  discuss  non-linear  dynamic  simulation  of  a  1  1/2  twisting  front  somersault. 
Simulation  results  suggest  that  maneuvers  using  prescribed  limb  motion  will 
not  be  reliable  to  off-nominal  initial  conditions.  If  instead,  a  compliant  passive 
model  is  used  with  feed  forward  torques  the  maneuver  can  be  made  more  reli¬ 
able.  I  also  describe  twisting  somersault  experiments  with  the  3D  Biped  robot. 
Using  simulation  we  found  that  it  was  important  that  the  robot  have  an  inertia 
tensor  more  like  that  of  a  human  to  produce  a  human-like  front  somersault  with 
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1/2  twist.  Laboratory  experiments  with  3D  Biped  robot  twisting  somersaults 
were  not  successful  at  least  in  part  due  to  insufficient  actuator  power. 

Ch.  7  Here  I  summarize  the  results  of  the  thesis  and  discuss  future  work. 

Appx.  In  the  Appendix  I  provide  Mathematica  code  for  deriving  non-linear  equations 
of  motion  of  a  human  model,  and  for  analytically  linearizing  this  model.  I 
provide  a  dehnition  of  the  parameters  used  in  the  linear  model.  I  also  provide 
a  derivation  of  the  analytic  solution  for  rigid  body  rotation. 
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Rigid  Body  Rotation 

2.1  Introduction 

Humans  are  multi-body  systems  and  gymnastic  maneuvers  involve  multi-body  rota¬ 
tions.  We  can  better  understand  multi-body  rotation  by  understanding  the  dynamics 
and  stability  of  rigid  body  rotation.  In  this  chapter,  I  briefly  discuss  some  properties 
of  rigid  body  rotation.  Rigid  body  rotational  stability  about  a  principal  axis  de¬ 
pends  only  on  the  relative  magnitude  of  the  principal  inertias.  This  result  provides  a 
valuable  reference  for  multi-body  rotational  stability.  A  linear  analysis  of  rigid  body 
rotation  provides  simple  stability  results,  and  identihes  important  non-dimensional 
parameters  that  will  be  useful  in  multi-body  analysis.  Visualization  tools  for  rigid 
body  rotation  will  also  prove  useful  in  understanding  how  to  control  rotational  motion 
in  multi-body  systems.  Material  for  this  chapter  is  based  on  the  text  of  [Hughes  86]. 

2.2  Nonlinear  Dynamic  Model 

The  rotational  motion  of  a  free  rigid  body  about  its  center  of  mass  is  governed  by 
Euler’s  equations. 


27 
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IlUJl  —  {I2  —  l3)L02i03  =  0 

I2O2  —  {I3  —  Ii)‘-^3‘^i  =  0  (2-1) 

I3O3  —  (/l  —  l2)>-^l>-^2  =  0 


where  Ij  refers  to  the  principal  moment  of  inertia  and  ujj  refers  to  the  angular 
velocity  about  that  principal  axis. 

The  orientation  of  the  rigid  body  with  respect  to  an  inertial  coordinate  frame  is 
described  by  a  set  of  three  Euler  angles.  The  Euler  angles  are  used  to  dehne  the 
3-by-3  matrix  of  direction  cosines,  Cbi^  that  transforms  a  vector  described  in  inertial 
coordinates  into  one  described  in  the  body  hxed  coordinate  system.  (Eor  the  rigid 
body  analyses  of  this  chapter  I  assume  that  the  principal  axis  frame  is  coincident 
with  this  body  axis  frame.)  I  use  a  ‘2-1-3’  Euler  angle  sequence  for  this  purpose. 
I  borrow  the  names  for  the  three  Euler  angles,  somersault,  tilt,  and  twist,  from 
Yeadon  [Yeadon  84].  I  use  the  letters,  $,0,T  to  refer  to  the  somersault,  tilt,  and 
twist  angles  respectively. 

To  describe  the  attitude  of  the  body  in  inertial  space,  a  coordinate  system  initially 
parallel  to  the  inertial  reference  frame  is  hrst  rotated  through  the  somersault  angle 
about  the  inertially  hxed  ’2’  axis,  then  rotated  through  the  tilt  angle  about  the 
intermediate  ’1’  body  axis,  and  hnally  rotated  through  the  twist  angle  about  the 
body  hxed  ’3’  axis.  Eigure  2-1  depicts  the  $,  0,  T  Euler  angle  sequence  applied  to  a 
rigid  human  form.  The  rotation  matrix  is  given  by 


Cbi  = 


C\siC%  +  SisiSsS^  SisiCs  —CisiS^  +  S\siSsC% 

—S^CA  +  CASqS^  C^iCq  SASA  +  CASqCA 

U0O<|)  U0U<|) 


(2.2) 


The  subscript  bi  refers  to  the  fact  that  this  rotation  matrix  rotates  a  vector  from  the 
mertial  system  into  the  body  axis  system.  S  and  C  are  the  sine  and  cosine  of  the 
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Figure  2-1:  Illustration  of  the  somersault  ($),  tilt  (0),  and  twist  (d^)  Euler  angle 
sequence  used  to  define  body  attitude. 


respective  angle.  This  Euler  angle  description  of  body  attitude  has  a  singularity,  as  do 
all  three  parameter  descriptions  of  body  attitude.  The  singularity  for  this  particular 
sequence  occurs  at  a  tilt  angle  of  ±7r/2.  The  kinematic  equations  governing  the 
evolution  of  the  Euler  angles  are  given  by 


T 

S^SqCq  CASqCq  1 

UJl 

0 

= 

o 

1 

iL>2 

$ 

S^iCq  C^iCq  0 

ids 

The  inverse  of  the  above  relationship  is  also  useful.  It  is  given  by 


UJl 

0  C\^  S\^Cq 

T 

UJ2 

= 

0  —5*1])  C\^Cq 

0 

LJo, 

1 

o 

1 

CD 

$ 

(2.3) 


(2.4) 
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The  dynamic  equations  of  motion  (2.1)  are  coupled  so  that  general  rotational 
motion  involves  all  three  degrees  of  freedom.  However,  three  simple  solutions  exist  to 
these  equations.  If  a  rigid  body  is  spinning  perfectly  around  a  principal  axis  such  that 
two  of  the  three  angular  velocities  are  zero  then  the  non-linear  coupling  terms  vanish. 
The  body  will  continue  to  spin  about  that  axis  without  coupling  to  the  other  degrees 
of  freedom.  An  analysis  of  these  spin  solutions  will  reveal  that  for  a  tri-inertial  body, 
a  body  with  three  different  principal  inertias,  only  two  of  the  spin  solutions  are  stable 
while  the  third  is  unstable.  Here  stability  means  if  the  spin  axis  of  the  system  is 
moved  away  from  the  nominal  solution  the  attitude  of  the  spin  axis  will  not  diverge 
radically  from  its  initial  orientation.  The  linear  stability  analysis  provided  in  the  next 
section  reveals  the  required  conditions  for  stable  rotation  about  a  principal  axis. 

2.3  Linear  Analysis 

In  this  section  we  linearize  Euler’s  equations  for  the  rigid  body  in  order  to  study  the 
stability  of  steady  rotation  about  a  principal  axis.  The  linear  equations  govern  the 
motion  of  the  body  relative  to  a  reference  frame  that  is  rotating  steadily  about  the 
principal  axis  of  inertia.  The  resulting  equations  will  be  valid  for  small  deviations  of 
the  system  from  the  reference  frame.  This  linearization  process  will  be  used  again 
when  we  linearize  multi-body  rotation  about  an  equilibrium  spin  condition.  The 
linearizing  condition  will  be  steady  somersaulting  rotation  about  the  ‘2‘  axis  with 
angular  velocity  $.  The  angular  velocity  of  the  steadily  rotating  reference  coordinate 
system  is  ujra- 

iO^a  —  (^*^) 

where  A2  is  the  unit  vector  along  the  ’2’  axis  in  the  reference  coordinate  system.  The 
angular  velocity  of  the  principal  axis  system,  is  comprised  of  the  angular  velocity 


CHAPTER  2.  RIGID  BODY  ROTATION 


31 


relative  to  the  reference  coordinate  frame,  and 


iOpa  iOpf  Y  ^ra 


(2.6) 


Let  the  attitude  of  the  principal  axis  system  relative  to  the  reference  frame  be  de¬ 
scribed  by  the  2-1-3  Euler  angle  sequence  of  ^,0,?/’.  (The  lower  case  notation  indi¬ 
cates  that  these  are  linearized  states.  These  Euler  angles  describe  the  deviations  of 
the  body  from  the  rotating  reference  frame.)  We  can  use  Equations  2.2  and  2.4  to 
express  the  components  of  Copa  in  terms  of  0,  if),  their  derivatives  and  $.  Eirst  using 
Equation  2.4 


0  C-y  h'yC'e 

0  -h'y  C^Ce 

e 

1  0  -Se 

(2.7) 


where  Upr  denotes  a  column  whose  components  are  the  elements  of  the  vector  Apr 
when  described  in  the  principal  axis  frame.  To  form  the  sum  in  Equation  2.6,  we 
need  to  express  all  the  vector  components  in  the  same  frame.  To  do  this,  using 
Equation  2.2  as  a  model,  form  the  rotation  matrix  Cpr  that  rotates  a  vector  from  the 
reference  frame  into  the  principal  coordinate  frame.  Then  solve  for 


iOpr  T  Cp 


When  Cpr  and  Upr  are  linearized  for  small  V’?  ^  we  get 


0 

p 

^pa  — 

1 

-9 

(2.8) 
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This  expression  can  be  differentiated  to  find  the  Euler  angle  expression  for  angular 
acceleration. 

9  ij 

Opa  =  ^  $  0  (2-9) 

^  -9 

Now,  substituting  Equations  2.8  and  2.9  into  Equation  2.1  and  eliminating  terms  that 
are  second  order  in  0,  and  tj)  results  in  the  following  equations 

Ij-^{l2-l3-ll)i^-^^{l3-l2)9  =  0 

=  0  (2.10) 

/3^  +  $(/2_/3_/i)0_$2(Ji_  J2)^  =  0 

These  equations  govern  the  behavior  of  the  non-linear  system  in  the  vicinity  of  the 
pure  spin  solution.  The  linear  states,  9,  (f),  tf),  describe  the  deviation  of  the  body  axes 
from  the  reference  coordinate  frame. 

Examination  of  Equations  2.10  reveals  that  the  system  is  unstable  in  the  sense 
that  perturbations  in  ^  will  result  in  unbounded  growth  in  cf).  Nevertheless,  a  reduced 
form  of  stability  called  directional  stability  [Hughes  86]  is  possible  in  the  subsystem 
of  0,V’-  Directional  stability  means  that  the  two  attitude  variables  9  and  tj)  will  not 
diverge  from  zero  if  the  system  is  perturbed  slightly  from  the  pure  spin  solutions. 
This  means  that  the  spin  axis  will  continue  to  point  in  roughly  the  same  inertial 
direction  when  disturbed  from  its  equilibrium  orientation. 

The  equations  for  9  and  tj)  decouple  from  cf).  They  can  be  written  as 

Mx  +  ^Gx  +  ^^Kx  =  0  (2-11) 

where 

9 

X  = 

.  ^  . 
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M 


h  0 

0  /s 


G 


0  -(/2-/3-/1) 

I2  —  I3  —  G  0 


K 


I2  —  I3  0 
0  I2  —  h 


2.3.1  Conservative  Gyric  Systems 


To  investigate  the  stability  of  the  9  —  if)  subsystem,  I  use  the  matrix  second  order 
stability  theory  described  in  [Hughes  86].  Hughes  describes  a  second  order  system  of 
equations  of  the  form  of  Equation  2.11  as  a  conservative  gyric  system  when 


=  M  >0 
G^  =  -G 
=  K 

The  coefficient  matrices  (Af,  G  and  K)  are  respectively  associated  with  inertial,  gyric, 
and  stiffness  forces.  Conservative  refers  to  the  fact  that  the  system  energy  is  con¬ 
served,  while  gyric  reflects  the  fact  that  these  terms  often  arise  in  spinning  systems. 
This  form  of  equations  will  be  present  in  a  multi-body  analysis  of  rotating  systems. 

Hughes  proves  that  asymptotic  stability  for  a  conservative  gyric  system  is  impos¬ 
sible  by  showing  that  if  s  is  a  root  of  the  characteristic  equation  of  2.11  then  —  s  must 
also  be  a  root.  Strictly  left  half  plane  poles  will  always  have  their  right  half  plane 
counterparts.  Therefore,  stability,  as  opposed  to  asymptotic  stability,  is  the  strongest 
result  possible  for  a  conservative  gyric  system.  A  stable  system  will  have  all  of  the 
roots  of  its  characteristic  equation  on  the  imaginary  axis. 

A  sufficient  condition  for  stability  of  a  conservative  gyric  system  is  that  it  be 
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statically  stable,  i.e.  the  stiffness  matrix  must  be  positive  definite,  K  >  0.  The 
stiffness  matrix  for  the  two-by-two  system  above  is  positive  dehnite  if  R  is  the  major 
principal  axis  of  inertia.  That  is,  rotation  about  the  major  axis  is  directionally  stable 
by  virtue  of  its  static  stability.  However,  K  >  0  is  only  a  sufficient  condition  for 
stability.  Rotation  about  the  minor  principal  axis  is  also  stable.  Since  in  this  case  the 
stiffness  matrix  is  not  positive  dehnite  the  system  is  considered  gyrically  stabilized. 
To  test  for  gyric  stability  we  can  check  the  roots  of  the  characteristic  equation  of  the 
second  order  system.  The  system  will  be  considered  stable  if  the  roots  are  purely 
imaginary. 

These  results  do  not  preclude  the  possibility  that  asymptotic  stability  could  be 
obtained  by  adding  damping  to  a  conservative  gyric  system.  However,  while  damp¬ 
ing  tends  to  make  statically  stable  systems  become  asymptotically  stable  systems, 
damping  also  tends  to  destabilize  gyrically  stabilized  systems.  Hughes  proves  that 
statically  unstable  systems  are  destabilized  if  they  have  a  positive  dehnite  damping 
matrix. 


2.3.2  Rigid  Body  Inertia  Ratios 

Before  we  solve  for  the  roots  of  the  characteristic  equation  of  Equation  2.11  we  should 
note  that  we  can  simplify  our  analysis  by  recognizing  that  only  the  ratios  of  inertia 
are  important  to  the  analysis  rather  than  the  individual  values  of  inertia.  Dividing 
each  equation  of  Equation  2.11  by  the  corresponding  diagonal  term  of  M  results  in 
the  following  form  of  Af,  G  and  K . 


M 


1  0 
0  1 


G 


0  1-h 

(1  -  h)  0 
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K 


ki  0 
0  ks 


where 


A  physical 
constraints 


,  _  h  —  h  ,  _  h  —  h 

M  —  - } - 5  j - 

body  can  not  have  arbitrary  values  of  Ii^R  and  I3  [Hughes  86]. 
are  captured  by  the  fact  that 


These 


-1  <  <  1 


-1  <  A;3  <  1 

Therefore,  all  possible  rigid  body  conhgurations  can  be  represented  on  a  plot  of  the 
parameter  space  of  ki  and  ks  restricted  to  the  unit  square.  Figure  2-2  shows  how  ki 
and  ks  depend  upon  the  relative  size  of  Ii,  I2,  and  I3.  Also  included  in  the  hgure  are 
schematic  drawings  of  rectangular  prisms  that  would  have  approximately  the  correct 
inertia  ratios  for  selected  points  around  the  diagram. 

The  characteristic  equation  of  this  system  is  formed  with  the  following  determinant 

det[Ms‘^  +  $(75  +  A"]  =  0 

which  results  in  the  following 

(5) 


bi  —  1  kik3 

62  =  ki  k3 


where 
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Figure  2-2:  The  inertia  ratio  diagram  shows  how  ki  and  ks  depend  upon  the  relative 
value  of  the  principal  inertias.  ki  =  ,  k^  =  ■  The  different  rectangular 

prisms  located  near  the  axes  and  corners  of  the  diagram  indicate  an  example  shape 
that  would  correspond  to  the  local  inertia  ratios. 
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For  stability  we  require 

1.  bi  =  1  +  kiks  >  0 

2.  1)2  =  kiks  >  0 

3.  (Gf  -  ib2  =  (1  -  k^ksf  >0 

The  only  stability  condition  not  met  automatically  is  kiks  >  0.  This  requirement 
specihes  that  ki  and  ks  must  have  the  same  sign.  On  a  plot  of  the  {A;l,  k3}  parameter 
space  stability  occurs  in  the  hrst  and  third  quadrants. 

Figure  2-3  shows  the  regions  of  stability  for  the  simple  spin  solutions  of  a  rigid 
body. 

Fquation  2.12  is  simple  enough  that  we  can  solve  for  the  roots  of  the  characteristic 
equation  in  closed  form.  Two  roots  are  located  at  s  =  This  is  the  stroboscopic 

mode  of  rotation.  This  mode  of  rotation  occurs  when  the  body  is  still  spinning 
perfectly  about  the  principal  axis  but  this  axis  is  offset  from  the  original  orientation. 
The  other  two  roots  are  located  at  s  =  ±z/,  where 

//  =  [k^ksY^^^ 

In  the  hrst  and  third  quadrants  of  the  inertia  ratio  diagram,  p  is  positive  and  therefore 
these  roots  are  purely  oscillatory.  In  the  second  and  fourth  quadrants  of  Figure  2-2, 
p  is  imaginary  forcing  the  roots  of  the  characteristic  equation  to  have  real  positive 
and  negative  values.  In  these  quadrants,  the  unstable  mode  of  motion  is  governed  by 
the  following  equation 


X  =  where  s  =  {  —  kiksY^^^t 


(2.13) 


To  get  an  idea  of  how  unstable  the  system  is  we  compute  the  change  in  the  nominal 
somersault  angle,  $,  required  for  the  the  unstable  rotational  mode  to  grow  by  a  factor 
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Unstable 


Directionally  and  Statically  Stable 


Directionally  Stable,  Statically  Unstable 
but  Gyrically  Stabilized 


Figure  2-3:  Stability  diagram  for  the  simple  spin  solutions  of  a  rigid  body.  Rigid 
body  rotational  stability  depends  only  upon  the  two  non-dimensional  inertia  ratios 
ki  and  ks.  We  assume  the  body  is  rotating  about  the  principal  axis  corresponding  to 
I2  for  this  diagram. 
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Figure  2-4:  All  rigid  bodies  with  inertia  ratios  on  a  single  curve  have  an  unstable 
mode  that  grows  at  the  same  rate.  The  value  of  a  curve,  indicated  in  the  hgure,  is 
the  somersault  angle  (in  radians)  the  body  must  execute  before  the  unstable  twist 
mode  grows  by  a  factor  of  ten.  The  most  unstable  systems  are  those  with  inertia 
ratios  in  the  corners  {ki^ks}  =  {1,-1}  or  {  —  1,1}.  This  plot  is  symmetric  about  the 
origin. 


of  N. 


ln{N) 

{-kiksYD 


Figure  2-4  shows  curves  of  constant  $  on  the  {ki^  ks}  axes  for  N  =  10. 


(2.14) 


2.4  Rotational  Stability  of  the  Rigid  Human  Body 

For  any  particular  conhguration  of  the  human  body,  we  can  solve  for  the  orientation 
and  magnitude  of  the  principal  axes  of  inertia.  This  allows  us  to  compute  the  cor¬ 
responding  inertia  ratios  which  in  turn  provide  a  valuable  reference  for  the  stability 
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4.  L 


Figure  2-5:  Tuck,  pike,  and  layout  positions  for  a  human. 

of  the  human  body  in  rotation  about  the  principal  axes.  Figure  2-5  shows  from  left 
to  right  a  human  in  the  tuck,  pike  and  layout  positions.  Figure  2-6  shows  the  inertia 
ratios  of  a  human  performer  moving  from  a  tuck  position  through  a  pike  position  to 
a  layout  position.  To  make  this  hgure,  I  assumed  that  R  was  the  principal  inertia 
along  the  somersault  axis  of  the  human.  This  analysis  shows  that  for  a  rigid  body, 
tuck  and  pike  somersault  conhgurations  are  stable  in  rotation  about  the  somersault 
axis  while  steady  rotation  in  the  layout  position  is  unstable.  In  the  pike  and  tuck 
positions  the  inertia  ratios  are  in  the  stable  upper  right  quadrant  of  the  inertia  ratio 
diagram  while  the  inertia  ratios  of  the  layout  somersault  are  in  the  unstable  lower 
right  quadrant. 


Figure  2-6:  Inertia  ratios  of  a  human  performer  for  a  sequence  of  configurations 
connecting  tuck,  pike  and  layout  positions.  The  tuck  and  pike  positions  are  stable 
(for  a  rigid  body)  while  the  layout  is  unstable. 
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2.5  A  Map  of  Rigid  Body  Rotation 

The  linear  analysis  we  have  used  so  far  is  restricted  to  three  solutions  for  rigid  body 
rotation,  pure  spin  about  each  of  the  principal  axes.  A  general  solution  for  the  torque- 
free  motion  of  a  rigid  body  tumbling  in  space  would  be  useful  for  understanding  the 
range  of  behavior  a  rigid  body  can  exhibit  between  these  three  special  solutions. 
Solving  for  such  a  solution  is  one  of  the  classic  problems  of  dynamics.  The  con¬ 
tribution  that  this  solution  offers  today  is  a  concise  description  of  the  states  and 
non-dimensional  parameters  that  govern  rigid  body  rotation  (see  Appendix  A. 4).  In 
addition,  the  analytic  solutions  give  rise  to  elegant  geometric  interpretations  of  rigid 
body  rotation  that  help  to  provide  intuition  of  this  otherwise  complex  movement.  In 
this  section,  I  present  one  form  of  geometric  tool  derived  from  the  analytic  solution 
that  is  useful  in  visualizing  rigid  body  rotation.  I  think  of  this  tool  as  a  map  of  rigid 
body  rotation  because  it  shows  graphically  two  of  the  three  Euler  angle  trajectories 
involved  in  non-linear  rotational  motion  of  a  rigid  body.  This  map  not  only  makes 
clear  the  stable  and  unstable  axes  of  rotation  but  also  shows  two  distinct  regions 
of  qualitatively  different  motion.  A  performer  can  exert  control  over  his  rotational 
motion  by  moving  within  and  inbetween  these  two  regions. 

Two  integrals  of  motion  are  used  to  dehne  a  map  of  rigid  body  rotation.  During 
flight,  a  rotating  rigid  body  must  conserve  angular  momentum,  h,  and  kinetic  energy, 
T.  The  magnitude  of  both  of  these  quantities  can  be  written  as  the  equation  for  an 
ellipsoid  in  angular  velocity  space 

=  I^iol  +  I^t^l  +  /3CU3  (2.15) 

2T  =  -^3^3-  (2.16) 

Assume,  with  no  loss  of  generality,  that  the  angular  momentum  vector  is  aligned  with 
the  inertially  hxed  ’2’  axis.  Using  the  direction  cosine  matrix,  (2.2),  we  can  solve  for 
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the  components  of  h  in  the  principal  coordinate  frame. 


h  —  hciipi  +  hc2ip2  +  hc^ips 


where  {ci,  C2,  C3}  are  the  direction  cosines  of  h  and  ipj  is  the  unit  vector  in  the  prin¬ 
cipal  axis  system.  Alternatively,  the  principal  axes  components  of  angular  momentum 
can  also  be  written  as  follows 


h  —  IiLOiipi  +  /2^2*p2  +  -^3^3*p3 


Thus  we  have  established  that 


hci  =  IiUi,  hc2  =  121^2,  hcs  =  I3UJ3, 

The  ellipsoid  describing  angular  momentum  (2.15)  can  now  be  written  as  the  equation 
of  a  sphere  in  direction  cosine  space. 

2,2,2  1 

Cl  +  C2  +  C3  —  1 


Similarly  the  kinetic  energy  can  be  written  as  the  equation  of  an  ellipsoid. 


c?  c?  c? 

_i  _2  _3 

/l  /2  h 


2T 


For  a  hxed  kinetic  energy  and  angular  momentum,  the  direction  cosines  must  si¬ 
multaneously  lie  on  the  surface  of  both  the  momentum  sphere  and  energy  ellipsoid. 
Therefore,  the  intersection  of  these  two  shapes  describes  a  trajectory  in  direction 
cosine  space  that  the  body  must  ‘follow’. 

Figure  2-7  shows  a  sample  map  for  the  possible  rotational  trajectories  of  a  ‘rigid’ 
human  body.  Fach  trajectory  shown  on  the  sphere  corresponds  to  a  different  kinetic 
energy  of  rotation.  The  axes  of  the  sphere  are  the  direction  cosine  axes.  The  sphere 
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Figure  2-7:  A  map  of  rigid  body  rotation  for  a  human  in  the  layout  position.  The 
axes  of  the  sphere  and  the  principal  axes  of  the  body  remain  parallel  as  the  body 
rotates.  The  inertially  hxed  angular  momentum  vector  (black)  paints  trajectories 
onto  the  surface  of  the  sphere  as  the  sphere  rotates.  Each  trajectory  on  the  sphere 
corresponds  to  a  different  rotational  energy.  The  trajectories  indicate  the  tilt  and 
twist  angles  of  the  body  as  it  rotates.  This  map  does  not  include  somersault  angle 
and  does  not  show  the  time  dependence  of  the  tilt  and  twist  Euler  angles. 

rotates  with  the  body  so  that  its  axes  remain  parallel  to  the  principal  axes  of  the 
body.  The  angular  momentum  vector  is  shown  protruding  from  the  sphere.  The 
sphere  (and  body)  must  move  so  that  the  (inertially  hxed)  angular  momentum  vector 
remains  in  the  ’slot’  that  is  appropriate  for  the  given  kinetic  energy. 

The  stable  and  unstable  axes  of  rotation  are  obvious  from  this  map.  The  stable 
principal  axes  are  surrounded  by  trajectories  that  enclose  the  axes  while  the  unstable 
middle  principal  axis  shows  trajectories  that  converge  then  diverge  from  that  axis. 
This  map  also  shows  that  the  trajectories  are  divided  into  two  qualitatively  different 
regions.  This  division  is  based  on  the  rotational  kinetic  energy. 
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For  a  fixed  angular  momentum,  the  energy  of  rotation  is  bounded  above  and  below 
by  the  rotational  energy  in  pure  spin  about  the  minimum  and  maximum  principal 
axes  respectively. 
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Between  these  extremes  exists  a  continuum  of  energy  levels  of  rotational  motion. 
The  dividing  point  among  these  trajectories  is  the  energy  required  to  spin  about  the 
middle  principal  axis,  r  =  7^ — .  Those  trajectories  with  higher  energy  involve  monot- 
ically  increasing  (or  decreasing)  twist  while  those  with  less  energy  involve  oscillatory 
twisting.  In  Figure  2-7,  the  trajectories  which  enclose  the  minimum  principal  axis 
(head-to-toe  axis)  of  the  performer  involve  monotonic  twist  while  those  that  enclose 
the  maximum  principal  axis  of  inertia  involve  twist  angles  that  oscillate  between  0 
and  ±7r  depending  on  which  region  the  trajectory  is  located  in. 

We  can  read  the  twist  and  tilt  Fuler  angles  of  a  trajectory  from  this  map  in  an 
intuitive  way.  Using  Fquation  2.2  we  can  derive  an  explicit  relationship  between  the 
direction  cosines  and  the  Fuler  angles  as  follows 


Cl 

S^Ce 

C2 

= 

C^Ce 

C3 

-Se 

Drawing  the  components  of  the  angular  momentum  vector  as  shown  in  Figure  2-8,  we 
realize  that  the  spherical  coordinates  of  the  unit  angular  momentum  vector  are  defined 
by  the  tilt  and  twist  Fuler  angles,  {0,  T}.  The  tilt  angle,  0  is  the  (negative)  latitude 
of  the  energy  curve  and  the  twist  angle,  T,  is  the  (negative)  longitude.  Therefore  we 
can  simply  read  the  tilt  and  twist  Fuler  angles  of  the  body  from  the  polar  coordinates 
of  the  angular  momentum  vector  on  this  map.  For  example,  the  attitude  of  the  body 
shown  in  Figure  2-7  is  approximately  0  =  3°,  T  =  0. 
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Figure  2-8:  The  spherical  coordinates  used  to  dehne  orientation  of  the  unit  angular 
momentum  vector  in  direction  cosine  space  are  the  negative  of  the  tilt  and  twist  Euler 
angles. 

2.6  Dynamic  Simulation  Environment 

The  analytic  and  experimental  results  presented  in  this  thesis  are  complemented  with 
numeric  results  from  non-linear  dynamic  simulation  of  the  creature  in  question.  This 
section  briefly  describes  the  simulation  environment  used  to  compute  the  motion  and 
produce  graphic  output  from  dynamic  simulations. 

The  simulation  environment  consists  of  three  parts  that  allow  us  to  compute  the 
movement  of  a  creature,  make  a  movie  of  its  motion  using  computer  graphics,  and 
analyze  its  motion  using  time  histories  of  simulated  data.  Using  software  developed 
in  this  laboratory,  the  Leg  Laboratory  at  MIT,  we  have  integrated  the  simulation  en¬ 
vironment  so  that  the  three  parts  are  all  generated  automatically  from  a  single  input 
hie  that  describes  the  shape,  topology,  and  mass  properties  of  the  system.  The  part 
that  computes  the  motion  does  so  by  numerically  integrating  the  non-linear  equations 
of  motion  for  a  given  creature.  The  non-linear  equations  of  motion  are  derived  us¬ 
ing  the  commercially  available  software  package,  SD-FAST.  This  package  derives  the 
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equations  of  motion  for  a  rigid  link,  multi-body  system  using  the  description  file  as 
input.  The  equations  of  motion  are  automatically  incorporated  into  a  dynamic  simu¬ 
lation  that  handles  user  interface  to  the  software.  The  computer  graphic  and  analysis 
software  is  also  generated  automatically  from  the  description  of  the  creature  allow¬ 
ing  us  to  animate  movement  from  simulation  data  using  either  a  simple  “working” 
picture  or  a  higher  quality,  computer  graphic  image.  The  only  part  of  the  simulation 
environment  that  is  not  automatically  generated  is  the  control  software.  This  soft¬ 
ware  is  used  to  specify  the  desired  behavior  of  the  system.  It  does  so  by  computing 
joint  torques  or  other  actuator  inputs  that  we  assume  the  creature  to  have.  Passive 
forces  like  those  due  to  springs,  dampers  and  gravity  are  input  to  the  equations  of 
motion  in  a  manner  similar  to  active  joint  torques.  The  ability  to  quickly  produce 
dynamic  simulations  of  a  variety  of  creatures  has  made  this  simulation  software  a 
useful  research  tool  in  the  investigation  of  machine,  human  or  animal  movement. 

2.7  Summary 

Material  in  this  chapter  is  derived  from  the  text  of  [Hughes  86].  In  this  chapter,  I 
reviewed  the  dynamic  and  kinematic  equations  of  rigid  body  rotation.  Linear  equa¬ 
tions  of  motion  describe  rigid  body  rotation  in  the  neighborhood  of  the  equilibrium 
solution  of  pure  spin  about  a  principal  axis.  A  stability  analysis  of  these  equations 
shows  that  rigid  body  rotation  is  stable  about  the  major  and  minor  principal  axes  but 
is  unstable  about  the  intermediate  axis.  The  linear  stability  analysis  can  be  simpli- 
hed  by  the  use  of  the  non-dimensional  inertia  ratios,  ki  and  ks.  A  stability  diagram 
simply  demarcates  stable  and  unstable  regions  in  the  ki^  k^  parameter  space.  I  show 
the  values  of  ki  and  k^  that  a  (rigid)  human  body  would  assume  in  the  tuck,  pike, 
and  layout  positions  in  order  to  provide  a  stability  reference  for  the  rotating  human 
body.  A  geometric  map  derived  from  the  closed-form  solution  of  rigid  body  rotation 
captures  the  range  of  possible  rotational  modes.  This  map  intuitively  shows  the  stable 
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and  unstable  axes  of  rotation.  It  also  demarcates  between  two  regions  of  qualitatively 
different  rotational  motion. 


Chapter  3 


Robot  Tucked  Somersaults 

3.1  Introduction 

The  somersault  is  a  maneuver  in  which  a  performer  jumps  into  the  air  and  rotates 
once  about  a  side-to-side  axis  before  landing  on  the  ground.  The  main  requirement  of 
a  successful  somersault  is  a  balanced  landing  which  in  turn  requires  that  the  performer 
hnish  the  somersault  with  a  pre-specihed  body  attitude.  In  this  chapter,  I  explore 
attitude  control  techniques  for  producing  stable  landing  conhgurations  for  the  tucked 
somersault  of  a  3D  biped  running  robot. 

The  tucked  somersault  is  distinguished  by  maintaining  the  tuck  position  during 
most  of  the  maneuver.  Humans  ‘tuck’  by  holding  their  knees  close  to  the  chest 
with  the  knee  joints  flexed  to  fold  the  lower  legs  under  the  body.  An  important 
dynamic  feature  of  the  tuck  somersault  is  that  in  most  humans  it  involves  rotation 
about  the  major  principal  axis  of  inertia.  Since  a  rigid  body  rotating  about  its  major 
principal  axis  is  directionally  stable,  we  might  expect  the  tuck  somersault  to  be  stable 
in  the  sense  that  the  axis  of  rotation  will  tend  to  maintain  its  inertial  orientation. 
This  stability  in  turn  simplihes  the  control  of  body  attitude  at  landing  in  the  tuck 
somersault. 

Directional  stability  of  the  spin  axis  during  a  tucked  somersault  means  that  a 
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Figure  3-1:  Photograph  of  the  3D  Biped  used  for  experiments.  The  body  is  an 
aluminum  frame  on  which  are  mounted  hip  actuators  and  computer  electronics.  Each 
ball  and  socket  hip  joint  has  three  low  friction  hydraulic  actuators  that  can  position 
the  leg  forward  and  aft,  side-to-side,  and  can  rotate  the  leg  along  the  axis  of  the  leg. 
A  hydraulic  actuator  within  each  leg  changes  its  length,  while  an  air  spring  makes 
the  leg  springy  in  the  axial  direction.  Sensors  measure  the  lengths  of  the  legs,  the 
positions  of  the  hip  actuators,  pressure  in  the  hip  actuators  and  contact  between 
the  foot  and  the  floor.  Gyroscopes  measure  the  inertial  attitude  of  the  body.  An 
umbilical  cable  connects  the  machine  to  hydraulic,  pneumatic,  and  electrical  power 
supplies.  Control  computations  are  done  my  microprocessors,  some  located  on  board 
and  some  nearby  in  the  laboratory.  Communications  cables  connect  all  processors. 

somersault  control  strategy  need  deal  primarily  with  avoiding  over-rotation  or  under¬ 
rotation  about  the  somersault  axis.  Rotation  rate  about  the  somersault  axis  can  be 
controlled  by  changing  rotational  inertia.  Since  angular  momentum  must  be  conserved 
during  flight,  increasing  inertia  will  slow  down  the  somersault  rate  while  decreasing 
inertia  will  increase  somersault  rate.  If  the  performer  knows  the  time  of  flight,  then 
control  of  somersault  rate  provides  a  means  of  controlling  the  somersault  angle  at 
landing. 

In  this  chapter,  I  discuss  somersault  experiments  with  a  3D  biped  robot,  (Figure  3- 
1).  Figure  3-2  shows  a  sequence  of  photographs  of  the  3D  Biped  taken  while  the  robot 
performed  a  successful  somersault.  The  somersault  axis  of  the  robot  is  coincident 
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with  the  maximum  principal  axis  of  inertia  making  the  robot  somersault  dynamically 
similar  to  the  tuck  somersault  in  humans.  To  initiate  the  somersault  the  biped  runs 
forward,  jumps  to  attain  a  double  stance  phase,  then  thrusts  with  both  legs  while 
pitching  its  body  forward.  Once  airborne,  the  robot  shortens  its  legs  (tucks)  to 
accelerate  the  forward  rotation  and  swings  its  legs  to  a  predetermined  position  with 
respect  to  the  body.  During  flight  the  robot  uses  a  feedback  algorithm  that  changes 
the  leg  length  to  produce  a  rotation  rate  that  will  yield  the  desired  somersault  angle 
at  landing.  To  accommodate  errors  in  the  estimated  time  of  landing,  the  robot 
moves  its  feet  to  track  the  desired  landing  conhguration  as  the  system  approaches  the 
ground.  The  robot  does  not  use  any  active  control  of  out-of-plane  rotation  during 
flight.  Rather  it  uses  a  broad  stance  during  takeoff  to  minimize  tilt  rotation  at  the 
beginning  of  flight  and  it  uses  a  broad  stance  during  landing  to  minimize  the  effect 
of  landing  tilt  errors.  The  robot  has  successfully  performed  the  somersault  in  the 
laboratory.  On  its  best  day,  the  robot  regained  balance  on  landing  to  continue  stable 
running  on  seven  out  of  ten  somersault  attempts. 

3.2  The  Mechanics  of  the  Somersault 

The  most  basic  requirement  of  a  somersault  is  that  the  performer  neither  over-rotate 
nor  under-rotate  the  landing.  Accurate  control  of  the  landing  attitude  allows  careful 
placement  of  the  foot  relative  to  the  center  of  mass  of  the  robot  which  is  a  requirement 
for  stable  dynamic  running  [Raibert  84].  Considering  only  the  somersault  degree  of 
freedom,  the  attitude  requirement  is  expressed  by  equating  the  time  of  flight  and  the 
time  to  rotate  through  the  desired  change  in  somersault  attitude, 

A$  2^0  +  \/^lT  2g{zo  -  ztd) 

$0  9 


3.1 


where 
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Figure  3-2:  A  sequence  of  photographs  (arranged  in  left-right  order  starting  upper 
left)  taken  during  the  execution  of  a  somersault.  The  robot  is  running  from  left  to 
right.  Approximate  relative  time  of  each  image:  upper  left-0.0  s,  upper  right-0.15  s, 
middle  left-0.33  s,  middle  right-0.66  s,  lower  left-0.80  s,  lower  right-1.02  s. 
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A$  is  the  change  in  somersault  angle  required  during  flight 
$0  is  the  somersault  rate  of  the  body  at  lift-off 
2;o  is  the  height  of  the  center  of  mass  (c.o.m.)  at  lift-off 
Zo  is  the  vertical  velocity  of  the  c.o.m.  at  lift-off 
Zfd  is  estimated  height  of  the  c.o.m.  at  touchdown 
g  is  the  acceleration  of  gravity 

Equation  3.1  relies  on  several  simplifying  assumptions:  1)  the  somersault  dynamics 
are  governed  by  the  planar  equation  /$  =  0,  implying  that  the  external  torques  due  to 
supply  cables  or  wind  resistance  are  negligible,  and  2)  only  the  rotation  in  somersault 
is  signihcant  and  the  somersault  axis  is  a  maximum  principal  axis  of  inertia  so  that 
tilt  and  twist  angles  will  stay  small  if  they  start  small  thus  allowing  us  to  ignore  them, 
and  3)  the  legs  do  not  swing  with  respect  to  the  body  during  the  flight  phase,  so  $o 
represents  the  angular  rates  of  both  the  body  and  the  legs. 

When  the  3D  Biped  robot  somersaults  it  rotates  about  its  major  principal  axis. 
Figure  3-3  shows  that  the  region  of  stable  rotation  about  the  somersault  axis  is  large. 
We  may  then  expect  that  as  long  as  the  somersault  is  initialized  with  the  angular 
momentum  vector  close  to  the  major  principal  axis  then  it  will  remain  close  to  that 
axis.  This  in  turn  means  that  the  tilt  and  twist  angles  of  the  robot  at  landing  will 
be  small  and  the  somersault  dynamics  simplify  as  indicated  above.  In  the  remainder 
of  this  chapter  we  assume  that  this  simplihcation  is  valid  in  computations  involving 
rotational  dynamics  of  the  robot. 

3.3  Somersault  Control  Strategies 

The  goal  of  the  somersault  control  strategy  is  to  produce  a  landing  attitude  that  allows 
the  robot  to  maintain  balance.  To  regain  balance  on  the  landing,  it  is  important  that 
the  robot  achieve  a  desired  horizontal  displacement  of  the  foot  relative  to  the  center 
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Figure  3-3:  The  somersault  axis  of  the  3D  Biped  is  coincident  with  the  major  principal 
axis  of  inertia.  With  the  legs  in  the  fully  extended  position,  as  shown  here,  the  minor 
principal  axis  is  the  ’head-to-toe’  axis. 
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of  mass.  In  the  plane  of  the  somersault  this  horizontal  displacement  is  given  by 

d  =  I  sin  fi  (3-2) 

where  /  is  the  leg  length  and  fi  is  the  leg  inclination  angle,  the  angle  the  leg  makes 
with  vertical  in  the  plane  of  rotation  (Figure  3-4).  If  the  legs  are  held  hxed  relative  to 
the  body,  then  the  desired  landing  attitude,  or  equivalently  A$,  can  be  found  from  d. 
If  the  initial  conditions  of  flight,  $o,  io,  and  Zo,  are  established  accurately  so  that  for 
a  desired  Ztd  and  A$  Equation  3.1  is  satished,  then  the  desired  landing  attitude  will 
be  achieved.  More  generally,  one  approach  for  generating  somersaults  is  to  establish 
accurately  the  initial  conditions  of  flight  to  a  state  that  is  empirically  determined, 
then  to  execute  an  open  loop  pattern  of  actuator  signals  to  produce  components  of  the 
desired  behavior.  Hodgins  and  Raibert  used  this  approach  in  programming  a  planar 
biped  robot  to  perform  front  somersaults  with  a  90%  success  rate.  The  success  of 
such  an  approach  depends  upon  how  precisely  one  can  reproduce  the  state  of  the 
robot  and  how  sensitive  the  desired  movement  is  to  variations  in  the  state. 

For  a  running  robot  with  less  regular  and  repeatable  motion,  such  as  the  3D 
Biped,  precise  initialization  is  more  difficult.  Therefore,  reliable  production  of  a 
desired  landing  attitude  might  be  improved  with  an  in-flight  feedback  strategy  that 
modihes  the  performance  of  the  somersault  based  on  the  state  of  the  robot  in  flight. 

To  the  extent  a  system  with  non-zero  angular  momentum  can  change  its  inertia, 
it  can  also  change  its  rotation  rate.  For  a  somersault,  if  the  time  until  landing  is 
known  then  control  over  rotation  rate  amounts  to  control  over  the  landing  attitude. 
The  basis  of  the  somersault  control  strategy  is  to  change  the  robot  somersault  inertia 
to  produce  a  rotation  rate  that  will  yield  the  desired  body  attitude  at  the  time  of 
landing. 

The  robot  can  change  its  inertia  by  extending  or  retracting  (untucking  or  tucking) 
its  prismatic  legs.  In  flight,  if  the  angular  rate  and  moment  of  inertia  of  the  robot  in 
one  conhguration  are  and  R  then  with  constant  angular  momentum  the  angular 
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Figure  3-4:  In  order  to  maintain  balance  at  landing  the  leg  inclination  angle,  /3,  must 
be  near  the  desired  value.  /3  is  the  angle  the  leg  makes  with  the  local  vertical.  It  is  a 
function  of  the  the  hip  angle,  a  and  the  body  attitude,  $. 

rate  in  a  conhguration  with  inertia  R  is  $2  =  (/i//2)^i-  Joint  limits  restrict  the 
range  of  inertia  and  thus  limit  control  over  the  rotation  rate.  For  the  3D  Biped 
robot,  the  inertia  about  the  somersault  axis  ranges  from  1.22  kg  rrU  with  legs  retracted 
to  lAQkgrrU  with  legs  fully  extended.  In  moving  from  a  fully  retracted  to  a  fully 
extended  position  the  robot  can  reduce  its  somersault  rate  to  81%  of  its  initial  value. 
We  will  refer  to  the  regulation  of  somersault  rate  and  landing  attitude  via  leg  length 
as  the  tuck  servo. 

3.3.1  Control  of  Somersault  Angle 

To  implement  the  tuck  servo  with  state  feedback,  we  hrst  pose  the  requirement  of 
Fquation  3.1  as  a  function  of  the  state  during  flight  rather  than  at  lift-off. 

-  $  %  +  +  2g{zo  -  cos  I3td  ltd) 

$  g 


{t-O  (3.3) 


where 
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^td  is  the  desired  somersault  attitude  at  touch  down, 

$  is  the  current  somersault  attitude, 

$  is  the  current  somersault  rate, 

fitd  is  the  desired  leg  inclination  angle  at  touch  down, 

ltd  is  the  leg  length  at  touch  down, 

t  is  the  current  time, 

to  is  the  lift  off  time. 

The  robot  inertia  and  thus  $  are  functions  of  ltd,  so  we  solve  for  the  ltd  that  will  make 
Equation  3.3  an  equality. 

We  nominally  require  the  hip  angle  to  be  zero  on  landing.  If  this  is  the  case  then 
fCd  =  ^td  and  the  landing  attitude  and  leg  length  determine  the  height  of  the  hip. 
The  center  of  mass  is  coincident  with  the  hip,  so  Ztd  is  given  by 

Ztd  =  cos  ^td  ltd  (3.4) 


We  have  also  assumed  that  the  tilt  angle  at  touchdown  is  zero. 

Dehne  a  function,  /,  as  the  difference  between  the  left  and  right  hand  sides  of 

(3,3). 


J  ~  U  ~  ^o)  ~ 


$ 


(3.5) 


If  the  desired  landing  attitude  is  to  be  achieved  then  /  =  0.  Otherwise,  the  pitch 
rate  needs  to  be  increased  or  decreased  depending  on  the  sign  of  /.  The  control 
algorithm  we  use  changes  the  leg  to  a  length  that  will  make  /  =  0.  The  function,  /, 
is  a  non-linear  function  of  the  leg  length,  so  we  use  a  Newton  search  [Strang  86]  to 
recursively  solve  for  the  desired  leg  length.  The  recursion  uses  the  hrst  order,  Taylor 
series  expansion  of  /  =  0, 

o  =  /(/)  +  |a; 


(3.6) 
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This  yields  the  following  equation  for  A/ 


A/ 


-/ 

df/dl 


(3.7) 


Ideally,  once  the  desired  leg  length  is  achieved,  /  remains  zero  for  the  rest  of  the 
somersault  and  the  desired  leg  length  becomes  the  leg  length  at  touch  down,  /  — )•  ltd- 
During  each  control  cycle  the  tuck  servo  executes  the  following  process 


1.  measure  the  state  of  the  robot 

2.  compute  /  and  df  / dltd 

3.  compute  Altd 

4.  estimate  the  next  value  of  /  using  ltd  =  hd  +  Altd- 

5.  if  /  sa  0  go  to  (6),  else  go  to  (2)  using  new  values  of  /  and  ltd 

6.  servo  the  leg  length  to  Itd- 

To  perform  step  2  above  we  need  to  solve  for  df  /dltd-  Taking  the  derivative  of 
(3.5)  with  respect  to  fd  results  in 


df  /dltd  — 


—  cos  $ 


td 


\jzlT  ‘ig{zo  -  cos  ltd) 


,  —  $ 
$2 


(3,8) 


To  obtain  d^ / dfd  we  use  the  fact  that  the  angular  momentum  is  a  constant  so  that 


where 


$ 


h 

[lo  +  2miC) 


(3.9) 


h  angular  momentum 

r  distance  from  the  lower  leg  c.o.m.  to  the  robot  c.o.m. 


nil  lower  leg  mass 
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lo  robot  pitch  inertia  about  c.o.m.  when  r  =  0. 


The  total  inertia  of  the  robot  as  a  function  of  leg  length  is  (/o  +  2m;r^).  The  2miC 
term  represents  that  inertia  due  to  the  distance  of  the  lower  legs  from  the  robot  c.o.m. 
This  is  the  component  of  inertia  that  we  control  as  the  legs  change  length.  From  (3.9) 
we  get 

d^/dltd  =  -  ,  2\2  (3.10) 

[lo  +  niir^Y 

The  robot  center  of  mass  moves  very  little  when  the  legs  are  extended  or  retracted 
so  we  assume  dr/dlfd  =  1-  Substituting  for  h  from  (3.9)  results  in 


d^jddtd 


2m;r$ 

[lo  +  miC) 


(3.11) 


Now  by  substituting  (3.11),  (3.8),  and  (3.5)  into  (3.7),  we  can  compute  the  required 
Altd- 

The  computation  of  the  desired  leg  length  depends  upon  knowledge  of  robot  pa¬ 
rameters  such  as  inertia  and  leg  mass.  However,  since  the  process  is  repeated  each 
control  cycle  based  upon  the  sensed  state  of  the  robot,  sensitivity  to  precise  knowl¬ 
edge  of  these  parameters  is  reduced.  In  exchange  for  this  robustness  to  uncertainty 
we  give  up  the  ability  to  pre-specify  both  the  landing  attitude  and  the  leg  length  at 
touchdown. 


3.3.2  Accommodating  Landing  Time  Errors 

A  limitation  of  the  tuck  servo  strategy  is  its  dependence  upon  accurate  knowledge 
of  the  time  until  landing.  With  somersault  rates  on  the  order  of  500  deg / small 
errors  in  the  predicted  landing  time  can  result  in  intolerable  landing  attitude  errors. 
Because  we  have  no  measurement  of  the  vertical  position  or  speed  while  airborne  these 
quantities  must  be  derived  from  estimates  of  the  lift-off  conditions.  To  accommodate 
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errors  in  these  estimates,  we  have  modihed  the  landing  strategy  to  allow  the  robot 
to  land  with  the  desired  conhguration  anywhere  within  a  ‘window’  of  the  predicted 
landing  time. 

The  leg  inclination  angle  in  (3.2)  is  the  difference  between  the  somersault  angle 
and  the  hip  angle,  /3  =  $  —  a.  By  changing  the  hip  angle  the  control  algorithm 
makes  the  foot  track  the  desired  displacement  for  a  short  time  just  prior  to  or  after 
the  expected  landing  time.  In  this  way,  the  robot  maintains  the  desired  landing 
conhguration  during  a  ‘landing  window’  that  is  centered  about  the  nominal  landing 
attitude.  In  order  to  maximize  the  landing  window  the  legs  are  moved  ‘back’  relative 
to  the  body  [a  >  0)  early  in  the  hip.  Prior  to  landing,  the  foot  is  swept  forward  in 
order  to  track  the  desired  foot  position  until  touchdown. 

3.3.3  Control  of  Tilt  and  Twist  Angles 

Achieving  a  takeoff  attitude  with  no  tilt  and  twist  is  challenging  with  the  3D  Biped. 
As  the  robot  runs  its  body  somersault  angle  stays  close  to  zero  but  its  body  tilts  and 
twists  in  phase  with  the  stepping  cycle  as  it  runs.  During  a  normal  running  cycle, 
the  robot  uses  hip  actuators  during  stance  to  apply  torques  to  the  body  in  order  to 
control  body  attitude.  However,  we  found  that  it  was  difhcult  to  achieve  small  tilt 
and  twist  angles  and  rates  at  takeoff  using  only  the  hip  servos.  We  found  that  the 
best  method  for  keeping  tilt  and  twist  angles  and  angular  rates  small  during  takeoff 
was  to  use  a  wide  double  stance  during  the  hip  initiation  in  combination  with  the 
normal  attitude  control  used  during  stance  ([Raibert  84]).  This  wide  double  stance 
helps  stabilize  the  robot’s  tipping  motion  during  takeoff.  An  important  component  of 
this  approach  was  to  achieve  simultaneous  touch  down  of  both  feet  at  the  beginning 
of  the  double  stance  phase.  Similarly,  a  wide  double  stance  during  somersault  landing 
and  simultaneous  touchdown  of  both  feet  were  observed  to  minimize  the  effect  of  tilt 
and  twist  errors  on  the  somersault  recovery. 
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3.4  Experiments  with  3D  Biped  Somersaults 

The  3D  Biped  is  a  two-legged  robot  that  is  free  to  translate  and  rotate  in  3D  space. 
Each  leg  is  mounted  adjacent  to  the  center  of  gravity  of  the  body  with  a  ball  and 
socket  hip  joint.  The  hip  joint  allows  leg  rotation  about  the  x,  y,  and  2;  axes 
(±20°,  ±60°,  ±15°  see  Fig.  1).  Hydraulic  actuators  control  each  of  these  degrees 
of  freedom.  The  robot’s  telescoping  legs  contain  a  fourth  hydraulic  actuator  that  acts 
in  series  with  a  pneumatic  compression  spring.  Mass  properties  of  the  3D  Biped  are 
included  in  Table  3.1  Some  of  the  kinetic  energy  of  the  machine  is  stored  in  com¬ 
pression  of  the  air  spring  during  each  bounce,  and  returned  to  power  the  subsequent 
flight  phase.  Energy  is  added  to  the  hopping  oscillation  by  actively  compressing  the 
air  spring  with  the  hydraulic  piston  during  stance.  The  3D  Biped  maintains  balance 
while  running  by  performing  three  control  tasks  [Raibert  84]: 

1.  during  stance,  the  robot  maintains  body  posture  by  applying  hip  torques  be¬ 
tween  the  legs  and  the  body, 

2.  during  stance,  the  robot  adds  energy  to  the  air  spring  to  maintain  the  hopping 
oscillation,  and 

3.  during  flight,  the  robot  positions  the  foot  in  anticipation  of  the  next  stance 
phase  in  order  to  control  forward  velocity. 

To  execute  a  somersault,  the  3D  Biped  modihes  three  steps  in  an  otherwise  normal 
running  sequence.  The  robot  performs  a  hurdle  step  during  which  it  hops  higher  than 
normal  as  it  prepares  to  land  on  both  feet  for  the  flip  step.  The  flip  step  is  initiated 
by  thrusting  with  both  legs  while  pitching  the  body  forward.  During  the  landing  step 
the  robot  lands  on  both  feet  then  resumes  a  normal  running  gait.  The  control  actions 
used  to  execute  the  flip  are  summarized  in  Table  3.2. 

In  laboratory  experiments,  the  3D  Biped  has  successfully  performed  the  forward 
somersault  and  regained  a  stable  running  cycle  afterwards  (Figure  3-2). 
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Table  3.1 

3D  Biped  Parameters 


Symbol 

Description 

Quantity 

m 

total  mass 

31.44  {kg) 

nil 

lower  leg  mass 

0.652  {kg) 

lo 

somersault  inertia 

1.02  {kg  rrU) 

^  min 

min.  lower  leg  radius 

0.404  (m) 

^min 

min.  leg  length 

0.647  (m) 

^max 

max.  leg  length 

0.862  (m) 

Three  sets  of  data  from  a  successful  somersault  and  a  nearly  successful  somersault 
are  shown  in  Figures  3-5,  3-6,  and  3-7.  Figure  3-5  shows  data  for  the  approach,  flip, 
landing,  and  continuation  of  running  for  a  successful  somersault.  The  nominal  desired 
landing  attitude  was  set  to  350°.  The  robot  ran  steadily  until  the  hurdle  step  at  which 
time  it  hopped  higher  than  normal  as  it  prepared  to  land  on  both  feet.  During  the 
flip  step  the  body  is  thrust  upwards  and  accelerated  in  somersault.  The  desired  leg 
inclination  angle  was  set  to  —5.7°  based  on  the  forward  speed  in  flight.  The  robot 
lands  0.080s  earlier  than  anticipated  with  a  somersault  attitude  of  325°  and  with  a 
leg  inclination  angle  of  —11°.  The  actual  leg  inclination  angle  is  much  closer  to  its 
desired  value  than  the  corresponding  values  in  somersault  because  of  the  feet-back 
position  in  flight.  The  foot  positioning  servo  was  not  used  in  this  somersault  because 
the  desired  foot  position  was  always  slightly  in  back  of  the  actual  foot  position,  and 
the  legs  were  already  in  a  swept  back  conhguration.  Balance  is  regained  on  landing, 
but  since  the  robot  lands  with  a  slight  backward  lean,  forward  speed  is  lost.  Forward 
speed  and  posture  are  quickly  restored  during  the  following  steps. 

Note  the  oscillation  in  tilt  angle  during  flight.  This  oscillation  occurs  because  the 
robot  took  off  with  a  non-zero  tilt  rate.  Since  rotation  about  the  somersault  axis 
is  passively  stable,  the  oscillation  does  not  grow.  We  found  that  in  order  to  regain 
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Table  3.2 

Control  Summary  for  the  Somersault 


Step 

Action 

Approach 

Run  forward  (sa  Im/s)  with  alternating  gait 

Hurdle 

Pitch  up  slightly 

Hop  with  increased  leg  thrust 

Prepare  to  land  simultaneously  on  both  feet 

Flip 

Jump  with  maximum  thrust 

Pitch  body  forward  with  maximum  torque 
Shorten  legs  once  airborne 

Servo  hips  to  feet  back  position 

Fngage  tuck  servo 

Prepare  to  land  simultaneously  on  both  feet 
Track  desired  foot  position 

Landing 

Dissipate  energy  on  landing 

Return  somersault  rate  to  zero,  restore  posture 
Adjust  nominal  leg  length  based  on  Sd 

Following 

Resume  running  with  alternating  gait 
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balance  after  the  somersault,  the  tilt  angle  on  landing  must  be  kept  moderately  small, 
|0|  <  15°.  To  do  this  we  increased  the  stance  width  during  the  flip  step,  thereby 
providing  a  passively  stable  stance  conhguration  in  tilt  during  initiation. 

The  data  in  Fig.  4  shows  the  somersault  action  on  a  larger  scale  to  illustrate  the 
function  of  the  tuck  servo.  Between  the  beginning  of  the  flip  step  and  lift  off,  the 
magnitude  of  the  somersault  rate,  hip  angle  rate,  and  leg  length  all  increase.  At  lift 
off,  the  somersault  rate  of  the  body  declines  rapidly  as  the  legs  are  accelerated  to 
the  rotation  rate  of  the  body.  The  conservation  of  momentum  constraint  produces 
the  symmetry  between  the  absolute  angular  rate  of  the  body  and  the  relative  angular 
rate  of  the  hips.  During  this  time,  the  legs  are  tucked  to  the  shortest  possible  length. 

The  tuck  servo  is  engaged  as  the  hips  reach  the  desired  feet-back  position  and 
come  to  rest  relative  to  the  body.  The  robot  has  a  somersault  rate  of  606  deg  js  at 
the  time  the  tuck  servo  is  engaged.  At  this  somersault  rate,  it  is  estimated  that  the 
robot  will  over-rotate  by  84°.  This  error  is  illustrated  by  the  third  graph  of  Fig.  4 
which  compares  the  estimated  time  until  touchdown  and  the  estimated  time  until 
the  desired  somersault  attitude  is  achieved.  The  tuck  servo  extends  the  legs  to  the 
maximum  possible  length  to  slow  down  the  somersault  rotation  to  4:63  deg  js  at  which 
point  it  is  estimated  that  the  robot  will  land  at  nearly  the  desired  attitude.  The  robot 
maintains  this  conhguration  until  landing. 

Data  from  another  nearly  successful  somersault  is  shown  in  Fig.  5.  In  this  somer¬ 
sault,  the  robot  uses  the  foot  positioning  servo  to  keep  from  over-rotating.  Once  again 
the  desired  somersault  landing  attitude  was  350°.  The  desired  leg  inclination  angle 
was  set  to  —5.9°.  The  robot  lands  0.065s  later  than  anticipated  with  a  somersault 
attitude  of  391°  and  with  a  leg  inclination  angle  of  —3.2°.  As  the  robot  detected 
that  it  was  over-rotating  it  swept  the  feet  forward  quickly  to  track  the  desired  foot 
position.  This  increased  the  somersault  rate  because  of  the  conservation  of  angular 
momentum  and  contributed  to  the  somersault  attitude  error  on  landing.  A  velocity 
measurement  error  during  stance  after  landing  led  to  the  loss  of  forward  speed  of  the 
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Table  3.3 

3D  Biped  Flip  Attempts 


File 

Outcome 

l^td 

deg 

deg 

$ 

deg/s 

io  (est) 
m/s 

io  (meas) 
m/s 

io%err 

speed;„ 

m/s 

92.346.4 

success 

0.12 

-10.8 

547 

3.08 

3.04 

-1.2 

1.04 

92.346.5 

success 

-0.14 

-6.29 

529 

3.48 

3.33 

-4.3 

0.77 

92.346.6 

success 

0.33 

-8.47 

524 

3.41 

3.32 

-2.8 

0.89 

92.346.7 

success 

-0.31 

-13.7 

549 

3.53 

3.48 

-1.4 

1.10 

92.346.8 

6  steps 

-1.95 

-18.9 

553 

3.69 

3.70 

1.6 

1.02 

92.346.9 

success 

-0.39 

-12.7 

551 

3.45 

3.47 

0.7 

1.04 

92.346.10 

success 

-0.04 

-11.6 

545 

3.54 

3.48 

-  1.7 

0.96 

92.346.11 

success 

-0.28 

-15.3 

541 

3.15 

3.13 

-  0.5 

0.90 

92.346.12 

fall 

-3.42 

-1.04 

529 

3.28 

3.04 

-  7.4 

0.89 

92.346.13 

fall 

-1.87 

-17.7 

503 

3.33 

3.20 

-3.9 

0.88 

robot  so  that  it  was  momentarily  supported  by  safety  ropes  before  resuming  running 
and  therefore  not  considered  a  complete  success. 

Table  3.3  shows  data  compiled  for  ten  somersault  attempts  performed  by  the  3D 
Biped  in  the  laboratory.  The  robot  successfully  regained  balanced  running  on  seven 
of  these  attempts. 


3.5  Summary 


In  this  chapter,  I  discuss  a  strategy  for  robot  somersaults  that  combines  elements  of 
feed  forward  control,  feedback  control,  and  passive  dynamic  stability.  I  also  presented 
results  from  somersault  experiments  done  in  the  laboratory  on  a  3D  biped  running 
robot.  The  somersault  is  initialized  using  pre-programmed  patterns  of  action.  In 
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Figure  3-6:  Somersault  angle,  absolute  somersault  rate  (solid)  and  relative  hip  rate 
(dashed).  The  third  graph  shows  the  estimated  time  until  touch  down  (solid)  and 
estimated  time  until  desired  somersault  attitude  is  achieved  (dashed).  The  last  graph 
is  the  measured  leg  length  (solid)  and  desired  leg  length(dashed)  for  one  leg.  Vertical 
lines  indicate  initiation  of  flip,  lift-off,  tuck  servo  initiation,  and  landing.  Data  hie 
B92.181.3. 
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Figure  3-7:  This  somersault  was  over-rotating  so  the  feet  were  swept  forward  at  the 
end  to  track  a  desired  foot  position.  Data  shows  absolute  somersault  rate  (solid)  and 
relative  hip  rate  (dashed)  and  estimated  time  until  touch  down  (solid)  and  estimated 
time  until  desired  somersault  attitude  is  achieved  (dashed).  Vertical  lines  indicate 
initiation  of  flip,  lift-off,  tuck  servo  initiation,  and  landing.  Data  hie  B92.181.4. 
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flight,  a  feedback  strategy  changes  robot  inertia  to  control  the  landing  attitude  of 
the  somersault.  The  robot  actively  positions  its  feet  to  maintain  a  desired  landing 
conhguration  during  an  interval  surrounding  the  predicted  landing  time.  The  passive 
tilt  stability  inherent  in  a  wide  double  stance  is  used  to  reduce  tilt  angle  and  rate  on 
takeoff.  The  passive  stability  of  a  rigid  body  rotating  about  its  maximum  principle 
axis  of  inertia  accounts  for  moderate  tilt  angles  on  touch  down  given  moderate  tilt 
angles  and  rates  on  lift-off. 


Chapter  4 


Passively  Stable  Layout 
Somersaults 

4.1  Introduction 

The  layout  somersault  is  an  airborne  maneuver  in  which  the  performer  rotates  about 
a  side-to-side  axis  while  maintaining  an  erect  body  conhguration.  The  layout  somer¬ 
sault  is  often  considered  to  be  inherently  unstable  because  it  involves  rotation  about 
the  middle  principal  axis  of  inertia,  an  unstable  mode  of  rotation  for  a  rigid  body. 
I  found  that  passive  arm  movement  could  neutrally  stabilize  the  layout  somersault. 
These  passive  movements  are  generated  by  dynamic  forces  that  arise  from  body  mo¬ 
tion  and  from  compliant  shoulders.  Figure  4-1  shows  a  sketch  of  the  type  of  three-body 
model  used  for  simulation,  analysis,  and  experimental  tests  in  the  laboratory.  This 
simple  three-body  model  of  a  human  can  exhibit  passively  stable  layout  somersaults 
if  an  appropriately  tuned  spring  is  used  at  the  shoulder. 

Layout  somersaults  are  considered  stable  if  steady  somersaulting  is  achieved  with¬ 
out  exhibiting  any  ’large’  oscillations  in  the  tilt  or  twist  angles.  For  example.  Fig¬ 
ure  4-2  shows  simulation  data  of  the  simple  human  model  for  two  cases: 

1.  Rigid  body  (dashed  lines):  this  model  has  very  stiff  shoulder  springs  and 
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Figure  4-1:  Diagram  of  a  dynamic  model  used  to  study  layout  somersaults.  The  head, 
torso,  and  legs  comprise  a  single  rigid  body.  The  arms  are  connected  to  the  body 
with  pin  joints  that  allow  rotation  about  the  1  body  axis.  Joint  torques  are  provided 
by  torsional  springs  and  (possibly)  dampers. 
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dampers  to  approximate  the  rigid  body  case.  As  with  the  true  rigid  body,  the 
tilt  and  twist  angles  periodically  move  through  large  excursions  from  zero  as 
the  body  somersaults. 

2.  Passively  stable  body  (solid  lines):  this  model  has  tuned  shoulder  springs 
without  dampers  that  allow  considerable  passive  arm  movement.  The  tilt  and 
twist  angles  of  the  body  stay  near  zero. 

The  rigid  body  tilt  and  twist  angles  shown  in  Figure  4-2  are  shown  in  another 
form  in  the  map  of  Figure  2-7  which  is  repeated  here  in  Figure  4-3.  The  rigid  body 
instability  is  evident  by  the  trajectories  that  converge  and  diverge  upon  the  somersault 
axis.  The  twisting  oscillation  of  the  simulation  data  corresponds  to  a  trajectory  on 
the  map  that  is  centered  about  the  maximum  principal  axes  (front-to-back). 

In  this  chapter,  I  present  results  of  studies  on  the  passive  layout  somersault  using 
non-linear  dynamic  simulation  and  linear  stability  analysis.  By  deriving  a  linearized 
model  of  the  dynamics  of  a  three-body  model,  I  show  that  passive  layout  stability 
depends  upon  both  the  arm  angle  and  the  shoulder  spring  constant.  I  show  condi¬ 
tions  under  which  passive  layout  stability  can  theoretically  be  achieved.  Using  the 
simplest  possible  model  of  the  layout  somersault.  I  explain  the  fundamental  dynam¬ 
ics  of  passive  stabilization.  I  also  present  summarized  results  of  non-linear  dynamic 
simulations  that  support  the  linear  stability  analysis. 

4.2  Simple  Human  Model 

Figure  4-1  shows  a  simplihed  human  body  model  used  to  study  the  dynamic  stability 
of  somersaulting  motion.  From  the  non-linear  equations  of  motion  of  this  model  I 
analytically  derive  linearized  equations  that  govern  the  motion  of  the  body  relative  to 
pure  somersaulting  rotation.  The  linear  equations  of  motion  decouple  into  two  distinct 
subsystems  which  allow  further  model  simplihcations.  When  inertial  parameters  for 
a  human  are  used  in  this  model,  a  linear  stability  analysis  indicates  that  stability  of 
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Figure  4-2:  Number  of  complete  somersault  revolutions,  twist  and  tilt  Euler  angles 
describing  body  attitude,  and  left  and  right  shoulder  angles  for  two  simulations.  The 
solid  lines  indicate  the  passively  stable  case.  The  dashed  lines  indicate  the  rigid  body 
case.  The  data  for  the  passively  stable  case  shows  nearly  eight  stable  somersaults 
about  the  middle  principal  axis. 


CHAPTER  4.  PASSIVELY  STABLE  LAYO  UT  SOMERS  A  ULTS 


74 


Figure  4-3:  A  map  of  rigid  body  rotation  for  a  human  in  the  layout  position.  The 
inertially  hxed  angular  momentum  vector  (black)  paints  trajectories  onto  the  surface 
of  the  sphere  as  the  sphere  rotates.  The  axes  of  the  sphere  and  the  principal  axes  of  the 
body  remain  parallel  as  the  body  rotates.  Each  trajectory  on  the  sphere  corresponds 
to  a  different  rotational  energy.  The  trajectories  indicate  the  tilt  and  twist  angles 
of  the  body  as  it  rotates.  This  map  does  not  include  somersault  angle  and  does  not 
show  the  time  dependence  of  the  tilt  and  twist  Euler  angles. 
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the  layout  somersault  can  be  achieved  for  a  range  of  nominal  arm  angles  and  shoulder 
spring  constants. 

The  human  body  model  has  shoulder  pin  joints  that  allow  the  arms  to  be  raised 
and  lowered  in  the  ‘2-3’  plane  of  the  body.  The  head,  torso,  and  legs  are  modeled  as 
a  single  rigid  body.  Torques  at  the  shoulder  joints  are  provided  by  torsional  springs 
and  dampers.  Neglecting  translation  of  the  center  of  mass,  the  model  has  hve  degrees 
of  freedom  (d.o.f.):  three  for  rotation  of  the  body  with  respect  to  inertial  space  and 
two  for  relative  movement  of  each  arm.  The  orientation  of  the  body  with  respect 
to  an  inertial  coordinate  frame  is  described  using  a  2-1-3  sequence  of  Euler  angles 
(Figure  2-1.)  To  describe  the  attitude  of  the  body  relative  to  the  inertial  frame, 
a  coordinate  system  initially  parallel  to  the  inertial  reference  frame  is  hrst  rotated 
through  the  somersault  angle  about  the  inertially  hxed  ‘2’  axis,  then  rotated  through 
the  tilt  angle  about  the  intermediate  ‘1’  body  axis,  and  hnally  rotated  through  the 
twist  angle  about  the  body  hxed  ‘3’  axis. 

The  rotation  of  the  right  and  left  arms  relative  to  the  body  are  given  respectively 
by  the  angles  and  7;.  An  arm  angle  is  equal  to  zero  when  the  hand  is  held  next 
to  the  thigh.  Positive  rotation  of  the  arm  about  the  ‘1’  axis  is  given  by  the  right 
hand  rule.  The  right  and  left  arm  angles  are  rewritten  in  terms  of  symmetric  (7^) 
and  asymmetric  (7^)  components  as  follows: 

7s  =  1/2(7/  -  7r) 

7a  =  1/2(7/  +7r) 

The  parameters  necessary  to  describe  the  three-body  system  in  the  non-linear 
equations  of  motion  are  given  below.  The  body  is  assumed  to  be  symmetric  from  left 
to  right  when  the  arms  are  held  in  a  symmetric  orientation. 
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Figure  4-4:  Diagram  of  a  five  d.o.f.  model  with  body  axes  and  vectors  labeled. 


nib 


mi 

/6i,  Ib2,  lbs 
III,  Ih,  II3 
rbi,  rbs,  rbs 
rli,  r/2,  rls 
ksh  ■)  bsh 


mass  of  the  body 
mass  of  arms 

body  principal  inertias  in  body  hxed  axes 
arm  principal  inertias  in  arm  hxed  axesF^o  =  0.0 
vector  components  from  body  c.g.  to  left  arm  pin  joint 
vector  components  from  left  arm  c.g.  to  pin  joint 
shoulder  spring  and  damping  constants 


Figure  4-4  shows  vectors  and  body  axes  dehned  in  the  above  list. 

In  order  to  study  the  layout  somersault  as  performed  by  humans,  I  will  use  data 
for  a  human  performer,  Carl  Furrer,  1982  World  Trampoline  Champion  [Yeadon  84]. 
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Table  4.1 

Inertia  Data  for  Carl  Furrer 


Parameter 

Value 

Units 

nib 

58.3 

kg 

nil 

3.743 

kg 

IW 

10.19 

kg  rrC 

9.93 

kg  rrC 

lbs 

0.52 

kg  rrC 

Ik 

0.128 

kg  rrC 

Ik 

0.128 

kg  rrC 

Ik 

0.0037 

kg  rrC 

rbi 

0.0 

m 

rbs 

0.18 

m 

rbs 

0.479 

m 

rk 

0.0 

m 

rk 

0.0 

m 

rk 

0.262 

m 

The  data  for  Furrer ’s  body  parameters  are  included  in  Table  4.1. 

4.2.1  Nonlinear  Equations  of  Motion 

I  derive  the  non-linear  equations  of  motion^  as  a  starting  point  for  a  linear  analysis 
that  will  follow.  I  used  Kane’s  method  [Kane  85]  to  derive  the  equations.  Rather  than 
include  the  lengthly  equations  here  I  include  the  Mathematica  code  used  to  produce 
them  in  Appendix  A.l.  The  equations  are  derived  using  a  reference  coordinate  frame 
located  at  the  center  of  mass  of  the  three  body  system.  Therefore,  the  rotational 
equations  are  decoupled  from  the  translational  equations. 

^Non-linear  equations  used  for  dynamic  simulation  were  derived  using  a  commercially  available 
package,  SD-FAST.  This  version  of  the  equations  was  not  amenable  to  analytic  linearization. 
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4.2.2  Linearized  Equations  of  Motion 

The  linearized  equations  will  describe  the  motion  of  the  system  with  respect  to  a 
reference  frame  that  is  steadily  rotating  about  the  2  axis  with  rate  Pure  somersault 
rotation  about  a  principal  axis  of  inertia  is  an  equilibrium  solution  for  this  system  if 
a  constant  feed  forward  shoulder  torque  is  applied  to  cancel  the  centrifugal  forces  due 
to  steady  rotation.  This  steady  torque  can  be  achieved  by  pre-tensioning  the  shoulder 
spring.  The  magnitude  of  this  torque  will  be  a  function  of  the  nominal  symmetric 
arm  angle,  Tgo,  and  the  rotation  rate,  $.  Any  deviation  of  the  arm  from  T^o  will 
result  in  additional  shoulder  torques  from  the  springs  and  dampers. 

The  linearized  equations  of  motion  for  this  model  follow: 


where 


M£  +  ^(G  +  B)i  +  +  K')  X  =  0 


M 


G 


mil  0  mi3  0  0 

0  m22  0  0  0 

mm  0  m33  0  0 

0  0  0  mil  0 

0  0  0  0  mss 

0  gi2  0  0  0 

-gi2  0  5-23  0  0 

0  ~5'23  0  0  0 

0  0  0  0  5'4s 

0  0  0  — t/rs  0 


(4.1) 


(4.2) 


(4.3) 
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0  0  0  0  0 

0  0  0  0  0 

D  =  0  0  2Uh/^  0  0 

0  0  0  0  0 

0  0  0  0  2Uh/^ 

kii  0  ki3  0  0 

0  7^22  0  0  0 

K  =  0  A;33  0  0 

0  0  0  0  0 

0  0  0  0  7^55 

0  0  0  0  0 

0  0  0  0  0 

K' =  0  0  2kshl^^  0  0  (4.6) 

0  0  0  0  0 

0  0  0  0  2kshl^^ 

The  state  vector,  x,  is  comprised  of  (in  order):  9  -  the  tilt  angle,  tj)  -  the  twist 
angle,  7^  -  the  asymmetric  deviation  of  the  two  arms  from  the  nominal  arm  angle, 

(f)  -  the  deviation  of  the  somersault  angle  from  the  frame  steadily  rotating  at  rate  $, 
and  7s  -  the  symmetric  deviation  of  the  two  arms  from  the  nominal  arm  angle. 

The  expressions  for  individual  components  of  Af,  G,  and  K  are  listed  in  Appendix‘A.3. 
In  Equation  4.1  the  nominal  rotation  rate,  $,  scales  the  rate  terms  and  scales  the 
spring  terms.  This  is  evidence  of  the  fact  that  the  coriolis  forces  involving  somersault 
rate  give  rise  to  rate  dependent  terms  and  the  centrifugal  forces  involving  somersault 
rate  give  rise  to  spring  terms.  The  terms  in  D  and  K'  arise  from  the  spring  and 
damper  model  of  the  shoulder  torques  for  each  arm. 


(4.4) 


(4..5) 


Ti  =  -kskir/i  -  Fjo)  -  b.kii 


(4.7) 
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'^r  —  -j-  Tso)  bsh'^r 


(4,8) 


Non-dimensional  states  (z^z/z)  are  derived  by  scaling  i,  and  i  by  $  or  as 
follows 


Z  =  X 

z  =  i/d) 
i  =  i/(i)2 


Rewriting  the  equations  of  motion  as  functions  of  the  non-dimensional  states  results 
in  the  following 

{MS  +  {G  +  D)i  +  {K  +  K')  z)  =  0  (4.9) 

It  is  now  clear  that  $  serves  only  as  a  time  scale  of  Equation  4.9.  We  may  consider 
$  =  1  for  simplicity  in  studying  stability  as  long  as  we  are  willing  to  work  with 
scaled  shoulder  spring  and  damping  terms.  In  the  following  stability  analysis  we  use 
a  non-dimensional  shoulder  spring  dehned  as  follows 


where  Rh  is  the  inertia  of  the  arm  about  the  shoulder.  For  Furrer’s  data  Rh  = 
0.3849^(5'  irA. 

The  reason  for  describing  arm  angles  in  terms  of  symmetric  and  asymmetric  com¬ 
ponents  is  clear  upon  examination  of  the  matrices  of  Equations  4.2  -  4.6.  This  de¬ 
composition  allows  us  to  see  that  in  general  the  linearized  system  decouples  into  two 
distinct  subsystems;  the  upper  3-by-3  system  that  couples  tilt  and  twist  with  asym¬ 
metric  arm  movement  and  the  lower  2-by-2  system  that  couples  rotation  rate  with 
symmetric  arm  movement.  We  consider  these  two  subsystems  separately  in  a  linear 
stability  analysis.  It  will  turn  out  that  the  parameters  which  stabilize  the  composite 
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system  are  nearly,  but  not  exactly,  identical  to  those  that  stabilize  the  3-by-3  system. 

4.3  Stability  Analysis 

In  this  section,  I  present  a  linear  stability  analysis  of  the  layout  somersault.  To  start, 
I  present  the  stability  of  the  equivalent  rigid  body,  that  is  the  rigid  body  that  would 
result  if  the  arms  were  rigidly  hxed  to  the  body.  As  expected,  if  the  body  were  rigid 
the  layout  somersault  would  be  unstable.  Then  I  proceed  to  an  analysis  of  the  three- 
body  system.  The  goal  is  to  determine  the  shoulder  spring  constants  and  arm  angles 
that  produce  stable  somersaults  when  the  remaining  body  parameters  are  hxed.  A 
root  locus  plot  shows  how  the  location  of  the  linear  system  poles  move  as  the  shoulder 
spring  constant  is  systematically  varied.  The  root  locus  plot  shows  that  under  some 
conditions,  all  poles  of  the  linear  system  are  simultaneously  on  the  imaginary  axis, 
implying  (neutral)  stability.  The  results  of  the  root  locus  are  also  presented  in  the 
form  of  stability  diagrams  that  show  under  what  values  of  shoulder  spring  and  arm 
angle  the  passive  stability  can  be  realized. 

4.3.1  Rigid  Body 

As  a  start,  let  us  examine  the  characteristics  of  the  equivalent  rigid  body  for  this 
system.  I  consider  the  equivalent  rigid  body  to  be  the  rigid  body  that  has  the  same 
principal  inertia  as  the  whole  three  body  system.  If  we  consider  the  arms  to  be  hxed 
rigidly  to  the  body  then  we  can  solve  for  the  equivalent  rigid  body  inertias,  and  thus 
the  inertia  ratios,  ki  and  ks,  (see  Section  2.3.2)  as  a  function  of  the  nominal  arm 
angle,  Tso-  The  inertia  ratios  ki  and  ks  capture  the  inertia  of  the  equivalent  rigid 
body  and  therefore  determine  its  stability. 

For  a  true  rigid  body  rotating  about  its  intermediate  axis,  the  unstable  mode  of 
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motion  will  grow  at  a  rate  given  by  Equation  2.13 

2;  =  where  s  =  (  — 

Since  is  simply  the  nominal  somersault  angle,  the  change  in  somersault  angle,  $, 
required  for  the  unstable  mode  to  grow  by  a  factor  of  N  is  given  by  Equation  2.14 


Using  the  data  for  Eurrer  (Table  4.1)  in  the  three  d.o.f.  model  with  a  symmetric 
arm  angle  of  r^o  =  1.3  rad.  the  resulting  rigid  body  inertia  ratios  are  ki  =  0.65  and 
ks  =  —0.87.  Using  these  numbers  with  =  10  results  in  $  =  3.1  rad,  or  the  unstable 
mode  will  grow  by  a  factor  of  10  in  less  than  one  half  somersault.  Eigure  4-5  shows 
a  plot  of  the  inertia  ratios  of  the  hve  d.o.f.  model  for  symmetric  arm  angles  between 
0.0  and  3.0  radians.  These  inertia  ratios  are  superposed  upon  curves  of  constant  $. 

4.3.2  Multi-Body 

Some  useful  information  regarding  the  stability  of  this  system  is  available  by  a  simple 
examination  of  the  form  of  the  matrices  in  Equation  4.1.  If  shoulder  damping,  bsh^  is 
zero  then  D  =  0  and  Equation  4.1  describes  a  conservative  gyric  system.  As  discussed 
in  Section  2.3.1,  the  strongest  stability  result  one  may  expect  from  a  conservative  gyric 
system  is  neutral  stability,  i.e.  the  system  will  be  a  perfect  oscillator.  A  sufficient 
condition  for  stability  is  for  the  net  stiffness  matrix  to  be  positive  dehnite,  K-\-K'  >  0. 
This  is  not  possible  since  K  +  K'  is  not  even  full  rank  for  the  hve  d.o.f.  system.  This 
happens  because  the  system  is  unstable  to  perturbations  in  somersault  rate  as  was 
the  case  for  the  rigid  body  system  in  Chapter  2.  However,  even  the  three-by-three 
subsystem  of  tilt,  twist,  and  asymmetric  arm  movement  does  not  have  a  positive 
dehnite  stiffness  matrix  for  the  hxed  body  parameters  used  here.  Therefore,  if  the 
system  is  to  be  stable  it  must  be  gyrically  stabilized. 
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Figure  4-5:  Rigid  body  stability  diagram  for  Furrer  s  data.  The  dark  lines  show  ki 
and  ks  for  0  <  F^o  <  3.0.  The  light  lines  are  lines  of  constant  the  change  in 
somersault  angle  required  for  the  unstable  rigid  body  mode  to  grow  by  a  factor  of 
ten.  The  plot  shows  that  the  human  form  is  relatively  unstable  in  the  rigid  body 
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Figure  4-6:  Root  locus  plot  for  the  3-by-3  subsystem  of  tilt,  twist,  and  asymmetric 
arm  movement.  =  0.0  and  the  nominal  symmetric  arm  angle,  r^o,  is  1.3  rad. 

Light  dots  indicate  root  locations  for  Ah/i  ~  —2.0.  Dots  get  progressively  darker  as 
Ksh  increases  to  4.0. 

Root  Locus 

In  this  section  I  present  plots  of  the  roots  of  the  characteristic  equation  of  4.9  as  the 
shoulder  spring  constant  and  nominal  symmetric  arm  angle  are  systematically  varied. 
Figure  4-6  is  a  root  locus  for  the  3-by-3  subsystem  of  tilt,  twist,  and  asymmetric 
arm  angle.  This  plot  shows  how  the  3-by-3  system  roots  move  as  the  shoulder  spring 
stiffness  is  gradually  increased  when  F^o  =  1.3  rad.  Note  the  symmetry  about  the 
imaginary  axis  as  predicted  for  a  conservative-gyric  system. 

When  the  nondimensional  shoulder  spring,  Ah/i,  is  set  to  -2.0  the  four  symmet¬ 
rically  located  roots  have  both  real  and  imaginary  components.  As  the  shoulder 
stiffness  is  increased  the  symmetric  right  and  left  half  plane  roots  converge  on  the 
imaginary  axis  where  they  split.  One  set  moves  away  from  the  origin  along  the  imag¬ 
inary  axis  for  as  long  as  shoulder  stiffness  increases.  The  other  pair  converge  at  the 
origin  where  they  again  split  one  traveling  along  the  positive  real  axis  the  other  along 
the  negative  real  axis.  At  very  high  shoulder  stiffness  these  two  real  roots  approach 
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the  real  (stable  and  unstable)  rigid  body  roots.  At  high  shoulder  stiffness  the  high 
frequency  imaginary  pair  correspond  to  an  arm-body  oscillation. 

There  are  two  poles  of  this  system  at  that  do  not  move  with  changing  shoulder 
spring  constant.  These  are  the  roots  of  the  stroboscopic  mode  of  rotation.  This  mode 
is  a  rigid  body  mode  of  motion  that  does  not  include  arm  movement.  Examination 
of  the  eigenvectors  associated  with  this  mode  reveals  that  the  attitude  and  velocity 
vectors  are  parallel.  This  indicates  that  this  mode  involves  pure  spin  about  the  middle 
principal  axis,  the  equilibrium  condition,  when  that  axis  is  perturbed  slightly  from  the 
original  orientation.  This  rigid  body  mode  will  exist  for  all  values  of  the  parameters. 

For  some  values  of  shoulder  spring  all  of  the  roots  of  this  system  are  located  on 
the  imaginary  axis.  These  marginally  stable  conditions  suggest  that  for  these  values 
of  the  parameters  the  physical  system  may  be  stable.  Note  that  even  though  equal 
size  steps  in  K^h  were  used  to  numerically  evaluate  the  root  locus,  the  loci  do  not 
move  in  even  steps.  In  particular  the  loci  ‘jump’  to  the  imaginary  axis  from  the  right 
and  left  half  planes.  This  indicates  that  the  edge  of  the  stable  region  in  parameter 
space  is  steep. 

Figure  4-7  shows  a  root  locus  for  the  same  system  with  the  addition  of  a  viscous 
damper  at  the  shoulder  joint.  Note  the  loss  of  symmetry  about  the  imaginary  axis 
and  the  lack  of  a  set  of  spring  values  that  stabilize  the  system.  This  result  agrees 
with  the  spirit  of  Hughes  proof  that  damping  tends  to  destabilize  gyrically  stabilized 
systems.  However,  this  model  does  not  exactly  ht  the  requirements  of  his  proof  since 
the  damping  matrix  is  not  positive  dehnite.  Finally,  Figure  4-8  shows  how  the  root 
loci  plots  change  as  T^o  is  varied  between  0.0  and  2.8  radians. 

The  root  loci  of  this  section  suggest  that  neutral  stability  of  the  layout  somer¬ 
sault  is  possible  for  certain  combinations  of  nominal  arm  angle  and  shoulder  spring 
constant.  It  appears  that  asymptotic  stability  is  unlikely  as  shoulder  damping  tends 
to  destabilize  the  system.  In  the  next  section  we  will  examine  more  precisely  the 
dependence  of  stability  on  arm  angle  and  shoulder  spring  constant. 
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Figure  4-7:  Root  locus  plot  for  the  3-by-3  subsystem  of  tilt,  twist,  and  asymmetric 
arm  movement.  =  0.1  and  the  nominal  symmetric  arm  angle,  r^o,  is  1.3  rad. 

Light  dots  indicate  root  locations  for  Ah/i  ~  —2.0.  Dots  get  progressively  darker  as 
Ksh  increases  to  4.0. 

Stability  Diagrams 

In  this  section,  I  present  plots  that  show  precisely  how  stability  depends  on  the 
shoulder  spring  constant  and  nominal  arm  angle. 

Figure  4-9  shows  the  results  of  a  search  in  the  Fso-Rsh  parameter  space  for  condi¬ 
tions  that  stabilize  the  hve  d.o.f.  linearized  system.  This  collection  of  points  suggests 
a  region  in  the  parameter  space  that  stabilizes  the  system.  The  stable  points  in  Fig¬ 
ure  4-9  are  an  intersection  of  the  stable  points  of  the  two  decoupled  linear  systems 
that  comprise  the  hve  d.o.f.  system.  Figures  4-10  and  4-11,  are  the  corresponding 
plots  of  stable  points  for  the  3-by-3  subsystem  of  tilt-twist-asymmetric  arm  and  the 
2-by-2  subsystem  of  somersault-symmetric  arm.  Comparison  of  Figures  4-9  and  4-10 
shows  that  stability  of  the  hve  d.o.f.  system  is  well  represented  by  the  stability  of 
the  simpler  3-by-3  system,  i.e.  few  conditions  that  stabilize  the  3-by-3  system  fail  to 
stabilize  the  5-by-5  system. 
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Figure  4-8:  Root  loci  for  the  3-by-3  subsystem  of  tilt,  twist,  and  asymmetric  arm 
movement  as  the  nominal  symmetric  arm  angle,  r^o,  is  varied  between  0.0  (top  left) 
and  2.8  rad  (bottom  center)  ,  bsh  =  0,  —2.0  <  Kgh  <  4.0. 
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Figure  4-9:  Each  darkened  point  in  this  plot  corresponds  to  a  specihc  choice  of  Ksh 
and  the  nominal  symmetric  arm  angle  (Fso)  that  stabilize  the  full  hve  d.o.f.  linear 
system,  =  0. 


Figure  4-10:  Each  darkened  point  in  this  plot  corresponds  to  a  specihc  choice  of 
Ksh  and  the  nominal  symmetric  arm  angle  (Fso)  that  stabilize  the  three  d.o.f.  linear 
system  of  tilt,  twist,  and  asymmetric  arm  movement,  bshl^  =  0. 
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Figure  4-11:  Each  darkened  point  in  this  plot  corresponds  to  a  specihc  choice  of  Ksh 
and  the  nominal  symmetric  arm  angle  (Fso)  that  stabilize  the  two  d.o.f.  linear  system 
of  somersault  and  symmetric  arm  movement,  =  0.  Higher  values  of  Ah/i  than 

are  shown  on  this  plot  are  suspected  to  also  stabilize  the  system. 
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4.3.3  Non-Linear  Stability  via  Dynamic  Simulation 

The  linear  stability  analysis  (Figure  4-9)  predicts  that  passive  layout  somersault  sta¬ 
bility  can  be  achieved  for  nominal  symmetric  arm  angles  between  0.8  and  2.0  radians 
(assuming  only  positive  shoulder  spring  constants  are  allowed.)  We  would  normally 
expect  the  linear  analysis  results  to  be  valid  only  in  some  region  of  the  linearizing 
condition.  Stability  to  perturbations  away  from  the  linearizing  condition  will  depend 
upon  the  size  of  the  perturbations.  Since  in  the  case  of  passive  layout  somersaults 
the  linear  analysis  can  only  predict  neutral  stability  we  can  not  conclude  that  the 
non-linear  system  is  stable  even  for  arbitrarily  small  perturbations.  While  non-linear 
dynamic  simulation  can  not  prove  non-linear  stability  it  can  provide  a  fast  and  easy 
check  on  whether  or  not  the  non-linear  system  is  likely  to  be  stable.  In  this  section, 
I  use  non-linear  dynamic  simulation  as  a  check  on  the  linear  stability  results. 

Figure  4-12  shows  results  summarizing  a  series  of  dynamic  simulations  of  the  hve 
d.o.f.  human  model.  Fach  three  dimensional  plot  shows  how  many  stable  somersaults 
the  hve  d.o.f.  model  exhibited  in  simulation  as  a  function  of  the  nominal  symmetric 
shoulder  angle,  F^o,  and  the  shoulder  spring,  Ksh-  Fach  point  in  the  grid  represents 
a  separate  simulation.  A  somersault  was  considered  stable  if  it  did  not  exhibit  a 
twist  angle  of  at  least  ±7r/2  rad.  There  was  little  ambiguity  between  stable  and 
unstable  layout  somersaults.  Fither  the  twist  angle  remained  small  while  the  model 
somersaulted  or  the  twist  angle  would  grow  signihcantly  beyond  one-quarter  twist 
making  it  easy  to  distinguish  stable  and  unstable  motion.  The  plots  show  that  the 
region  in  parameter  space  that  the  linear  analysis  predicted  to  be  stable  is  also  the 
region  that  produces  the  most  stable  layout  somersaults  in  non-linear  simulation. 

The  height  of  the  plots  in  Figure  4-12  has  been  limited  to  eight  somersaults. 
Several  simulations  exceeded  this  limit.  For  selected  cases,  the  non-linear  dynamic 
simulations  indicated  that  the  maneuver  may  remain  stable  indehnitely  as  long  as 
energy  is  conserved.  Two  simulations  starting  from  an  initial  tilt  angle  of  0.01  radians 
and  zero  twist  angle  were  stopped  after  running  for  2000  simulated  seconds.  During 
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these  simulations  the  performer  produced  over  2750  stable  layout  somersaults.  The 
shoulder  springs  and  nominal  arm  angles  that  produced  these  long  running  maneuvers 
were  {Ah/i,rso}  =  {1.2, 1.4},  {1.0, 1.6}.  Other  conditions  may  also  have  produced 
stable,  long  running  maneuvers.  I  did  not  perform  an  exhaustive  search. 

4.4  How  Passive  Layout  Stability  Works 

In  the  last  section  we  showed  that  under  certain  conditions  a  hve  d.o.f.  dynamic 
model  of  a  human  can  exhibit  passively  stable  layout  somersaults.  We  used  numer¬ 
ical  searches  to  hud  combinations  of  shoulder  springs  and  nominal  arm  angles  that 
stabilized  the  layout  somersault.  This  search  process  does  not  provide  much  insight 
into  how  or  why  passive  stabilization  works.  In  this  section,  I  describe  what  I  consider 
to  be  the  fundamental  dynamic  processes  that  make  passive  layout  stability  possible. 
First,  I  look  at  the  differences  in  state  trajectories  between  the  passively  stable  case 
and  the  unstable  rigid  body  case.  The  differences  in  trajectories  between  the  two 
cases  make  evident  that  arm  motion  can  change  the  orientation  of  the  principal  axes 
in  a  way  that  stabilizes  the  layout  somersault.  We  then  look  more  closely  at  how 
this  arm  movement  arises.  To  do  this  I  focus  on  the  simplest  possible  model  that 
captures  passive  layout  stability.  I  dehne  the  important  non-dimensional  parameters 
that  govern  this  model  and  show  how  stability  depends  upon  these  parameters.  Then 
I  provide  an  intuitive  explanation  of  the  passive  dynamic  process  that  stabilizes  the 
layout  somersault. 

4.4.1  Stabilization  Via  Principal  Axes  Reorientation 

Batterman  [Batterman  68]  and  later  Frohlich  [Frohlich  79]  showed  how  an  athlete  can 
control  the  orientation  of  his  principal  axes  during  a  rotating  maneuver  using  arm 
movement.  This  control  allows  an  athlete  to  initiate  twist  in  a  somersault  without 
using  any  external  forces  or  torques.  Using  symmetry,  Yeadon  [Yeadon  84]  has  argued 
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Tilt  =  .03,  Twist  =  0  Tilt  =  .05,  Twist  =  0 


Figure  4-12:  These  plots  show  that  the  regions  (in  F^o  -  R,sh  space)  of  stable  som¬ 
ersaulting  motion  of  a  non-linear  dynamic  simulation  are  similar  to  those  regions 
predicted  to  be  stable  by  the  linear  analysis  (Figure  4-9).  The  height  of  each  grid 
point  indicates  the  number  of  stable  layout  somersaults  performed  during  a  separate 
dynamic  simulation  of  the  hve  d.o.f.  model.  The  axes  of  each  grid  indicate  the  sym¬ 
metric  arm  angle,  0.4  rad  <  F^o  <  2.2  rad,  and  the  non-dimensional  shoulder  spring, 
0.0  <  Ksh  <  1.8,  used  for  each  simulation.  The  different  plots  show  that  the  number 
of  stable  somersaults  performed  is  sensitive  to  the  initial  tilt  and  twist  angles. 
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Figure  4-13:  The  orientation  of  the  principal  axes  of  inertia  can  be  controlled  with 
arm  rotation.  The  tilt  of  the  principal  axes  of  the  hgure  on  the  right  is  a  result  of 
body  axis  tilt  due  to  airborne  reorientation  of  the  body  plus  principal  axes  tilt  due 
to  body  asymmetry. 

that  an  athlete  has  an  equal  opportunity  to  remove  twist  from  a  twisting  somersault 
using  arm  motion.  In  this  section  I  will  show  how  arm  movement  is  used  to  control 
the  orientation  of  the  principal  axes.  Then  I  will  show  how  principal  axis  orientation 
can  be  used  to  influence  rotational  maneuvers. 

The  orientation  of  a  human’s  principal  axes  relative  to  a  set  of  body  hxed  axes 
depends  upon  the  conhguration  of  the  arms  and  legs  of  the  performer.  Figure  4-13 
shows  a  human  hgure  in  two  different  body  conhgurations.  The  tilt  of  the  body  on 
the  right  is  the  result  of  the  airborne  reorientation  of  the  arms  from  the  symmetric 
conhguration  of  the  hgure  on  the  left  to  the  asymmetric  conhguration  of  the  hgure 
on  the  right.  Also  shown  are  inertially  hxed  axes  and  principal  axes  of  the  hgures. 
The  principal  axes  of  the  hgure  on  the  left  are  aligned  with  the  inertial  axes.  Notice 
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that  the  principal  axes  of  the  hgure  on  the  right  are  tilted  with  respect  to  the  inertial 
axes  and  also  inclined  relative  to  the  body  axes.  This  inclination  of  the  principal 
axes  from  the  body  hxed  axes  is  due  to  asymmetry  of  the  body  and  it  has  the  affect 
of  changing  the  ‘gain’  on  principal  axis  reorientation  due  to  arm  movement.  To  see 
how  principal  axis  orientation  can  be  used  to  control  twist  we  will  examine  a  plot  of 
rigid-body  rotational  trajectories. 

Figure  4-14  shows  the  twist  angles  plotted  as  a  function  of  tilt  angle  for  several 
rigid  body  trajectories  of  the  human  model  while  in  the  layout  conhguration.  This 
plot  is  essentially  a  close  up  view  of  the  {ci,  C3}  axes  of  the  map  shown  in  Figure  4-3. 
The  different  trajectories  correspond  to  different  energy  levels  of  rigid  body  rotational 
motion.  While  time  is  not  explicitly  shown  on  the  plot,  the  data  points  are  all  0.01  sec. 
apart.  Thus  dot  spacing  shows  that  twist  rate  increases  with  tilt  magnitude.  Arrows 
have  been  drawn  on  this  plot  to  show  the  direction  of  motion.  The  instability  of 
rotation  about  the  middle  principal  axis  is  evident  from  the  converging  and  diverging 
trajectories  near  the  origin.  These  trajectories  show  that  if  body  attitude  deviates 
slightly  from  the  origin  then  it  will  depart  the  origin  in  a  direction  that  depends  upon 
which  ‘side’  of  the  equilibrium  it  is  located. 

Tilting  the  principal  axes  using  arm  movement  has  the  effect  of  moving  the  body 
from  one  energy  trajectory  to  another.  While  an  instantaneous  change  in  tilt  can 
not  change  the  twist  angle,  it  can  change  the  twist  rate  by  putting  the  body  on  a 
trajectory  with  a  different  twist  rate.  Stabilization  of  the  layout  somersault  can  be 
viewed  as  a  process  of  selectively  adding  or  subtracting  twist  over  time  using  arm 
movement  so  that  the  trajectory  remains  in  the  vicinity  of  the  origin  of  Figure  4-14. 
Figure  4-15  shows  a  trajectory  of  a  non-linear  dynamic  simulation  of  the  passive,  hve 
d.o.f.  model  of  the  human.  This  trajectory  of  tilt  and  twist  angles  is  superposed  upon 
the  rigid  body  trajectories  of  the  same  model.  The  passively  stable  trajectory  circles 
the  origin  in  a  counter-clockwise  fashion  periodically  increasing  then  decreasing  the 
twist  and  tilt  angles.  The  counterclockwise  rotation  of  the  trajectory  means  that 
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Figure  4-14:  Tilt  angles  plotted  as  a  function  of  twist  angle  for  several  different  energy 
levels  of  rigid  body  rotation  of  a  human.  Arrows  indicate  the  direction  of  motion  with 
increasing  time.  Dots  are  all  0.01  seconds  apart.  The  arms  of  the  rigid  human  model 
were  abducted  to  a  symmetric  angle  of  1.3  radians. 
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Figure  4-15:  Tilt  angles  plotted  as  a  function  of  twist  angle  for  a  non-linear  dynamic 
simulation  of  the  five  d.o.f.  passive  model  of  the  human.  The  simulated  trajectory 
is  superposed  upon  rigid  body  trajectories  of  the  same  model.  This  simulation  was 
initialized  with  a  tilt  angle  of  0.01  rad,  zero  twist  angle  and  symmetric  arm  angle  of 
1.3  rad. 

while  the  body  twists  ‘with  the  flow’  in  the  upper  and  lower  halves  of  the  plane  it  has 
to  move  ‘upstream’  with  regard  to  the  tilt  angle.  How  is  this  accomplished? 

The  system  uses  arm  movement  to  accomplish  this  ‘upstream’  tilting  action.  In 
Figure  4-15  we  see  that  by  starting  with  an  initial  condition  of  positive  tilt  angle,  zero 
twist  angle,  and  arms  symmetric  the  body  begins  to  twist  in  the  negative  direction. 
From  this  initial  condition  the  tilt  of  the  principal  axes  apparently  decreases  despite 
the  rigid  body  tendency  to  increase  the  tilt  angle.  This  decrease  in  principal  axis  tilt  is 
a  result  of  the  arms  moving  asymmetrically  (left  arm  up,  right  arm  down)  in  response 
to  the  negative  twist  rate.  The  arm  tilting  (rotation  about  the  ‘1’  axis)  arises  from 
the  combination  of  somersaulting  (rotation  about  the  ‘2’  axis)  and  twisting  (rotation 
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Figure  4-16:  A  simple  model  that  captures  the  essential  dynamics  for  passive  layout 
stability.  The  single  arm  is  allowed  to  pivot  about  its  center.  The  arm  and  body 
centers  of  mass  are  coincident. 

about  the  ‘3’  axis.)  This  gyroscopic  affect  will  tend  to  move  the  arm  in  a  direction 
that  reduces  the  tilt  angle  of  the  principal  axes  as  long  as  the  body  is  twisting. 

This  oscillation  is  a  passive  behavior  that  emerges  from  the  mechanical  properties 
of  the  system  interacting  with  the  environmental  forces  due  to  rotation.  The  gyro¬ 
scopic  forces  which  provide  the  inherent  correction  of  this  model  depend  upon  the 
shape  and  position  of  the  arm.  To  see  how  this  passive  behavior  emerges  we  will  look 
at  a  simplihed  model  of  the  dynamics. 

4.4.2  The  Simplest  Model 

The  simplest  model  that  I  can  think  of  that  captures  the  essential  dynamics  of  passive 
layout  stability  is  shown  in  Figure  4-16  This  model  has  the  feature  that  only  asym¬ 
metric  arm  motion  is  allowed.  Recall  from  Section  4.3.2  that  symmetric  arm  motion 
could  be  ignored  when  looking  for  stability  of  the  hve  d.o.f.  model.  Also,  since  the 
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center  of  mass  of  the  arm  is  coincident  with  the  center  of  mass  of  the  body,  the 
equations  of  motion  will  be  signihcantly  simplihed  as  compared  to  the  more  complex 
human  model. 

The  linear  equations  of  motion  for  this  four  d.o.f.  model  will  show  that  the 
somersault  degree  of  freedom  again  decouples  from  the  remaining  dynamics.  Ignoring 
this  degree  of  freedom  and  assuming  no  damping  at  the  pin  joint  of  the  arm  results 
in  the  following  three  d.o.f.  linear  model. 


M  i  +  $  G  i  +  (Jh  +  K')  X  =  0 


M 


G 


K 


K' 


mil  0  mi3 

0  m22  0 

mm  0  m33 

0  gi2  0 
-gi2  0  5-23 

0  —523  0 

^11  0  km 

0  k22  0 
^13  0  km 

0  0  0 
0  0  0 

0  0  k,hl^'^ 


(4.10) 


where  the  state  vector,  x,  is  comprised  of  (in  order)  tilt,  twist,  and  arm  angle  and 
ksh  is  the  spring  between  the  arm  and  the  body.  The  model  parameters  are  given  in 
Table  4.2. 

This  dynamic  model  still  requires  us  to  dehne  the  six  inertia  parameters  of  the 
body  before  we  consider  stability.  Recall  from  Chapter  2  that  the  rotational  dynamics 
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Table  4.2 

Simplest  Model  Parameters 


Parameter 

Dehnition 

mil 

Ihi  +  III 

Ik 

m22 

lbs  +  II3 

mss 

Ik 

§12 

(Ibi  +  III)  —  {lbs  +  II2)  +  {lbs  +  II3) 

§23 

—Ik  +  II2  —  II3 

kii 

(lbs  +  II2)  —  {lbs  +  II3) 

kis 

II2  —  II3 

^22 

-{Ibi  +  Ik)  +  {Ib2  +  Ik) 

ks3 

Ik  —  II3 

of  a  rigid  body  depend  only  upon  two  non-dimensional  inertia  parameters.  Following 
the  spirit  of  this  result,  in  the  next  section  we  reduce  the  model  to  non-dimensional 
form  in  order  further  simplify  the  analysis. 

Non-Dimensional  Linear  Equations  of  Motion 

I  now  present  a  non-dimensional  version  of  the  linear  dynamic  equations  of  motion  of 
the  model  shown  in  Figure  4-16  for  the  case  of  steady  somersaulting  about  the  ‘2’  axis. 
The  purpose  of  this  step  is  to  discover  the  important  non-dimensional  parameters  of 
this  model.  The  states  are  made  non-dimensional  through  scaling  by  the  nominal 
somersault  rate.  Casting  the  equations  in  terms  of  non-dimensional  parameters  is 
accomplished  by  dividing  each  equation  by  the  corresponding  inertia  term  on  the 
diagonal  of  the  mass  matrix.  A  few  more  simple  algebraic  manipulations  will  bring 
the  system  of  equations  into  the  following  form.  (In  this  form  these  equations  do  not 
reflect  the  symmetry  that  we  expect  of  a  conservative  gyric  system.  However,  they 
can  easily  be  made  to  ht  this  canonical  form.) 


$2  (MS  +  Gi  +  Kz)  =  0 


(4.11) 
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M 


1  0  Irel 

0  1  0 

1  0  1 


G 


0  1-h 

(1  -  h)  0 

0  l-kh 


0 

-{l-kli){l-k3)Irel 
1  —  ki 

0 


K 


k\  0  kl\  Irel 

Oh  0 

kli  0  kli  +  Ksh 


The  five  independent  non-dimensional  parameters  are  defined  as  follows: 


ki 

ks 


Irel 

kh 

Ksh 


{Ib2  +  II2)  —  {lbs  +  II3) 
Ibi  +  III 

{Ib2  +  II2)  -  {Ibi  +  III) 

lbs  +  II3 

Ik 

Ibi  +  III 
II2  —  II3 
III 

ksh 

//l$2 


Equation  4.11  depends  upon  only  five  non-dimensional  parameters.  The  param¬ 
eters  ki  and  k3  are  the  principal  inertia  ratios  for  an  equivalent  rigid  body  system. 
They  describe  the  inertia  properties  of  the  system  that  would  result  if  the  arm  were 
hxed  in  place.  Notice  that  if  the  arm  degree  of  freedom  were  eliminated  from  Equa¬ 
tion  4.11  then  the  remaining  equations  are  identical  to  the  rigid  body  equations  (2.11). 
Irel  dehnes  how  big  the  inertia  of  the  arm  is  relative  to  the  inertia  of  the  rigid  body.  It 
is  important  in  determining  how  much  the  body  moves  in  response  to  arm  movement. 
The  next  term,  kk  is  an  inertia  ratio  analogous  to  ki  involving  only  arm  parameters. 


CHAPTER  4.  PASSIVELY  STABLE  LAYO  UT  SOMERS  A  ULTS 


102 


Table  4.3 

Furrer’s  Non-dimensional  Parameters 


Parameter 

Value 

ki 

0.65 

ks 

-0.87 

^rel 

0.06 

kh 

-0.99 

It  is  important  in  determining  the  tendency  of  the  arm  to  tilt  in  response  to  twist 
rate. 

We  use  data  from  the  human  model  (Table  4.1)  to  derive  sample  values  for  these 
dimensionless  parameters.  These  parameter  values  are  shown  in  Table  4.3.  These 
parameters  are  computed  assuming  the  arms  are  held  straight  out  to  the  side  of  the 
body.  To  compute  the  arm  inertia  of  the  model  in  Figure  4-16  I  let  III  =  II3  = 
2.0  X  0.3849kg  the  sum  of  the  arm  inertias  about  their  shoulder  axes  in  the  given 
conhguration. 

With  these  numbers  in  mind,  we  can  inspect  the  individual  terms  of  the  matrices 
of  Equation  4.11  to  help  us  develop  an  understanding  of  the  relevant  dynamic  forces. 
Without  arm  movement  the  system  reduces  to  the  rigid  body  equations  of  motion 
which  are  unstable  for  the  set  of  parameters  in  Table  4.3.  For  this  reason  I  will 
focus  on  the  arm  related  terms  only.  First  consider  the  3-2  term  of  G.  This  term 
is  nearly  equal  to  2.0.  If  the  twist  rate  of  the  model  is  negative  then  this  term  will 
contribute  to  a  positive  acceleration  of  the  arm  angle  which  is  precisely  the  corrective 
tendency  that  we  observed  in  the  plot  of  Figure  4-15.  Additional  arm  accelerations 
come  from  the  centrifugal  forces  associated  with  somersault  rotation.  These  forces 
are  represented  by  the  3-1  and  3-3  terms  in  the  stiffness  matrix,  K .  If  K^h  =  0  then 
the  3-3  term  of  the  stiffness  matrix  is  negative.  This  reflects  the  fact  that  if  the 
arm  is  deviated  from  its  equilibrium  position  then  it  will  experience  centrifugal  forces 
from  somersault  rotation  that  tend  to  drive  the  arm  further  away  from  zero.  This 
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is  because  the  centrifugal  forces  on  either  side  of  the  arm  are  opposite  one  another 
due  to  asymmetry  and  the  forces  increase  with  deflection  of  the  arm.  Without  a 
shoulder  spring  this  term  could  overwhelm  the  helpful  gyroscopic  accelerations  from 
the  3-2  term  in  G  to  destabilize  the  system.  However,  the  shoulder  spring  offsets 
this  destabilizing  force.  In  fact,  a  rule  of  thumb  for  choosing  the  shoulder  spring 
constant  to  stabilize  the  passive  system  is  to  choose  the  spring  that  exactly  cancels 
the  destabilizing  term  due  to  centrifugal  forces. 

Performing  a  root  locus  search  for  stabilizing  spring  constants  of  this  simple  3 
d.o.f.  system  results  in  the  following  range  of  stabilizing  shoulder  springs: 

0.98  <  Ksh  <  1.08 

which  approximately  cancels  kli  =  —  1. 

This  rule  for  picking  stabilizing  shoulder  springs  generalizes  to  the  more  compli¬ 
cated  model  of  the  human.  The  most  important  difference  in  the  human  model  is  that 
the  nominal  symmetric  angle  of  the  arms  can  range  from  0.0  radians  (hands  next  to 
legs)  to  TT  radians  (hands  overhead).  The  non-dimensional  shoulder  spring  constants 
that  stabilize  the  human  model  are  plotted  as  a  function  of  nominal  symmetric  arm 
angle  in  Figure  4-17  (repeated  from  Figure  4-10).  Also  shown  in  this  plot  is  the  value 
of  Ksh  that  would  exactly  cancel  the  spring  term  due  to  centrifugal  forces  on  the  arm. 
Again,  the  region  of  stabilizing  spring  constants  seems  to  be  dehned  by  the  negative 
of  this  centrifugal  force  term.  Finally,  it  is  interesting  to  note  that  for  a  nominal 
symmetric  arm  angle  of  7r/2  the  non-dimensional  stabilizing  shoulder  spring  constant 
is  approximately  Kgh  =  1.0. 

Stability  Diagrams  for  the  Simplest  Model 

The  stability  analyses  performed  thus  far  have  all  assumed  that  we  have  known  values 
of  the  body  inertia  parameters.  The  small  number  of  non-dimensional  parameters  of 
Fquation  4.11  make  it  possible  to  study  how  linear  stability  depends  upon  these 


Figure  4-17:  Plot  of  the  stabilizing  K^h  for  a  the  3-by-3  subsystem  of  the  human 
model.  The  solid  curve  shows  the  value  of  K^h  that  will  exactly  cancel  the  spring 
term  due  to  centrifugal  forces  on  the  arm. 
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inertia  parameters  as  well.  In  this  section  I  present  two  stability  diagrams  that  show 
how  linear  stability  depends  upon  the  parameters  Reh  s-iid  Ksh-  The  stability 

diagrams  consist  of  a  sequence  of  plots  of  the  {ki^ks}  parameter  space.  The  darkened 
regions  of  these  plots  represent  values  of  ki  and  ks  that  are  stable  for  the  specihed 
values  of  Rei  and  Ksh- 

Figure  4-18  shows  that  the  value  of  Kgh  that  stabilizes  the  largest  region  of  the 
unstable  fourth  quadrant  of  the  {ki^ks}  parameter  space  is  Kgh  ~  1-1.  Figure  4-19 
shows  that  larger  the  value  of  Rei  the  easier  it  is  to  stabilize  the  bottom  two  quadrants 
of  the  {ki^kR  parameter  space.  Large  Rei  corresponds  to  relatively  larger  arms. 

The  parameter  kli  is  a  shape  parameter  for  the  arm  that  describes  the  tendency 
of  the  arm  to  tilt  in  response  to  twisting  motion.  Figures  4-9,  4-10,  and  4-11  indicate 
that  the  effectiveness  of  the  arm  to  stabilize  the  layout  somersault  decreases  as  the 
arms  are  moved  away  from  F^o  =  7r/2  (hands  held  straight  out  to  the  side.)  I  believe 
this  change  in  effectiveness  is  due  the  change  in  the  parameter  kli.  To  check  this,  I 
allowed  the  shape  of  the  arm  in  Figure  4-16  to  reflect  the  symmetric  arm  angle  F^o- 
Then  I  searched  for  the  Kgh  that  would  stabilize  the  simple  three  d.o.f.  model  with 
hxed  body  parameters  representative  of  those  for  a  human.  The  stable  F^o  -  Kgh 
conhgurations  are  shown  in  Figure  4-20.  The  shape  of  this  plot  is  suggestive  of  the 
shape  in  the  stability  diagrams  for  the  hve  d.o.f.  model.  An  important  difference 
between  the  diagrams  for  the  simple  three  d.o.f.  model  and  the  hve  d.o.f.  model  is 
that  the  center  of  mass  of  the  arm  for  the  three  d.o.f.  model  is  coincident  with  the 
pin  joint.  I  think  this  may  lead  to  the  extra  stable  conditions  in  Figure  4-20. 

4.5  Summary 

In  this  chapter  I  presented  results  from  a  study  of  the  stability  of  several  simple 
models  of  the  layout  somersault.  A  rigid  body  analysis  of  the  layout  somersault 
predicts  that  the  maneuver  is  inherently  unstable.  A  linear  analysis  of  a  passive  three- 
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Figure  4-18:  These  plots  show  that  1.0  <  K^h  <  1-1  produces  the  largest  region  of 
stable  conhgurations  for  rotation  about  the  middle  principal  axis  for  the  “simplest 
model”.  Each  dark  point  represents  a  stable  conhguration  in  ki  -  space  for  rotation 
about  a  principal  axis,  kli  =  — 1,/re/  =  0.1  Each  plot  is  for  a  different  value  of  Ksh- 
Ksh  =  0.0,  0.5,  0.8,  0.9, 1.0, 1.1, 1.2, 1.5  and  2.0  for  plots  positioned  from  left  to  right 
and  down  the  page. 
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Figure  4-19:  These  plots  show  that  stabilization  gets  easier  with  larger  arms.  Each 
dark  point  represents  a  stable  conhguration  in  ki  -  space  for  rotation  about  a 
principal  axis,  kli  =  — l,Ah/i  =  1.0.  Each  plot  is  for  a  different  value  of  Rei-  Rei  = 
0.01,  0.02,  0.05,  0.1,  0.2,  0.3,  0.4,  0.5  and  0.6  for  plots  positioned  from  left  to  right  and 
down  the  page. 
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Figure  4-20:  This  plot  shows  that  the  stability  of  the  “simplest  model”  can  reflect 
the  dependence  on  arm  position  that  we  saw  in  the  hve  d.o.f.  model  (Figure  4-9.) 
Fach  dark  point  represents  a  stable  conhguration  of  the  simple  three  d.o.f.  model. 
The  shape  of  the  arm  of  the  “simplest  model”  is  changed  with  the  value  of  F^o  to 
resemble  the  nominal  symmetric  arm  angles  of  the  hve  d.o.f.  model. 
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body  human  model  of  the  layout  somersault  predicts  that  neutral  stability  can  be 
achieved  for  a  range  of  nominal  symmetric  arm  angles  and  shoulder  spring  constants. 
The  model  is  considered  passive  because  the  joint  and  body  motion  is  the  result 
of  the  passive  dynamic  interaction  of  the  bodies  and  the  environmental  forces.  No 
active  control  or  sensing  is  used  in  the  model.  Non-linear  dynamic  simulation  of  the 
passive  three-body  model  shows  that  the  most  stable  conhgurations  conform  to  the 
linear  analysis  results.  The  number  of  stable  somersaults  that  the  three  body  model 
performed  in  dynamic  simulation  before  exhibiting  a  twist  instability  is  sensitive  to 
initial  conditions.  However,  for  selected  cases,  non-linear  dynamic  simulations  of  the 
layout  somersault  indicated  that  the  maneuver  may  remain  stable  indehnitely  as  long 
as  energy  is  conserved. 

In  order  to  provide  a  more  intuitive  understanding  of  how  passive  stabilization  of 
the  layout  somersault  works,  we  considered  the  impact  of  principal  axis  orientation  on 
layout  stability.  I  showed  how  the  natural  tendency  of  the  arms  to  tilt  in  response  to 
twisting  movement  of  the  body  provides  a  built-in  correction  to  the  twist  instability. 
The  arm  tilt  forces  the  principal  axes  of  the  system  to  move  in  a  direction  that 
compensates  for  tilt  and  twist  errors.  This  built-in  correction  eliminates  the  divergent 
tendency  of  the  system  as  long  as  the  compliance  of  the  shoulders  cancels  the  unstable 
centrifugal  forces  on  the  arms.  This  effect  is  clarihed  by  studying  a  very  simple  two 
body  model  of  the  layout  somersault.  This  simple  model  also  allows  us  to  derive  hve 
non-dimensional  parameters  that  are  important  in  studying  passive  layout  stability. 
I  present  stability  diagrams  that  show  how  stability  of  the  simplest  model  depends 
upon  these  parameters. 


Chapter  5 


Layout  Somersault  Experiments 

5.1  Introduction 

To  find  out  whether  or  not  passively  stable  layout  somersaults  are  physically  possi¬ 
ble,  we  built  a  somersaulting  ’doll’  that  can  be  tested  in  the  laboratory.  The  doll  has 
springy  shoulders  that  allow  arm  movement  in  the  frontal  plane.  The  non-dimensional 
inertia  parameters  of  the  doll  show  that  it  is  dynamically  similar  to  the  hve  d.o.f. 
model  of  the  human  performer  studied  in  Chapter  4.  The  goals  of  the  experiments 
were  1)  to  determine  if  passive  layout  stability  was  physically  possible,  and,  2)  em¬ 
pirically  determine  the  shoulder  spring  constants  that  best  stabilized  the  motion. 

To  initialize  the  somersaults,  we  built  a  mechanical  launching  device  that  throws 
the  doll  into  the  air  with  angular  rotation  about  the  somersault  axis.  The  launcher 
helped  us  to  achieve  consistency  in  the  experimental  conditions  and  helped  to  mini¬ 
mize  the  human  influence  of  the  experiment.  The  operation  of  the  launcher  and  doll 
is  similar  in  appearance  to  a  gymnast  swinging  around  a  horizontal  bar  then  releasing 
the  bar  to  perform  a  multiple  somersaulting  dismount. 

In  this  chapter  I  describe  the  human-like  doll  and  launching  device  used  during 
experiments.  Using  models  from  the  previous  chapter,  I  perform  an  analytic  stabil¬ 
ity  analysis  of  the  doll  to  determine  under  what  conditions  the  doll  is  theoretically 


no 
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stable  in  somersault  rotation.  Then,  I  describe  the  somersault  experiments  in  which 
we  repetitively  launched  the  doll  with  different  shoulder  springs  to  determine  which 
springs  best  stabilized  the  doll.  During  experiments  we  found  that  the  doll  could  con¬ 
sistently  perform  at  least  three  and  one  half  stable  layout  somersaults  during  flight. 
The  doll  may  have  performed  more  layout  somersaults  if  we  had  been  able  to  observe 
longer  flight  times.  The  laboratory  ceiling  height  limited  our  opportunity  to  observe 
more  somersaults.  Also,  we  found  good  agreement  between  theory  and  experiment  in 
finding  the  shoulder  springs  that  best  stabilized  the  maneuver. 

5.2  Experimental  Apparatus 

In  this  section  I  describe  the  human-like  doll  used  during  experiments.  I  compare  three 
non-dimensional  parameters  of  the  doll  and  the  human  model  studied  in  Chapter  4 
to  show  that  the  two  models  are  dynamically  similar.  I  also  describe  the  launching 
device  used  to  initialize  the  somersaults  during  experiments. 

5.2.1  Human-Like  Doll 

I  built  a  mechanical  version  of  the  five  d.o.f.  human  model  studied  in  Chapter  4.  This 
doll  is  a  fifteen  inch  tall,  scaled  version  of  the  human  performer  whose  anthropometric 
data  were  used  in  the  stability  analysis  of  Chapter  4.  A  photograph  of  this  doll  is 
shown  in  Figure  5-1.  The  body  is  made  of  wood.  The  shoulder  joint  axes  are  3/16 
in.  diameter  steel  shafts  that  are  fixed  to  the  body  via  aluminum  brackets  screwed  to 
the  front  and  back  of  the  doll  body.  The  arms  are  made  of  1/4"  threaded  steel  rod. 
The  arms  are  attached  to  the  shoulder  shaft  via  an  aluminum  block  fitted  with  ball 
bearings  that  allow  free  shoulder  rotation  in  the  frontal  plane.  Three  shoulder  pulleys 
with  1/2  in.,  3/8  in.,  and  5/16  in.  radii  are  mounted  to  the  aluminum  block  to  be 
concentric  with  the  shoulder  shaft.  Linear  springs  attached  between  the  body  and 
the  shoulder  pulleys  act  as  torsional  springs  about  the  shoulder  shaft.  The  different 


CHAPTER  5.  LAYO  UT  SOMERS  A  ULT  EXPERIMENTS 


112 


Figure  5-1:  Photograph  of  the  mechanical  doll  used  for  experiments  with  layout 
somersaults. 

radii  pulleys  allow  a  single  linear  spring  to  be  used  as  three  different  torsional  springs. 
Two  sets  of  springs  were  used  in  opposition  to  each  shoulder.  The  shoulder  springs 
could  be  pre-tensioned  so  that  the  arm  had  any  desired  equilibrium  angle.  The  doll 
had  a  1/4  in.  diameter  vertical  shaft  attached  to  its  head.  This  shaft  is  used  as  a 
handle  for  the  mechanical  launcher  to  grasp.  One  30^'  threaded  weight  was  added  to 
each  arm.  Moving  the  position  of  this  weight  on  each  arm  allowed  us  to  change  the 
inertia  of  the  arm  without  changing  the  total  mass  of  the  body. 

The  physical  parameters  of  the  hve  d.o.f.  doll  without  added  arm  weights  and 
with  the  arm  weight  located  at  a  radius  of  0.043m  from  the  shoulder  center  of  rotation 
are  given  in  Table  5.2.1.  The  arm  inertias  are  given  in  an  axis  system  located  at  the 
center  of  mass  of  the  arm  and  parallel  to  the  body  axis  system  when  the  arm  is  held 
straight  down  at  the  side,  $so  =  0.0.  Data  for  this  chapter  are  all  for  the  doll  with 
arm  weight  at  a  radius  of  0.043  m  from  the  shoulder. 
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Table  5.1 

Doll  Body  Parameters 


Parameter 

No  weight 

30g  weight 
r  =  0.043  m 

nib 

0.466 

kg 

0.466 

kg 

nil 

0.0275 

kg 

0.0575 

kg 

IW 

0.00411 

kg  rrC 

0.00411 

kg  rrC 

0.00391 

kg  rrC 

0.00391 

kg  rrC 

lbs 

0.000354 

kg  rrC 

0.000354 

kg  rrC 

Ik 

0.0000559 

kg  rrC 

0.0000789 

kg  rrC 

Ik 

0.0000559 

kg  rrC 

0.0000789 

kg  rrC 

Ik 

0.000001 

kg  rrC 

0.000001 

kg  rrC 

rbi 

0.0 

m 

0.0 

m 

rbs 

0.0413 

m 

0.0413 

m 

rbs 

0.0664 

m 

0.0664 

m 

rk 

0.0 

m 

0.0 

m 

rk 

0.0 

m 

0.0 

m 

rk 

0.0832 

m 

0.0572 

m 

Figure  5.2.1  shows  a  plot  of  the  rigid  body  inertia  ratios  of  the  doll  and  of  the 
human  performer  for  symmetric  arm  angles  between  0.0  and  3.0  radians.  These  are 
superposed  upon  curves  that  show  how  quickly  the  unstable  somersault  mode  grows 
for  a  rigid  body.  Comparison  of  the  inertia  ratios  shows  that  the  doll  has  nearly 
the  same  rigid  body  inertia  ratios  as  the  human  model.  Figure  5-3  shows  a  plot  of 
the  non-dimensional  parameter  foi"  both  the  doll  and  the  human  model  for  arm 
angles  between  0.0  and  3.0  radians.  Comparing  these  curves  shows  that  the  doll 
has  relatively  more  massive  arms  than  the  human  which  makes  stabilization  slightly 
easier. 

5.2.2  Launching  Device 

We  built  a  launching  device  so  that  we  could  initialize  the  somersaults  with  consis¬ 
tency  and  with  minimum  direct  human  influence.  A  photograph  of  the  launcher  is 
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Figure  5-3:  Plot  of  Rei  as  a  function  of  nominal  symmetric  arm  angle,  F^o,  for  the 
doll  (dotted  line)  and  the  human  model. 


shown  in  Figure  5-4.  The  launcher  functioned  by  accelerating  the  doll  in  rotation 
around  a  horizontal  bar  then  releasing  the  doll  at  a  hxed  angle  around  the  bar.  The 
launcher  consisted  of  three  components:  1)  a  horizontal  bar,  2)  a  bearing  supported 
toggle  clamp  that  was  free  to  rotate  around  the  bar,  and  3)  a  torsion  spring  used  to 
accelerate  the  clamp  relative  to  the  bar. 

To  initialize  the  somersaults,  the  steel  rod  hxed  to  the  head  of  the  doll  was  clasped 
in  the  toggle  clamp.  The  doll  and  toggle  clamp  would  then  be  wound  around  the  hor¬ 
izontal  bar  to  tension  the  launcher  spring.  When  the  doll  and  toggle  clamp  were 
released  they  would  accelerate  around  the  horizontal  bar.  When  the  toggle  clamp 
lever  hit  a  trigger  bar,  the  toggle  clamp  popped  open,  releasing  the  doll  into  a  somer¬ 
saulting,  parabolic  trajectory.  The  position  of  the  trigger  bar  around  the  horizontal 
bar  could  be  changed  so  that  the  toggle  clamp  released  the  doll  at  an  angle  that  pro¬ 
duced  the  desired  flight  trajectory.  The  toggle  clamp  had  adjustments  that  allowed 
the  doll  to  be  launched  with  initial  tilt  or  twist  angles,  however,  for  all  experiments 
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Figure  5-4:  Photograph  of  the  mechanical  launcher  used  to  initialize  the  doll  somer¬ 
saults  in  experiments. 
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Figure  5-5:  The  points  in  this  plot  indicate  the  values  of  K^h  that  theoretically  will 
stabilize  layout  somersaults  in  the  doll.  The  cross-hairs  indicate  the  spring  constant 
that  best  stabilized  the  doll  during  experiments. 

we  set  the  initial  tilt  and  twist  angles  to  be  approximately  zero. 

5.3  Theoretical  Predictions 

Figure  5-5  shows  the  theoretically  determined  stabilizing  values  of  K^h  as  a  function  of 
symmetric  arm  angle  for  the  doll.  Also  shown  in  this  hgure  is  the  spring  constant  that 
best  stabilized  the  doll  with  a  nominal  arm  angle  of  F^o  ~  1-57  during  experiments. 
Note  that  this  plot  of  the  stabilizing  K^h  is  very  similar  to  those  for  the  human  in 
Figure  4-9.  The  most  notable  difference  is  the  larger  range  of  stabilizing  K^h  for  the 
doll.  This  is  due  to  the  larger  Rei  of  the  doll. 
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5.4  Description  of  Experiments 

The  experiments  consisted  of  a  sequence  of  launches  of  the  doll  to  determine  which 
shoulder  springs  stabilized  the  doll  in  somersault  rotation.  For  each  different  shoulder 
spring  used  on  the  doll,  we  launched  the  doll  twenty  times  using  approximately  the 
same  launch  spring  pre-tension.  The  toggle  clamp  released  the  doll  from  the  bar  at 
a  somersault  angle  of  approximately  rad.  We  video  taped  each  experiment.  Using 
the  video  tape  for  analysis,  for  each  launch  we  recorded  the  number  of  somersaults 
performed  before  the  doll  exhibited  a  twist  angle  of  approximately  |  (one  quarter 
twist).  If  during  a  dismount  the  doll  never  exhibited  this  large  twist  angle  then  the 
dismount  was  considered  to  be  completely  stable.  For  comparison,  we  also  performed 
twenty  launches  with  the  arms  of  the  doll  hxed  rigidly  in  place. 

Recall  from  Chapter  4  that  the  nominal  symmetric  arm  angle  does  affect  the 
value  of  the  stabilizing  spring  constant.  During  experiments  we  attempted  to  keep 
the  nominal  symmetric  arm  angle  approximately  equal  to  1.6  radians  by  initializing 
the  somersault  with  the  appropriate  arm  angle  and  shoulder  spring  pretension.  Large 
centrifugal  forces  on  the  arms  during  acceleration  around  the  bar  necessitated  that  we 
devise  a  method  for  hxing  the  arm  position  until  the  moment  of  release  from  the  bar. 
In  order  to  initialize  the  somersault  with  the  desired  arm  angle  we  propped  up  the 
arms  with  thin  steel  rods  held  between  the  arms  and  small  indentures  near  the  feet  of 
the  doll.  The  rods  were  held  in  place  by  the  shoulder  spring  pretension.  As  the  doll 
rotated  about  the  bar,  centrifugal  forces  pushed  the  arms  towards  the  feet  helping  to 
hold  the  steel  rods  and  arms  in  place.  Once  free  of  the  bar,  the  doll  began  rotating 
about  its  center  of  mass  which  tended  to  lift  the  arms  off  the  steel  rods  allowing  the 
rods  to  fall  clear  of  the  doll.  The  nominal  arm  angle  was  difficult  to  set  precisely  as 
it  depended  not  only  on  initial  position  of  the  arms  but  also  on  the  pretension  of  the 
springs  necessary  to  hold  the  arms  at  this  angle  during  rotation. 

Figure  5-6  shows  the  beginning  of  two  separate  experiments  with  the  doll. 
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Figure  5-6:  The  photo  on  the  left  show  the  rigid  doll  with  arms  clamped  in  place. 
This  doll  exhibits  the  twist  instability.  The  photo  on  the  right  shows  the  doll  with 
flexible  arms  performing  a  stable  layout  somersault. 
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5.5  Experimental  Results 

There  was  little  ambiguity  between  stable  and  unstable  dismounts.  Either  the  twist 
angle  remained  small  while  the  doll  somersaulted  or  it  would  grow  signihcantly  beyond 
one-quarter  twist  making  it  easy  to  distinguish  stable  and  unstable  cases.  The  rotation 
rate  available  from  our  launcher  and  laboratory  ceiling  height  of  4.3  m  limited  the 
doll  to  approximately  three  and  one  half  complete  somersaults  from  release  to  landing 
in  the  hands  of  our  human  ’catcher’.  Figure  5-7  shows  the  average  number  of  stable 
somersaults  performed  as  a  function  of  the  six  different  shoulder  spring  constants 
tested.  The  average  is  computed  over  twenty  sequential  launches  for  a  single  shoulder 
spring.  The  standard  deviation  for  each  condition  is  shown  with  error  bars.  The 
average  somersault  rate  was  16.6  rad/sec,  and  Rh  =  3.021  •  10“‘^A^m^. 

In  addition  to  recording  the  number  of  stable  somersaults  exhibited  during  each 
launch  we  recorded  whether  or  not  the  doll  still  appeared  stable  at  the  end  of  the 
maneuver  (did  not  exhibit  a  quarter  twist).  If  the  doll  had  not  exhibited  a  quarter 
twist  by  the  end  of  the  maneuver  then  it  was  considered  to  be  a  completely  stable 
dismount.  The  percentage  of  completely  stable  dismounts  as  a  function  of  Ah/i  is 
shown  in  Figure  5-8. 

Of  the  six  conditions  tested,  (hve  springs  plus  rigidly  hxed  arms)  one  shoulder 
spring  value  clearly  outperformed  the  rest.  The  shoulder  spring,  Ah/i  =  0.82,  per¬ 
formed  18/20  dismounts  without  exhibiting  the  twist  instability.  This  high  percent¬ 
age  of  completely  stable  dismounts  led  to  the  very  small  variance  in  Figure  5-7.  If  we 
could  have  observed  more  somersaults,  the  doll  may  have  exhibited  a  variance  more 
in  accordance  with  the  remaining  data.  The  variance  in  number  of  stable  somersaults 
appears  to  increase  with  softer  shoulder  springs.  While  the  doll  became  more  erratic 
under  these  conditions,  it  never-the-less  occasionally  performed  a  completely  stable 
dismount.  At  the  other  extreme,  the  doll  with  hxed  arms  reliably  exhibited  the  twist 
instability  at  about  one  complete  somersault  from  release. 

The  experimental  results  contained  in  Figures  5-7  and  5-8  are  in  agreement  with 
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Figure  5-7:  Plot  of  the  average  number  of  stable  somersaults  performed  during  twenty 
launches  of  the  doll  as  a  function  of  the  shoulder  spring  constant,  Ksh-  The  standard 
deviation  is  indicated  by  error  bars.  The  data  for  the  rigidly  hxed  arms  is  indicated 
by  the  inhnite  value  of  Ksh-  The  nominal  arm  angle  was  approximately  1.6  radians. 
The  best  value  of  Kgh  =  0.82  falls  withing  the  range  of  stability  predicted  by  theory 
as  indicated  in  Figure  5-5. 
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Figure  5-8:  Plot  of  the  percentage  (of  twenty)  of  completely  stable  dismounts  as  a 
function  of  the  non-dimensional  shoulder  spring  constant,  Ksh-  The  two  test  condi¬ 
tions  on  the  right  side  of  the  plot  exhibited  zero  completely  stable  somersaults.  The 
best  value  of  K^h  =  0.82  falls  withing  the  range  of  stability  predicted  by  theory  as 
indicated  in  Figure  5-5. 
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the  theoretical  results  for  stable  layout  somersaults.  Figure  5-5  shows  that  for  a 
Fso  ~  1.57,  the  stabilizing  shoulder  spring  should  have  values  in  the  range  0.7  < 
Ksh  A  1-4.  The  best  experimental  value  for  the  doll,  K^h  ~  0.8  falls  in  this  range 
with  performance  falling  off  on  either  side.  Furthermore,  this  best  value  also  falls 
in  the  range  of  stabilizing  shoulder  springs  for  the  human  model  data  as  shown  in 
Figure  4-9. 

5.6  Summary 

In  this  chapter  I  presented  the  results  of  layout  somersault  experiments  on  a  mechani¬ 
cal  human-like  doll.  Comparison  of  the  non-dimensional  model  parameters  of  the  doll 
and  of  a  human  indicate  that  the  doll  is  dynamically  similar  to  the  human  although 
slightly  easier  to  stabilize  due  to  relatively  more  massive  arms.  The  experiments 
demonstrate  that  the  doll  can  consistently  perform  at  least  three  and  one  half  sta¬ 
ble  layout  somersaults.  The  consistency  of  the  dismounts  depends  strongly  upon  the 
value  of  the  shoulder  spring.  The  best  shoulder  spring  exhibited  18/20  completely  sta¬ 
ble  dismounts,  somersaulting  dismounts  without  any  evident  twist  instability.  While 
softer  shoulder  springs  make  the  dismount  more  erratic,  they  never-the-less  allow  the 
doll  to  occasionally  perform  a  completely  stable  dismount  and  regularly  perform  two 
stable  layout  somersaults.  On  the  other  hand,  stiff  springs  make  the  doll  more  reliable 
but  less  stable.  The  best  value  of  the  non-dimensional  shoulder  spring  constant  was 
Ksh  =  0.82  which  is  within  the  region  of  stability  predicted  by  the  linear  stability 
theory. 


Chapter  6 


Twisting  Somersaults 

6.1  Introduction 

The  twisting  somersault  is  a  maneuver  in  which  the  performer  simultaneously  ro¬ 
tates  about  the  somersault  and  twist  axes  of  the  body.  Multiple  twisting,  multiple 
somersaulting  maneuvers  are  among  the  most  exciting  and  complex  aerial  maneuvers 
performed  by  gymnasts  and  other  athletes.  Unlike  the  pure  somersault,  the  twisting 
somersault  must  include  non-linear  rotational  coupling  between  the  different  body 
axes.  One  effect  of  this  non-linearity  is  that  the  effect  of  the  performer’s  control 
actions  will  change  during  a  maneuver  making  cause  and  effect  relationships  more 
complex  than  in  the  pure  somersault.  Navigation  and  feedback  control  of  twisting 
maneuvers  are  challenging  tasks.  Does  the  accurate,  reliable  performance  of  this 
maneuver  necessitate  a  feedback  control  strategy? 

In  this  chapter  I  briefly  discuss  the  open  loop  control  of  twisting  somersaults.  Open 
loop  control  means  the  performer’s  control  actions  are  simply  replayed  from  memory. 
We  would  like  to  know  if  twisting  somersaults  could  exhibit  passive  dynamic  stability 
when  performed  with  an  open  loop  control  strategy.  We  discuss  one  test  we  performed 
to  look  for  evidence  of  passive  dynamic  stability  in  the  twisting  somersault. 

I  present  results  from  the  non-linear  dynamic  simulation  of  a  1  1/2  twisting  front 
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somersault.  This  maneuver  requires  that  the  performer  execute  a  sequence  of  limb 
motions  during  flight.  We  found  that  when  a  simulated  performer  used  a  prescribed 
set  of  motions  for  executing  this  maneuver,  the  landing  attitude  of  the  performer  was 
sensitive  to  initial  conditions  leading  to  poor  landing  attitudes.  However,  when  the 
performer’s  control  movements  were  compliant,  the  reliability  of  the  landing  attitude 
was  signihcantly  improved.  There  appears  to  be  an  optimal  choice  for  the  performer’s 
compliance  that  leads  to  the  most  reliable  landings. 

I  also  discuss  the  control  of  twisting  somersaults  by  the  3D  Biped  robot.  We 
programmed  a  simulated  3D  Biped  robot  to  perform  a  front  somersault  with  half 
twist.  We  found  that  in  order  to  make  the  robot  maneuver  look  like  a  front  somersault 
with  twist  as  performed  by  a  human,  we  had  to  add  weight  to  the  robot  to  make  its 
moments  of  inertia  more  like  those  of  a  human.  We  also  had  to  use  stronger  actuators 
than  available  for  the  physical  robot.  Our  experiments  to  make  the  real  3D  Biped 
robot  perform  the  maneuver  in  the  laboratory  were  unsuccessful.  The  physical  robot 
actuators  had  insufficient  actuator  power  to  perform  the  maneuver. 

6.2  The  Tilt  of  Twisting  Somersaults 

Frohlich  [Frohlich  80]  described  two  techniques  for  performing  twisting  somersaults. 
In  the  torque  twist  the  athlete  derives  rotation  about  the  somersault  and  twist  axes 
from  external  forces  as  he  or  she  leaves  the  ground,  diving  board  or  other  apparatus. 
In  the  torque-free  twist  with  angular  momentum  the  athlete  initiates  twist  from 
an  airborne  somersault  with  an  asymmetric  movement  of  the  limbs  (Figure  f-f).  The 
net  effect  of  either  technique  is  to  tilt  the  principal  axes  of  the  body  relative  to  the 
angular  momentum  vector.  A  map  of  rotation,  such  as  that  in  Figure  2-7,  shows  that 
even  a  small  amount  of  tilt  of  the  principal  axes  from  the  layout  somersault  position 
will  result  in  twisting.  (This  is  what  makes  layout  somersaults  challenging.)  The 
greater  the  tilt  angle,  the  greater  will  be  the  twist  rate. 
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6.3  One  and  One  Half  Twisting  Front  Somersault 

Yeadon  [Yeadon  84]  discusses  a  “torque-free”  twist  technique  based  on  the  hula  move¬ 
ment  for  initiating  twist  from  a  piked  front  somersault.  The  hula  movement  involves 
a  swiveling  of  the  hips  not  unlike  that  required  to  swing  a  hula  hoop  about  the  hips. 
A  quarter  cycle  of  hula  movement  performed  during  a  pike  front  somersault  will  tilt 
the  performer’s  principal  axes  relative  to  the  angular  momentum  vector.  This  effect 
is  increased  if  the  arms  are  held  in  an  abducted  position  during  the  movement.  After 
the  hula  movement  the  performer  extends  from  the  pike  and  the  arms  are  adducted 
to  decrease  the  inertia  about  the  twist  axis.  The  extension  from  the  hula  movement 
should  happen  between  the  1/4  and  3/4  twist  positions  to  maximize  the  net  tilt  of 
the  body.  However,  between  these  limits,  Yeadon  claims  the  timing  of  the  extension 
is  not  critical  to  the  hnal  tilt  angle. 

6.3.1  The  Nominal  Case 

I  created  a  simulation  of  a  human  performing  a  one  and  one  half  twisting  front  somer¬ 
sault.  The  simulated  athlete  used  the  technique  described  in  the  previous  paragraph 
to  perform  the  maneuver  in  a  weightless  environment  (Figure  6-1).  The  human  model 
has  thirteen  joint  degrees  of  freedom,  including  three  in  each  shoulder,  one  at  each 
elbow,  two  in  each  hip,  and  one  in  the  torso  allowing  the  upper  body  to  twist  relative 
to  the  lower  body.  The  maneuver  was  initialized  from  an  upright  piked  somersault 
position  with  the  arms  held  straight  out  from  the  side  of  the  body.  The  maneuver 
hnished  in  a  layout  body  position  with  the  arms  held  straight  out  from  the  sides. 

The  control  movements  used  to  produce  this  maneuver  were  hand  programmed 
so  that  the  model  looked  natural  during  the  movement  and  hnished  with  the  desired 
attitude  and  body  conhguration.  The  control  movements  consisted  of  a  sequence  of 
desired  positions  of  the  joints  written  as  functions  of  time.  The  purpose  of  designing 
this  maneuver  was  to  empirically  hud  a  sequence  of  body  conhgurations  and  joint 
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Figure  6-1:  Images  arranged  in  right-to-left,  top-to-bottom  order  from  a  dynamic 
simulation  of  a  1  1/2  twisting  front  somersault.  The  maneuver  was  initialized  from 
somersault  rotation  in  the  piked  position.  The  control  movements  for  this  maneuver 
are  a  hand  programmed  sequence  of  joint  angles  written  as  functions  of  time. 

torques  that  would  produce  the  desired  maneuver  under  nominal  conditions.  The 
nominal  somersault  rate  at  the  beginning  of  the  maneuver  was  12.6  radj sec  and  the 
maximum  twist  rate  during  the  maneuver  is  40.2  radj  sec. 

6.3.2  OfF-Nominal  Performance,  Prescribed  Motion 

If  the  simulated  performer  uses  the  prescribed  set  of  control  movements  from  a  differ¬ 
ent  set  of  initial  conditions  than  the  nominal  case  then  the  trajectory  of  the  maneuver 
will  change.  How  sensitive  the  maneuver  is  to  variation  in  the  initial  conditions  is  im¬ 
portant  to  the  reliability  of  the  maneuver.  We  tested  the  sensitivity  of  this  maneuver 
to  initial  conditions  by  running  a  series  of  simulations,  each  starting  from  a  different 
set  of  initial  tilt  (0o)  and  twist  (To)  angles  of  the  body.  We  evaluated  the  reliability 
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of  the  maneuvers  by  comparing  the  landing  attitude  of  the  off-nominal  simulations 
to  that  of  the  nominal  case.  (The  landing  attitude  was  considered  to  be  the  attitude 
at  a  hxed  time  after  the  start  of  the  maneuver.)  The  error,  e  in  landing  attitude  was 
computed  as  follows: 

n  joints 

C  =  l/4(($,  -  +  (0,  -  0)2  +  (T,  -  T)2  +  1/13  ^  {qd,  -  qif) 

8  =  1 

where  0,  and  T  are  the  body  attitude  at  landing  qi  refers  to  the  joint  position 
and  the  subscript  d  refers  to  the  desired  value.  This  equation  for  the  error  emphasizes 
the  body  attitude  over  body  conhguration.  An  error  of  1.0  is  large;  it  could  mean  the 
twist  angle  or  somersault  angle  was  off  by  4.0  radians  or  about  270°. 

The  results  of  a  series  of  simulations  that  varied  the  initial  tilt  and  twist  attitude 
of  the  body  over  a  range,  —0.1  rad  <  0o  A  0.1  rad,  —0.1  rad  <  To  <  0.1  rad  are 
shown  in  Figure  6-2. 

The  nominal  maneuver,  0o  =  0,  To  =  0,  corresponds  to  the  center  grid  point 
of  this  hgure.  The  height  of  the  surface  there  is  zero.  Away  from  the  nominal  the 
landing  attitude  error  increases  except  for  a  narrow  valley  of  initial  tilt  and  twist 
attitudes  along  which  the  landing  error  remains  small.  It  appears  that  a  prescribed 
motion  strategy  would  not  produce  reliable  11/2  twisting  somersaults.  Is  there  a 
simple  open  loop  strategy  that  can  improve  this  performance? 

6.3.3  OfF-Nominal  Performance,  Compliant  Motion 

It  seems  unlikely  that  people  could  accurately  reproduce  prescribed  motions  in  a 
dynamic  movement  like  the  11/2  twisting  somersault.  People  use  springy  muscles 
and  tendons  to  position  their  limbs.  It  seems  likely  that  the  change  in  environmental 
forces  that  would  accompany  a  change  in  maneuver  trajectory  would  mean  that  the 
limb  movements  change  even  if  the  athlete  tried  to  execute  the  exact  same  motions. 
This  idea  is  the  basis  for  the  following  set  of  simulation  experiments. 
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Figure  6-2:  This  figure  shows  the  accuracy  of  a  1  1/2  twisting  somersault  in  the 
presence  of  off-nominal  initial  conditions.  The  simulated  performer  used  a  prescribed 
set  of  joint  angles  to  perform  the  maneuver.  The  two  grid  axes  correspond  to  initial 
tilt  and  twist  angles  of  the  body,  —0.1  rad  <  0o  <  0.1  rad,  —0.1  rad  <  To  <  0.1  rad. 
The  height  of  the  grid  indicates  the  accuracy  of  the  landing  attitude  of  the  simulation 
starting  from  the  corresponding  set  of  initial  conditions. 
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We  ran  a  series  of  simulations  of  the  11/2  twisting  front  somersault  that  used  the 
joint  torques  from  the  nominal  maneuver  as  feed  forward  commands.  We  also  used 
the  prescribed  joint  motion  from  the  nominal  maneuver  as  the  commanded  positions 
for  a  set  of  compliant,  position  plus  derivative  (pd)  servos  at  the  joints.  In  this 
open  loop  strategy,  the  net  joint  torques  would  be  a  combination  of  the  feed  forward 
torques  plus  pd  servo  torques.  Since  the  pd  servos  use  a  prescribed  set  of  positions 
as  desired  values,  they  act  like  passive  springs  and  dampers  at  the  joints.  If  the 
maneuver  started  from  the  nominal  initial  conditions  then  the  joints  would  follow  the 
prescribed  trajectories.  Then  the  torques  from  the  pd  servos  would  remain  zero  and 
the  nominal  maneuver  would  be  reproduced. 

Our  goal  in  this  part  of  the  experiment  was  to  hud  a  set  of  pd  servo  gains  (spring 
constants  and  damping  coefficients)  that  produced  the  most  reliable  performances. 
Rather  than  search  over  the  gain  parameters  of  pd  servos  at  thirteen  joints,  I  chose  to 
search  over  one  parameter.  Therefore,  I  compute  the  gains  of  all  joint  servos  according 
to  a  single  parameter,  the  body  clamped  natural  frequency,  cOn-  The  intent  behind 
this  choice  is  that  all  body  joints  have  a  similar  compliance  or  natural  frequency  of 
operation.  We  will  then  look  for  a  body  natural  frequency  that  produces  reliable  1 
1/2  twisting  somersaults. 

To  compute  the  servo  gains  at  a  given  joint  as  a  function  of  I  assume  that  the 
body  inboard  (towards  the  torso)  from  the  joint  is  inertially  hxed  and  all  out-board 
joints  are  immobilized.  This  way  the  model  simplihes  to  a  single  d.o.f.  joint  between 
the  limb  in  question  and  ground.  The  equations  reduce  to  a  simple  second  order 
system  of  equations  as  follows 

I X  c(x  —  Xd)  -\-  k{x  —  Xd)  =  0 

where  /  is  the  apparent  inertia  of  the  limb  at  the  joint,  c  is  the  damping  coefficient, 
and  k  is  the  spring  constant  of  the  joint.  This  simple  system  can  be  rewritten  in  the 
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canonical  form 

X  +  2C,LOn{x  -  Xd)  +  L0l{x  -  Xrf)  =  0 

therefore  c//  =  2C,iOn  and  kj I  =  We  choose  =  0.7  to  achieve  a  well  damped 
system  response  [Ogata],  Then  choosing  LOn  allows  us  to  compute  the  stiffness  and 
damping  constants  in  a  consistent  manner. 

We  systematically  varied  the  value  of  LOn  between  150  rad/s  and  20  rad/s.  At  each 
of  these  values  of  the  body  natural  frequency  we  performed  a  series  of  simulations 
starting  from  initial  conditions  just  as  in  the  case  of  prescribed  motion  at  the  joints. 
I  plotted  the  results  in  the  form  of  the  3D  plot  shown  in  Figure  6-2  for  each  frequency 
and  subjectively  evaluated  them.  There  was  a  clear  choice  for  the  best  natural  fre¬ 
quency.  Figure  6-3  shows  the  simulation  results  for  the  case  of  LOn  =  ‘iQradjs.  This 
plot  shows  that  the  reliability  of  the  11/2  twisting  somersault  performed  with  the 
open  loop  strategy  shows  marked  improvement  over  the  prescribed  motion  case.  Now, 
nearly  half  the  set  of  initial  conditions  results  in  small  hnal  attitude  errors.  Further¬ 
more,  the  landing  attitude  errors  were  worse  for  either  smaller  or  larger  values  of 
iOn-  These  basic  results  held  true  for  variations  in  the  angular  momentum,  h,  of  the 
maneuver  as  well  (0.9  h  <  h  <  1.1  h). 

Figures  6-4  and  6-5  show  data  from  two  simulations  using  the  nominal  value  of 
angular  momentum  and  with  0o  =  —0.1,  and  To  =  —0.1.  The  data  of  hgure  6-4 
shows  that  when  prescribed  joint  motion  is  used  for  an  off-nominal  maneuver,  the 
body  attitude  error  becomes  large  at  the  end  of  the  maneuver.  In  contrast,  the  data 
of  hgure  6-5  shows  that  while  the  joint  angles  incur  some  error  during  the  maneuver 
the  body  attitude  is  close  to  the  desired  value  at  the  end  of  the  maneuver. 

Since  the  timing  of  this  maneuver  will  scale  with  somersault  rate  it  is  instructive 
to  show  how  the  body  clamped  natural  frequency  compares  to  the  nominal  somersault 
rate,  ^  =  2.38.  It  appeared  in  simulations  that  the  compliance  that  allowed 

the  arms  to  open  prematurely,  thus  increasing  the  body  inertia  about  the  twist  axis, 
was  important  in  off-nominal  simulations.  For  this  reason,  one  might  consider  scaling 


CHAPTER  6.  TWISTING  SOMERSAULTS 


132 


Wn 


30 


Figure  6-3:  This  figure  shows  the  accuracy  of  a  1  1/2  twisting  somersault  in  the  pres¬ 
ence  of  off-nominal  conditions.  The  simulated  performer  used  feed  forward  torques 
plus  passive  pd  servos  at  the  joints  to  perform  the  maneuver.  Pd  servo  gains  were 
chosen  according  to  Un  =  ‘iWadj s.  The  pd  servos  used  the  prescribed  joint  angles 
from  the  nominal  maneuver  as  desired  values.  The  two  grid  axes  correspond  to  initial 
tilt  and  twist  angles  of  the  body,  —0.1  rad  <  0o  <  0.1  rad,  —0.1  rad  <  To  <  0.1  rad. 
The  height  of  the  grid  indicates  the  accuracy  of  the  landing  attitude  of  the  simulation 
starting  from  the  corresponding  set  of  initial  conditions. 
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Figure  6-4:  Simulation  data  from  all/2  twisting  somersault  that  used  prescribed 
control  motions  for  initial  conditions  0o  =  —0.1,  and  =  —0.1  are  shown  with  solid 
lines.  The  body  attitude  and  joint  angles  of  the  nominal  maneuver  are  shown  with 
dashed  lines.  While  the  prescribed  joint  positions  are  accurate  for  the  off-nominal 
case  (bottom  three  graphs),  the  body  landing  attitude  error  is  large.  (The  traces  stop 
at  the  landing  time.) 
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Figure  6-5:  Simulation  data  from  all/2  twisting  somersault  with  Un  =  30,  0o  = 
—0.1,  and  4^0  =  —0.1  are  shown  with  solid  lines.  The  desired  values  of  the  body 
attitude  and  joint  angles  are  shown  with  dashed  lines.  Joint  angles  incur  signihcant 
tracking  errors  but  the  landing  attitude  closely  follows  that  of  the  nominal  maneuver. 
(The  traces  stop  at  the  landing  time.) 
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the  body  clamped  natural  frequency  by  the  maximum  twist  rate  of  the  nominal 
maneuver,  =  wrx  =  0.746. 

The  body  clamped  natural  frequency  was  chosen  to  be  a  simple  parameter  that 
described  the  natural  frequency  of  the  whole  body.  However,  since  the  computation 
of  this  parameter  assumed  that  parts  of  the  body  were  inertially  hxed,  the  actual 
eigenfrequencies  of  the  system  will  differ  from  this  body  clamped  frequency.  The 
actual  eigenfrequencies  of  the  body  for  the  initial  conhguration  (pike)  and  the  twisting 
conhguration  (wrap)  of  the  Rudi,  are  shown  in  Table  6.3.3. 

6.4  Twisting  Somersault  of  the  3D  Biped 

We  programmed  a  simulated  3D  Biped  robot  to  use  the  “torque  twist”  method  to 
initiate  a  twisting  somersault.  To  use  this  technique  the  simulated  3D  Biped  robot 
accelerated  its  torso  in  somersault  and  twist  during  the  stance  phase  just  prior  to 
take-off.  These  two  components  of  momentum  should  result  in  an  angular  momentum 
vector  that  is  tilted  with  respect  to  the  principal  axes  of  the  robot. 

Figure  3-3  shows  the  map  of  rotational  motion  of  the  3D  Biped.  When  compared 
to  the  map  for  the  layout  somersault  in  the  human  (Figure  2-7)  we  see  that  the  robot 
requires  signihcantly  more  tilt  of  the  principal  axes  (relative  to  the  angular  momentum 
vector)  to  achieve  a  twisting  somersault.  This  is  partly  due  to  the  fact  that  the  robot 
somersault  axis  is  the  major  principal  axis  but  also  due  to  the  signihcantly  different 
shape  of  the  regions  on  the  two  maps.  The  qualitative  difference  in  shape  of  the  maps 
is  due  to  the  difference  in  rigid  body  inertia  ratios  of  the  robot  and  of  the  human. 
To  put  it  simply  the  human  is  long  and  skinny  and  the  robot  is  short  and  fat.  This 
difference  makes  the  twist  harder  to  achieve  (more  tilt  required),  and  it  means  the 
maneuver  will  not  look  much  like  that  of  a  human. 

In  its  current  conhguration,  the  twisting  somersault  mode  of  the  3D  Biped  would 
force  the  robot  to  lay  on  its  side  some  time  during  hight.  This  body  orientation  is  not 
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Table  6.1 

Rudi  Eigenfrequencies 


Pike 

Wrap 

iOn  radj sec 

c 

LOn  radj  sec 

C 

921.0 

0 

210.0 

0 

709.2 

0 

188.3 

0 

359.4 

0 

159.2 

0 

161.2 

0 

107.7 

0 

160.6 

0 

93.1 

0 

88.0 

0 

39.7 

0.92 

69.0 

0 

39.7 

0.92 

36.7 

0.85 

38.9 

0.90 

36.7 

0.85 

38.9 

0.90 

31.3 

0.73 

30.1 

0.70 

31.3 

0.73 

30.1 

0.70 

31.0 

0 

30.0 

0.70 

30.1 

0.70 

30.0 

0.70 

30.1 

0.70 

27.8 

0 

30.0 

0.70 

26.8 

0.62 

30.0 

0.70 

26.8 

0.62 

28.3 

0 

26.7 

0 

24.7 

0 

26.5 

0.62 

24.7 

0 

26.5 

0.62 

22.9 

0.53 

24.9 

0.58 

22.9 

0.53 

24.9 

0.58 

22.7 

0 

24.8 

0.58 

22.0 

0 

24.8 

0.58 

21.9 

0 

24.7 

0 

16.0 

0.37 

24.1 

0 

16.0 

0.37 

23.8 

0 
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what  we  identify  with  a  twisting  somersault  in  a  human.  The  human’s  long  axis  stays 
closer  to  the  vertical  during  a  twisting  somersault.  In  addition,  this  inertial  orientation 
of  the  robot  is  undesirable  from  a  practical  standpoint  because  this  orientation  is 
coincident  with  gimbal  lock  in  the  gyros  used  to  measure  body  attitude.  During 
laboratory  experiments,  the  physical  robot  achieved  this  horizontal  position  which 
resulted  in  damage  to  the  gyroscopes.  To  correct  this  situation  we  changed  the  robot 
inertia,  and  thus  its  rotational  modes,  to  look  more  like  that  of  a  human. 

To  change  the  robot  inertia  we  added  weight  to  increase  the  major  and  interme¬ 
diate  principal  inertias  without  increasing  the  minor  principal  inertia.  We  did  this 
by  adding  weight  along  the  ’3’  axis  of  the  robot.  Figure  6-6  shows  the  new  map  of 
rotation  for  the  3D  Biped  with  a  I.Okg  weight  added  0.8  m  above  the  hips  of  the 
robot.  This  map  looks  much  more  like  that  of  a  human  now  and  the  twisting  somer¬ 
sault  should  bear  this  resemblance  as  well.  Figure  6-7  shows  a  sequence  of  computer 
graphic  images  of  a  simulated  3D  Biped  with  added  weight  (and  strong  actuators) 
performing  a  front  somersault  with  one  half  twist.  Notice  that  the  simulated  robot 
lands  the  twisting  somersault  facing  the  opposite  direction  it  started  from.  Data 
from  the  simulation  is  included  in  Figure  6-8.  The  simulated  running  robot  regained 
balance  on  landing  to  continue  stable  dynamic  running. 

We  tried  this  maneuver  with  the  physical  3D  Biped  robot  in  the  laboratory.  The 
added  weight  and  inertia  were  too  large  for  the  robot  to  achieve  sufficient  flight  time 
or  angular  momentum  to  produce  the  maneuver. 


6.5  Summary 

In  this  chapter  I  presented  simulation  results  of  a  1  1/2  twisting  front  somersault 
performed  by  a  model  gymnast  with  thirteen  joint  degrees  of  freedom.  The  maneuver 
is  initiated  from  a  piked  front  somersault.  The  twisting  maneuver  results  from  a 
sequence  of  movements  of  the  limbs  and  torso  made  during  flight.  We  found  that 
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Figure  6-6:  Spherical  coordinate  map  of  the  tilt  and  twist  Euler  angle  trajectories  for 
the  3D  Biped  with  I  kg  of  weight  added  at  a  distance  of  0.8  m  above  the  hip. 
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Figure  6-7:  Sequence  of  images  of  the  simulated  3D  Biped  somersault  with  twist. 
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Figure  6-8:  Hopping  height,  somersault,  tilt  and  twist  Euler  angles  of  the  3D  Biped 
during  a  simulated  front  somersault  with  1/2  twist.  The  simulated  robot  passes 
through  an  Euler  angle  singularity  at  approximately  1.2  sec.  causing  discontinuties 
in  the  data.  This  singularity  did  not  affect  the  dynamic  simulation  however  as  a 
different  set  of  attitude  parameters  were  used. 
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prescribed  limb  movements  could  produce  reliable  11/2  twisting  somersaults  only  for 
a  small  set  of  off-nominal  initial  conditions.  In  contrast,  a  control  strategy  that  used 
feed  forward  joint  torques  in  a  tuned  passive  dynamic  model  of  the  performer  could 
produce  reliable  maneuvers  for  a  much  larger  set  of  off-nominal  initial  conditions. 

We  do  not  yet  know  the  signihcance  of  this  result.  There  are  many  degrees  of 
freedom  of  this  model.  Too  many  to  easily  distinguish  between  important  and  unim¬ 
portant  effects.  Simplihcation  of  this  model  may  reveal  salient  features  that  can  be 
analytically  conhrmed.  These  simulation  results  from  the  rudi  are  potentially  inter¬ 
esting  because  it  suggests  that  a  passive  dynamic  approach  to  twisting  somersaults 
may  be  capable  of  producing  reliable  maneuvers  as  was  the  case  with  the  layout  som¬ 
ersault.  It  is  interesting  that  many  of  the  eigenfrequencies  of  the  most  reliable,  tuned 
compliant  system  for  the  rudi  were  in  the  neighborhood  of  the  maximum  twist  rate  of 
the  maneuver.  This  appears  to  be  analogous  to  the  ideal  choice  of  arm-body  oscilla¬ 
tion  in  the  passive  layout  somersault.  In  this  case  the  best  shoulder  spring  resulted  in 
an  eigenfrequency  of  the  arm-body  oscillation  that  was  equal  to  the  somersault  rate. 

In  this  chapter,  I  also  described  simulation  experiments  with  3D  Biped  twisting 
somersaults.  The  distribution  of  mass  in  the  3D  Biped  robot  makes  a  twisting  som¬ 
ersault  particularly  difficult.  In  order  to  enter  a  rotational  mode  that  could  produce 
the  twist  angles  desired,  the  robot  had  to  assume  body  attitudes  unlike  those  asso¬ 
ciated  with  a  human  somersault  with  twist.  This  extreme  body  attitude  could  lead 
to  mechanical  difficulties  with  inertial  instruments  on  the  physical  robot.  We  added 
weight  to  the  robot  so  its  natural  rotational  modes  were  more  like  those  of  a  human. 
A  simulated  3D  Biped  robot  was  able  to  perform  a  1/2  twisting  front  somersault.  The 
simulated  robot  landed  the  maneuver  and  continued  running  stably  afterwards.  The 
added  weight  made  twisting  somersaults  of  the  physical  3D  Biped  robot  impossible 
due  at  least  in  part  to  insufficient  actuator  power. 


Chapter  7 


Summary  and  Discussion 


Inspiration  for  this  thesis  comes  from  the  remarkable  performance  of  aerial  maneu¬ 
vers  by  gymnasts,  divers,  trampolinists,  and  skiers.  From  a  distance,  the  control  of 
aerial  maneuvers  appears  to  be  a  very  complex  task.  The  shortage  of  control  during 
flight  and  the  non-linear  rotational  mechanics  make  this  problem  more  than  a  little 
daunting.  However,  adversity  breeds  creativity  and  athletes  and  other  people  have 
been  very  creative  in  working  with  these  difficulties  to  produce  elegant  solutions  to 
movement  control  problems.  It  appears  that  we  can  learn  something  from  human  ca¬ 
pability  in  making  dynamic  flying  or  floating  machines  that  can  reorient  themselves 
quickly,  reliably,  and  accurately. 

This  thesis  is  concerned  with  how  people  can  incorporate  known  aerial  movement 
techniques  into  a  strategy  that  produces  reliable  performances.  People  can  influence 
aerial  maneuvers  through  the  relative  movement  of  their  limbs  and  torso  during  flight. 
While  a  hxed  sequence  of  control  movements  can  be  found  to  produce  a  desired 
maneuver,  prescribed  limb  motion  may  provide  little  reliability  when  a  maneuver  is 
subject  to  the  variety  of  conditions  that  would  accompany  any  real  performance.  It 
is  possible  that  people  produce  reliable  maneuvers  by  sensing  these  variations  and 
actively  computing  responses  to  compensate  for  them.  However,  the  complexity  of 
this  approach  warrants  the  search  for  a  simpler  strategy.  A  persistent  goal  of  this  work 
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has  been  to  find  control  strategies  that  do  not  depend  entirely  on  active  feedback, 
that  are  simple,  and  that  can  work  in  laboratory  machines. 

The  contribution  of  this  thesis  is  to  point  out  via  example  that  passive  dynamic 
solutions  to  movement  control  are  an  alternative  to  certain  problems  previously  con¬ 
sidered  to  require  active  feedback  control.  We  focus  on  an  open  loop,  passive  dynamic 
strategy  for  producing  gymnastic  maneuvers.  Open  loop  means  that  control  effort 
(torques  at  joints)  are  simply  replayed  from  memory  during  a  maneuver  (the  motor 
tape  model  of  biological  control).  Under  identical  conditions  open  loop  control  will 
produce  identical  maneuvers.  Using  open  loop  control  under  different  initial  condi¬ 
tions  will  produce  different  outcomes  due  to  the  inherent  dynamics  of  the  system. 
This  may  make  open  loop  control  by  itself  incapable  of  producing  reliable  maneuvers. 
However,  it  is  possible  that  the  passive  dynamic  behavior  of  the  performer’s  body 
could  automatically  compensate  for  off-nominal  conditions  in  a  way  that  produces 
reliable  maneuvers.  The  passive  dynamic  behavior  of  a  system  can  be  tuned  through 
the  selection  of  passive  elements  like  springs  (springy  muscles)  or  by  choosing  nominal 
body  configurations  during  the  maneuver.  Careful  tuning  of  the  passive  dynamics  of 
a  system  coupled  with  open  loop  control  may  provide  a  strategy  for  producing  reliable 
gymnastic  maneuvers. 

This  strategy  is  simple  in  the  sense  that  it  requires  no  computation  during  the 
maneuver.  Only  memory  playback  is  required.  Compensatory  movements  are  ’com¬ 
puted’  by  the  physical  system  as  part  of  its  natural  behavior.  Incorporation  of  this 
strategy  could  simplify  gymnastic  performances  by  reducing  the  amount  of  active 
control  and  computation  required  by  the  athlete. 

To  the  extent  possible,  we  have  insisted  upon  laboratory  experimentation  to  test 
our  ideas.  We  feel  that  laboratory  experimentation  on  physical  machines  forces  one 
to  discover  the  salient  features  of  a  problem  that  may  be  difficult  to  reveal  otherwise. 
(Another  reason  is  that  laboratory  experimentation  is  also  a  lot  of  fun!)  However, 
while  a  laboratory  experiment  can  be  used  to  validate  a  strategy  for  machine  control. 
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it  can  only  suggest  the  viability  of  a  strategy  for  human  use. 

We  studied  three  different  gymnastic  maneuvers  in  this  thesis,  the  tucked  somer¬ 
sault,  the  layout  somersault,  and  the  11/2  twisting  front  somersault. 

7.1  Robot  Tucked  Somersaults 

The  tucked  somersault  is  a  common  maneuver  in  athletic  events.  The  main  require¬ 
ment  of  this  maneuver  is  a  balanced  landing  which  in  turn  requires  a  precise  body 
attitude.  We  programmed  a  3D  Biped  robot  to  perform  front  somersaults  in  the 
laboratory.  The  control  strategy  relied  implicitly  upon  the  passive  dynamic  stability 
of  the  tucked  somersault.  We  used  an  active  feedback  system  to  control  somersault 
rotation  rate  and  foot  placement  at  landing. 

A  linear  analysis  of  rigid  body  rotation  showed  that  the  tilt  and  twist  attitude  of 
the  somersaulting  3D  Biped  robot  is  passively  stable.  However  somersault  attitude 
is  unstable.  This  well  known  result  led  us  to  use  active  control  of  the  somersault 
degree  of  freedom  in  the  3D  Biped  maneuver.  The  robot  controls  somersault  rotation 
rate  by  tucking  or  untucking  its  legs  during  flight  to  change  the  rotational  inertia. 
During  flight  we  depend  upon  a  stable  passive  dynamic  response  of  the  tilt  and  twist 
dynamics  to  insure  small  tilt  and  twist  angles  at  landing.  To  keep  the  initial  tilt 
and  twist  angles  small  we  use  a  wide  double  stance  of  the  robot  during  take-off  and 
insure  that  the  feet  touch  down  simultaneously.  Another  important  element  of  active 
control  of  the  3D  Biped  somersault  was  the  placement  of  the  feet  prior  to  landing. 
This  element  was  necessary  to  compensate  for  errors  in  the  estimated  landing  time. 
On  its  best  day  the  robot  has  performed  successful  front  somersaults  on  seven  out  of 
ten  attempts  in  laboratory  experiments. 

We  do  not  know  if  active  control  of  rotation  rate  in  the  3D  Biped  robot  is  necessary 
to  perform  the  front  somersault.  It  is  possible  that  open  loop  control  plus  active  foot 
placement  prior  to  landing  could  also  be  used  to  perform  reliable  front  somersaults. 
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This  was  the  case  with  the  planar  biped  somersault  work  of  Hodgins  and  Raibert.  A 
strategy  that  we  would  like  to  experiment  with  in  the  future  could  provide  a  middle 
ground  between  the  active  rotation  rate  control  used  in  this  study  and  the  (primarily) 
open  loop  approach  used  in  the  planar  biped.  This  strategy  would  use  tuned  leg 
springs  to  passively  influence  somersault  rotation  rate.  This  technique  could  correct 
rotation  rate  errors  by  allowing  the  legs  to  extend  at  higher  than  normal  rotation 
rates  (due  to  higher  centrifugal  forces)  and  vice  versa  for  slower  rotation  rates.  This 
passive  dynamic  approach  coupled  with  active  foot  placement  may  provide  reliable 
front  somersaults  with  less  active  feedback  control. 


7.2  The  Layout  Somersault 

The  layout  somersault  involves  rotation  about  the  middle  principal  axis  of  inertia, 
an  unstable  rotation  for  a  rigid  body.  Biomechanics  researchers  have  suggested  that 
human  athletes  use  active  feedback  control  during  flight  to  stabilize  the  layout  somer¬ 
sault.  We  found  that  the  layout  somersault  could  passively  be  made  neutrally  stable. 

The  layout  somersault  is  stabilized  by  passive  arm  movement  in  the  frontal  plane. 
The  inherent  tendency  of  the  arms  to  tilt  in  response  to  twisting  movement  of  the  body 
provides  a  built-in  correction  for  the  layout  instability.  The  arm  accomplishes  this  by 
changing  the  orientation  of  the  principal  axes  with  respect  to  the  angular  momentum 
vector.  This  built-in  correction  eliminates  the  divergent  tendency  of  the  system  as 
long  as  a  carefully  selected  shoulder  spring  is  used  to  cancel  unstable  centrifugal  forces 
on  the  arms.  These  results  are  conhrmed  with  linear  stability  analysis,  non-linear 
dynamic  simulation  and  laboratory  experiments  with  a  somersaulting  doll.  During 
experiments  the  doll  can  consistently  perform  at  least  three  and  one  half  stable  layout 
somersaults. 

Analysis  of  a  simplihed  model  of  the  layout  somersault  dynamics  revealed  several 
salient  features  of  the  movement.  The  simplihed  model  revealed  the  simple  rule 
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for  picking  stabilizing  non-dimensional  shoulder  springs,  K^h  =  1-0.  This  spring 
constant  insures  that  the  natural  frequency  of  the  arm-body  oscillation  is  matched 
to  the  rotation  rate  of  the  body.  It  also  is  the  exact  spring  required  to  cancel  the 
destabilizing  centrifugal  forces  on  the  arms.  The  empirically  determined,  best  value 
of  Ksh  for  the  experimental  doll  was  K^h  =  0.82.  This  value  is  in  agreement  with  the 
theoretical  results.  I  think  that  it  is  signihcant  that  the  stabilizing  shoulder  spring 
just  cancels  the  spring-like  action  of  centrifugal  forces.  This  balance  between  opposing 
forces  allows  the  arm  to  be  responsive  to  the  gyric  forces  which  ultimately  stabilize 
the  sytem. 

The  simple  model  also  helped  explain  why  arm  orientation  is  important  in  stabi¬ 
lizing  the  layout  somersault.  The  tendency  of  the  arms  to  tilt  in  response  to  twist  is 
greatest  when  the  arms  are  held  straight  out  to  the  side.  This  tendency  is  decreased 
as  the  arms  are  raised  or  lowered  from  that  position.  Another  feature  that  the  simple 
model  revealed  is  that  stabilization  of  the  layout  somersault  is  easier  with  bigger  and 
bigger  arms.  This  is  not  surprising  for  in  the  limit  when  the  arms  (held  straight  out 
to  the  side)  are  much  larger  than  the  body,  the  system  resembles  a  rigid  body  rotating 
about  its  minimum  principal  axis,  a  well  known  gyrically  stabilized  conhguration. 

We  ignored  the  control  of  somersault  angle  in  the  layout  somersault.  Stabiliza¬ 
tion  of  the  somersault  angle  will  be  important  for  practical  layout  somersaults  as  it 
was  for  the  3D  Biped  tucked  somersaults.  However,  decoupling  of  the  tilt,  twist, 
and  asymmetric  arm  movement  from  the  somersault  and  symmetric  arm  movement 
suggests  that  a  somersault  control  strategy  and  a  tilt-twist  control  strategy  could  be 
developed  separately  then  superposed  in  the  complete  system. 

7.3  Twisting  Somersaults 

Twisting  somersaults  are  among  the  most  interesting  and  complex  of  aerial  maneu¬ 
vers.  They  involve  simultaneous  rotation  about  the  twist  axis  and  the  somersault 
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axis.  Their  dynamic  complexity  makes  feedback  control  of  these  maneuvers  chal¬ 
lenging.  Is  it  possible  that  open  loop  control  can  be  used  to  reliably  produce  these 
maneuvers?  Perhaps  it  can  if  the  performers  body  acts  like  a  tuned  passive  dynamic 
system.  Our  results  on  the  open  loop,  passive  dynamic  control  of  twisting  somersaults 
are  of  a  preliminary  nature.  However,  they  point  to  a  potentially  interesting  subject 
for  future  work.  These  results  suggest  that  passive  dynamics  may  help  the  reliable 
performance  of  twisting  maneuvers  as  was  the  case  for  the  layout  somersault. 

We  used  a  thirteen  joint  human  model  to  simulate  the  performance  of  a  1  1/2 
twisting  front  somersault.  The  twisting  maneuver  was  initiated  from  a  piked  front 
somersault  using  asymmetric  movement  of  the  limbs.  We  found  that  using  prescribed 
limb  motion  during  the  maneuver  produced  inaccurate  landing  attitudes  when  the 
initial  body  attitude  was  allowed  to  vary.  On  the  other  hand,  an  open  loop,  passive 
dynamic  strategy  was  able  to  produce  reliable  landing  attitude  for  a  relatively  large 
set  of  initial  conditions.  In  this  strategy,  open  loop  torque  commands  from  the  nom¬ 
inal  trajectory  were  combined  with  torques  from  passive  springs  and  dampers  at  the 
joints.  We  tuned  the  passive  compliance  of  the  performer’s  joints  to  hud  the  value  of 
compliance  that  produced  the  most  reliable  maneuvers.  We  used  a  single  parameter 
to  characterize  the  compliance  of  all  the  body  joints.  There  was  a  clear  choice  for 
the  best  system  compliance.  More  or  less  compliance  at  the  joints  led  to  less  reliable 
performance.  Inspection  of  the  eigenfrequencies  of  the  model  revealed  that  the  most 
reliable  system  had  natural  frequencies  in  the  vicinity  of  the  maximum  twist  rate  of 
the  maneuver.  As  in  the  layout  somersault,  it  may  be  that  a  balance  between  the 
centrifugal  forces  of  twisting  and  the  spring  forces  of  the  joints  is  required  to  make 
the  maneuver  reliable. 

While  the  reliability  of  the  open  loop,  passive  dynamic  system  was  better  than  that 
of  the  prescribed  motion  system  there  is  still  room  for  improvement.  There  were  sets  of 
initial  conditions  that  were  ’close’  to  nominal  for  which  landing  attitude  performance 
was  seriously  degraded.  Could  a  different  passive  system  improve  reliability  in  this 
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region  as  well?  Perhaps  a  simple  active  controller  could  complement  the  passive 
dynamics  to  produce  an  even  more  reliable  system.  Yeadon  has  discussed  the  need 
for  good  timing  in  certain  maneuvers.  Perhaps  a  simple  active  controller  that  changed 
the  replay  rate  of  a  set  of  pre-recorded  actions  based  on  perceived  errors  could  get 
the  timing  ‘right’  to  further  improve  maneuver  reliability. 

7.4  Do  Humans  Use  Passive  Dynamics 

We  do  not  know  if  humans  use  a  passive  dynamic  approach  to  stabilization  of  gym¬ 
nastic  maneuvers.  This  thesis  can  not  prove  or  disprove  the  human  use  of  such  a 
strategy.  The  results  of  this  thesis  can  only  suggest  that  such  a  strategy  is  a  viable 
one  for  particular  maneuvers,  or  parts  of  maneuvers.  What  experiments  could  we  do 
to  learn  more  about  how  people  perform  these  maneuvers? 

It  would  be  useful  to  measure  the  range  of  initial  conditions  from  which  athletes 
can  reliably  perform  aerial  maneuvers,  and  also  measure  their  corresponding  limb 
movements.  If  the  variations  in  initial  conditions  and  control  movements  were  large 
then  one  may  argue  that  dynamic  compensation  rather  than  a  prescribed  motion 
strategy  was  at  work.  One  could  look  for  correlations  between  different  initial  condi¬ 
tions  and  different  control  movements.  Could  a  passive  dynamic  model  explain  any 
observed  correlations? 

Tests  of  the  compliance  of  human  muscle  could  reveal  if  the  effective  spring  con¬ 
stants  of  human  limbs  are  near  the  values  predicted  by  theory  to  produce  passively 
stable  maneuvers.  This  could  support  a  passive  dynamic  theory  but  not  prove  it’s  use. 
Measurements  of  the  electromyographic  signal  of  the  muscles  may  not  help  to  decide 
what  strategy  people  use.  Even  with  accurate  recordings  of  many  performances  of 
a  maneuver  it  would  be  difficult  to  distinguish  between  a  motor  tape  model  and  an 
active  control  model  of  muscle  activation. 

A  difficulty  in  making  a  distinction  between  an  active  control  strategy  and  an 
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open  loop,  passive  dynamic  strategy  is  that  both  approaches  could  produce  a  simi¬ 
lar  dynamic  response.  In  fact,  one  may  argue  that  if  there  is  one  clear  technique  to 
achieving  a  reliable  maneuver  then  both  strategies  would  have  to  use  it.  This  possibil¬ 
ity  is  illustrated  by  the  fact  that  the  active  controller  designed  by  Yeadon  to  stabilize 
the  layout  somersault  used  the  same  control  effect  as  the  passive  strategy  presented 
in  this  thesis.  Both  techniques  used  arm  tilt  to  control  body  twist.  Yeadon’s  tech¬ 
nique  depended  upon  sensing  the  twist  rate  of  the  body  and  using  it  to  compute  an 
appropriate  arm  tilt  response.  Whereas,  in  the  passively  stable  case  the  arm  tilt  was 
produced  from  the  inherent  dynamics  of  the  system.  Passive  dynamic  control  is  a 
subset  of  the  space  of  active  controllers  that  require  little  if  any  on-line  computation. 

7.5  A  Passive  Dynamic  Theory  of  Control  De¬ 
sign? 

Can  passive  dynamic  control  be  formulated  into  a  machine  design  and  control  the¬ 
ory?  I  did  not  use  nor  did  I  develop  a  consistent  theory  of  passive  dynamic  design 
for  controlling  gymnastic  maneuvers.  However,  beyond  the  educated  guesses  I  used 
to  hud  some  answers,  elements  of  linear  and  optimal  control  theory  are  central  to  the 
approach  used  to  analyze  the  movements  of  this  thesis.  Linear  analysis  allows  one  to 
examine  the  behavior  of  a  dynamic  system  in  the  neighborhood  of  a  known  solution. 
Biology  provides  us  with  many  examples  of  movements  that  we  know  work.  I  have 
made  explicit  use  of  these  known  solutions  to  search  for  passive  dynamic  stability. 
Tucked  and  layout  somersaults  can  be  analyzed  using  linear  methods.  Linear  analysis 
of  the  twisting  somersault  will  be  more  difficult  but  may  help  reveal  the  important 
features  of  this  maneuver.  Searching  for  passive  dynamic  solutions  to  known  move¬ 
ments  could,  I  think,  be  developed  into  a  theory.  In  fact,  one  may  argue  that  a  passive 
dynamic  theory  of  control  may  consist  of  a  restriction  of  the  optimal  control  theory  to 
passive  dynamic  elements.  Extension  of  a  passive  dynamic  theory  to  unknown  move- 
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ment  solutions  will  be  more  difficult,  reflecting  the  complexity  of  non-linear  systems. 
However,  perhaps  this  is  where  the  greatest  pay-off  from  a  passive  dynamic  theory 
could  lie. 

Why  is  it  so  easy  to  tell  animal  movement  from  machine  movement?  The  richness 
of  human  and  animal  behavior  reflects  the  complexity  of  the  dynamic  systems,  their 
bodies,  they  are  compelled  to  use.  But  this  richness  does  not  just  come  from  complex¬ 
ity,  it  comes  also  from  coordination.  Relatively  speaking,  machine  movement  is  still 
in  its  infancy.  It  would  be  a  grand  goal  of  any  designer  to  make  a  machine  that  could 
move  like  an  animal.  It  is  an  unfortunate  consequence  of  the  complexity  of  non-linear 
systems  that  non-linear  control  design  techniques  frequently  rely  upon  cancellation  of 
the  dynamics  in  order  to  reach  a  solution.  This  approach  risks  the  design  of  forced, 
uncoordinated  behavior.  Perhaps  we  should  search  for  design  techniques  that  em¬ 
brace  the  inherent  dynamics  of  a  system  rather  than  avoid  them.  Perhaps  passive 
dynamic  stability  could  act  as  a  guiding  principal  in  the  design  of  graceful  movement 
in  machines.  Perhaps  it  already  serves  this  purpose  in  animal  movement. 
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A.l  Mathematica  Code  for  Non-linear  Equations 
of  Motion 

Following  is  the  Mathematica  code  for  deriving  the  non-linear  equations  of  motion  of 
the  hve  d.o.f.  system. 

(*  File  for  defining  3D  equations  of  motion  of  a  three  link  body  with 
links  connected  separately  to  the  central  body  by  two  pin  joints 
allowing  rotation  about  the  x  body  axis.  *) 

(*  same  as  manbdof .m  except  arm  angles  defined  in  terms  of  symmetric 
and  antisymmetric  components  phil  =  phia  +  phis,  phir  =  phia  -  phis) 

(*  inertia  matrices  *) 

ib  =  {{ibll , 0 , 0} , {0 , ib22 , 0} , {0 , 0 , ib33}} ;  (*  main  body  *) 
il  =  {{illl,0,0},{0,il22,0},{0,0,il33}};  (*  left  limb  *) 
ir  =  {{irll , 0 , 0} , {0 , ir22 , 0} , {0 , 0 , ir33}} ;  (*  right  limb  *) 

(*  body  to  joint  vector  for  left  arm,  in  arm  coord.  *) 
bjl  =  {bjll,  bjl2,  bjl3}; 

(*  body  to  joint  vector  for  right  arm,  in  arm  coord.  *) 
bjr  =  {bjrl,  bjr2,  bjr3}; 

(*  inboard  body  to  joint  vector  for  left  arm,  in  body  coord.  *) 
ibjl  =  {ibjll,  ibjl2,  ibjl3}; 

(*  inboard  body  to  joint  vector  for  right  arm,  in  body  coord.  *) 
ibjr  =  {ibjrl,  ibjr2,  ibjr3}; 
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(*  transformation  from  body  coordinates  to  limb  left  coordinates  thru  phi 
a  pos.  phi  means  limb  left  rotates  positively  relative  to  body  *) 
cbtol  =  {{1,  0,  0}, 

{0,  Cos[phia[t]  +  phis[t]],  Sin[phia[t]  +  phis[t]]}, 

{0,  -Sin[phia[t]  +  phis[t]],  Cos[phia[t]  +  phis[t]]}}  ; 

cltob  =  Transpose [cbtol]  ; 

(*  transformation  from  body  coordinates  to  limb  right  coordinates  thru  ph 
a  pos.  phir  means  limb  right  rotates  positively  relative  to  body  *) 
cbtor  =  {{1,  0,  0}, 

{0,  Cos[phia[t]  -  phis[t]],  Sin[phia[t]  -  phis[t]]}, 

{0,  -Sin[phia[t]  -  phis[t]],  Cos[phia[t]  -  phis[t]]}}  ; 

crtob  =  Transpose [cbtor]  ; 

(*  inertia  matrices  of  arms  expressed  in  body  coord  *) 
ilb  =  cltob  .  il  .  cbtol; 
irb  =  crtob  .  ir  .  cbtor; 

(*  vector  from  net  eg  to  body  eg,  in  body  coord  *) 

rhob  =  -(ml  (ibjl  -  cltob. bjl)  +  mr  (ibjr  -  crtob . bj r) ) /mt ; 

(*  vector  from  net  eg  to  left  limb  eg,  in  body  coord  *) 
rhol  =  rhob  +  ibjl  -  cltob. bjl; 

(*  vector  from  net  eg  to  right  limb  eg,  in  body  coord  *) 
rhor  =  rhob  +  ibjr  -  crtob. bjr; 
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(*  angular  velocity  of  body  in  body  coord.*) 
wbi  =  {wbl[t],  wb2[t],  wb3[t]}; 

(*  angular  velocity  of  left  limb  w.r.t.  body  in  body  coord.*) 
wlb  =  {phia' [t]  +  phis' [t],  0,  0}; 

(*  angular  velocity  of  limb  in  body  coord.*) 
wli  =  wbi  +  wlb; 

(*  angular  velocity  of  right  limb  w.r.t.  body  in  body  coord.*) 
wrb  =  {phia' [t]  -  phis' [t]  ,  0,  0}; 

(*  angular  velocity  of  right  limb  in  body  coord.*) 
wri  =  wbi  +  wrb; 

(*  wbi  tilde  matrix  for  cross  products  *) 

wbitilde  =  {{0 , -wbi [  [3] ] , wbi [ [2] ] } , {wbi [ [3] ] , 0 , -wbi [ [1] ] } , {-wbi [ [2] ] ,wbi[[l]] ,0}}; 

(*  wli  tilde  matrix  for  cross  products  *) 

wlitilde  =  {{0 , -wli [ [3] ] , wli [ [2] ] } , {wli [ [3] ] , 0 , -wli [ [1] ] } , {-wli [ [2] ] ,wli[[l]] ,0}}; 

(*  wri  tilde  matrix  for  cross  products  *) 

writilde  =  {{0 , -wri [ [3] ] , wri [ [2] ] } , {wri [ [3] ] , 0 , -wri [ [1] ] } , {-wri [ [2] ] ,wri[[l]] ,0}}; 

(*  velocity  of  body  c.g.  *) 

vbi  =  Expand [D [rhob ,t]  +  wbitilde  .  rhob] ; 
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(*  velocity  of  left  limb  c.g.  *) 

(*  all  components  of  rhol  in  body  coord,  so  cross  with  wbi  *) 
vli  =  Expand [D  [rhol ,t]  +  wbitilde  .  rhol]; 

(*  velocity  of  right  limb  c.g.  *) 

(*  all  components  of  rhor  in  body  coord,  so  cross  with  wbi  *) 
vri  =  Expand [D  [rhor ,t]  +  wbitilde  .  rhor]; 

(*  Matrix  of  Partial  velocities  : 

columns  correspond  to  generalized  speeds 

1  -  wrt  wbl[t] 

2  -  wrt  wb2  [t] 

3  -  wrt  wb3[t] 

4  -  wrt  phia' [t] 

5  -  wrt  phis '  [t] 

rows  correspond  to  body  velocites 

1  -  vbi 

2  -  wbi 

3  -  vli 

4  -  wli 

5  -  vri 

6  -  wri 

*) 

pv  =  { 

{D [vbi , wbi [t] ] ,D [vbi , wb2 [t] ] ,D [vbi , wb3 [t] ] ,D[vbi,phia' [t]] ,D [vbi, phis' [t]]}, 
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{D [wbi , wbl [t] ] ,D [wbi , wb2 [t] ] ,D [wbi , wb3 [t] ] ,D[wbi,phia' [t]] ,D[wbi, phis' [t]]}, 

{D [vli , wbl [t] ] ,D [vli , wb2 [t] ] ,D [vli , wb3 [t] ] ,D[vli,phia' [t]] ,D[vli, phis' [t]]}, 

{D [wli , wbl [t] ] ,D [wli , wb2 [t] ] ,D [wli , wb3 [t] ] ,D[wli,phia' [t]] ,D[wli, phis' [t]]}, 

{D [vri , wbl [t] ] ,D [vri , wb2 [t] ] ,D [vri , wb3 [t] ] ,D[vri,phia' [t]] ,D[vri, phis' [t]]}, 

{D [wri , wbl [t] ] ,D [wri , wb2 [t] ] ,D [wri , wb3 [t] ] ,D[wri,phia' [t]] ,D[wri, phis' [t]]} 

}; 

(*  inertia  forces  and  inertia  torques  *) 

(*  all  vectors  are  expressed  in  body  coordinates.  See  pp26  of  notebook 
dated  6/4/93-  for  explanation  of  the  following,  ilb  and  irb  are 
inertia  of  arms  expressed  in  body  coord  therefore  they  are  not 
constant  in  body  coord,  system.  Therefore  in  taking  total 
time  derivative  of  ang.  mom.  need  to  take  time  derivative 
of  ilb  and  wbi  parts.  Since  everything  in  body  coord,  need 
to  take  cross  product  of  ang.  mom.  with  wbi. 

*) 

infb  =  mb  wbitilde.vbi  +  mb  D[vbi,t]; 

intb  =  D[ib,t]  .  wbi  +  ib  .  D[wbi,t]  +  wbitilde  .  ib  .  wbi; 

infl  =  ml  wbitilde. vli  +  ml  D[vli,t]; 

inti  =  D[ilb,t]  .  wli  +  ilb  .  D[wli,t]  +  wbitilde  .  ilb  .  wli; 

infr  =  mr  wbitilde. vri  +  mr  D[vri,t]; 

intr  =  D[irb,t]  .  wri  +  irb  .  D[wri,t]  +  wbitilde  .  irb  .  wri; 
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(*  generalized  inertia  forces  *) 

fil  =  pv [ [1 , 1] ] . inf b  +  pv [ [2 , 1] ] . intb  + 
pv [ [3 , 1] ] . inf 1  +  pv [ [4 , 1] ] . inti  + 
pv [ [5 , 1] ] . inf r  +  pv [ [6 , 1] ] . intr ; 

fi2  =  pv [ [1 , 2] ] . inf b  +  pv [ [2 , 2] ] . intb  + 
pv [ [3 , 2] ] . inf 1  +  pv [ [4 , 2] ] . inti  + 
pv [ [5 , 2] ] . inf r  +  pv [ [6 , 2] ] . intr ; 

fi3  =  pv [ [1 , 3] ] . inf b  +  pv [ [2 , 3] ] . intb  + 
pv [ [3 , 3] ] . inf 1  +  pv [ [4 , 3] ] . inti  + 
pv [ [5 , 3] ] . inf r  +  pv [ [6 , 3] ] . intr ; 

fi4  =  pv[[l,4]] .infb  +  pv [ [2 ,4] ] . intb  + 
pv [ [3 ,4] ] . inf 1  +  pv [ [4 ,4] ] . inti  + 
pv [ [5 , 4] ] . inf r  +  pv [ [6 , 4] ] . intr ; 

fib  =  pv [ [1 , 5] ] . inf b  +  pv [ [2 , 5] ] . intb  + 
pv [ [3 , 5] ] . inf 1  +  pv [ [4 , 5] ] . inti  + 
pv [ [5 , 5] ] . inf r  +  pv [ [6 , 5] ] . intr ; 

(*  External  forces  expressed  in  body  coordinates  *) 

(* 

positive  acting  on  left  limb,  neg  on  body  use  wli,wbi 
phi_rest  is  the  spring  rest  length,  should  be  less  than 
the  nominal  phil  in  order  to  make  taul  in  steady  rotation 
cancel  the  centrifugal  force. 
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*) 

ftaul  =  {-ksh  (phia[t]  +  phis[t])  -  bsh(phia' [t]  +  phis ' [t] ) , 0 , 0} ; 

ftaur  =  {-ksh  (phia[t]  -  phis[t])  -  bsh(phia' [t]  -  phis ' [t] ) , 0 , 0} ; 

(^generalized  external  forces  *) 

fel  =  -pv [ [2 , 1] ] . ftaul  +  pv [ [4 , 1] ] .f taul  -  pv [ [2 , 1] ] . ftaur  +  pv [ [6 , 1] ] . f taur ; 

fe2  =  -pv [ [2 , 2] ] . ftaul  +  pv [ [4 , 2] ] .f taul  -  pv [ [2 , 2] ] . ftaur  +  pv [ [6 , 2] ] . f taur ; 

fe3  =  -pv [ [2 , 3] ] . ftaul  +  pv [ [4 , 3] ] .f taul  -  pv [ [2 , 3] ] . ftaur  +  pv [ [6 , 3] ] . f taur ; 

fe4  =  -pv [ [2 ,4] ] . ftaul  +  pv [ [4 ,4] ] .f taul  -  pv [ [2 ,4] ] . ftaur  +  pv [ [6 ,4] ] . f taur ; 

fe5  =  -pv [ [2 , 5] ] . ftaul  +  pv [ [4 , 5] ] .f taul  -  pv [ [2 , 5] ] . ftaur  +  pv [ [6 , 5] ] . f taur ; 

(*  Expand  necessary  for  later  simplifications  *) 

efil  =  Expand[fil]; 
efi2  =  Expand [fi2]; 
efi3  =  Expand [fi3]; 
efi4  =  Expand [fi4]; 
efi5  =  Expand [fib]; 


efel  =  Expand [fel]; 
efe2  =  Expand [fe2]; 
efe3  =  Expand [fe3]; 
efe4  =  Expand [fe4]; 
efe5  =  Expand [fe5]; 
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A. 2  Mathematica  Code  for  Linearizing  Equations 
of  Motion 

Following  is  the  Mathematica  code  for  analytically  linearizing  the  non-linear  equations 
of  motion  of  the  previous  section. 

(* 

file  for  linearizing  manBdof _sym_antisym  model  for  following  case. 

here  the  arm  movement  is  decomposed  into  symmetric  and 
antisymmetric  components,  phil  =  x4  +  x5 ,  phir  =  x4  -  x5 
x4  is  the  antisymmetric  component 
x5  is  the  symmetric  component 

vector  from  left  limb  c.g.  to  shoulder  is  {0,0, bjS}  in  limb  coord, 
vector  from  right  limb  c.g.  to  shoulder  is  {0,0, bjS}  in  limb  coord, 
mass  of  left  and  right  limb  are  equal,  mr  =  ml  =  m 
inertia  of  left  and  right  limb  are  equal,  Il=Ir 
vector  from  body  c.g.  to  left  shoulder  is  {0 , ibj 2 , ibj 3}in  body  coord 
vector  from  body  c.g.  to  rt .  shoulder  is  {0 , -ibj 2 , ibj 3}in  body  coord 
nominal  body  rotation  rate  is  {0,W,0} 

left  arm  is  nominally  at  phiO,  right  arm  is  nominally  at  -phiO, 
deviations  in  both  arm  movements  are  of  equal  sign  and  magnitude. 

After  the  5  dof  linearized  model  is  complete,  a  constraint  equation 

of  the  form  x  =  B  z  can  be  applied  to  the  system  of  equations. 

B  represents  the  dependence  of  the  original  state  x  on  another 
state  vector  z.  z  can  be  of  lower  dimension  than  x.  Then  the 
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transformed  set  of  linear  equations  are  produced  as  follows: 

B^T  M  B  zddot  +  B^T  C  B  zdot  +  B^T  K  B  z 

This  transformation  can  be  applied  in  this  file  or  in  a  later 
one  that  evaluates  the  numerical  quantities. 

*) 

(*  for  yman  model  *) 

(* 

simplify  =  {bjll->0,  bjl2->0,  bjl3->bj3,\ 
bjrl->0,  bjr2->0,  bjr3->bj3,\ 
ibjll->0,  ibjl2->  ibj2,  ibj 13->ibj 3 , \ 
ibjrl->0,  ibjr2->  -ibj2,  ibj r3->ibj 3 , \ 
irll->illl,  ir22->il22,  ir33->il33,  mr->m,  ml->m,\ 
mt->mb+2m} ; 

*) 

(*  for  yman  model  with  no  shoulder  width  and  pvt  at  c.g.*) 

simplify  =  {bjll->0,  bjl2->0,  bjl3->bj3,\ 
bjrl->0,  bjr2->0,  bjr3->bj3,\ 
ibjll->0,  ibjl2->  0,  ibjl3->0,\ 
ibjrl->0,  ibjr2->  0,  ibjr3->0,\ 

irll->illl,  ir22->il22,  ir33->il33,  mr->m,  ml->m,\ 
mt->mb+2m} ; 


sfil  =  efil/ . simplify ; 
sfi2  =  efi2/ . simplify ; 
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sfi3  =  ef 13/ . simplify ; 
sfi4  =  ef 14/ . simplify ; 
sfi5  =  ef 15/ . simplify ; 

(*  Replace  angular  velocities  and  accelerations  with  corresponding 
euler  angle  expressions. 

*) 

S213  =  { 

{0  ,  Cos  [thl  [t]  ]  ,  Sin  [thl  [t]  ]  Cos  [th2  [t]  ]  }  , 

{0  ,  -Sin  [thl  [t]  ]  ,  Cos  [thl  [t]  ]  Cos  [th2  [t]  ]  }  , 

{1,  0,  -Sin[th2[t]]} 

}; 


C213  =  { 

{  Cos  [thl  [t]  ]  Cos  [th3  [t]  ]  +  Sin  [thl  [t]  ]  Sin  [th2  [t]  ]  Sin  [th3  [t]  ]  , 
Sin  [thl  [t]  ]  Cos  [th2  [t]  ]  , 

-Cos  [thl  [t]  ]  Sin  [th3  [t]  ]  +  Sin  [thl  [t]  ]  Sin  [th2  [t]  ]  Cos  [th3  [t]  ]  } , 
{  -Sin  [thl  [t]  ]  Cos  [th3  [t]  ]  +  Cos  [thl  [t]  ]  Sin  [th2  [t]  ]  Sin  [th3  [t]  ]  , 
Cos  [thl  [t]  ]  Cos  [th2  [t]  ]  , 

Sin  [thl  [t]  ]  Sin  [th3  [t]  ]  +  Cos  [thl  [t]  ]  Sin  [th2  [t]  ]  Cos  [th3  [t]  ]  }  , 

{  Cos  [th2  [t]  ]  Sin  [th3  [t]  ]  , 

-Sin  [th2  [t]  ]  , 

Cos  [th2  [t]  ]  Cos  [th3  [t]  ]  } 

}; 

thdot  =  {thl'[t],  th2'[t],  th3'[t]}; 
wba  =  {0,W,0}; 


wca  =  S213. thdot  +  C213.wba 
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wcadot  =  D[wca,t]; 

(*  See  Hughs  pp  116  and  my  notebood  6/14/93  pp  33  to  see  linearization 
of  a  steadily  rotating  system.  This  system  steadily  rotates 
about  the  2  axis  with  speed  W. 

the  linearized  states  are  not  consistent  with  Hughs  notation  so 
that  the  mass  matrix  retains  its  nominal  form. 

xl  tilt  angle  (th2) 

x2  dev.  from  nominal  sault  angle  (th3) 

x3  twist  angle  (thl) 

x4  dev  of  both  arms  in  anit-symmetric  mode  form  PhiO 

x5  dev  of  both  arms  in  symmetric  mode  form  PhiO 

*) 

eulerreplace  =  {wbl [t] ->wca [ [1] ] ,  wb2 [t] ->wca[ [2] ] ,  wb3 [t] ->wca [ [3] ] , 

wbl ’ [t] ->wcadot [ [1] ] ,  wb2 ’ [t] ->wcadot [ [2] ] ,  wb3 ’ [t] ->wcadot [ [3] ] } ; 

linname  =  {phia[t]  ->  x4,  phis[t]  ->  x5,  \ 
thl  [t]  ->  x3,  th2[t]  ->  xl,  th3[t]  ->  x2,  \ 

thl'Et]  ->  xdot3  W,  th2'[t]  ->  xdotl  W,  th3'[t]  ->  xdot2  W,  \ 
phia' [t]  ->  xdot4  W,  phis' [t]  ->  xdotB  W,  \ 

thl"[t]  ->  xddot3  (W  W),  th2"[t]  ->  xddotl  (W  W)  ,  th3"[t]  ->  xddot2  (W  W)  ,  \ 
phia' '  [t]  ->  xddot4  (W  W) ,  \ 
phis' '[t]  ->  xddotB  (W  W)}; 


lincond  =  {xdotl->0,  xdot2->0,  xdot3->0,  xdot4->0,  xdotB->0,\ 
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eulf il 
eulf i2 
eulf i3 
eulf 14 
eulf 15 

renf 11 
renf 12 
renf 13 
renf 14 
renf 15 

renf el 
renf e2 
renf e3 
renf e4 
renf e5 

mmll  = 
mml2  = 
mml3  = 
mml4  = 
mml5  = 

mm21  = 
mm2  2  = 


xl->0,  x2->0,  x3->0,  x4->0,  x5->phl0}; 

=  Expand [sf 11/ . eulerreplace]  ; 

=  Expand [sf 12/ . eulerreplace]  ; 

=  Expand [sf 13/ . eulerreplace]  ; 

=  Expand [sf 14/ . eulerreplace]  ; 

=  Expand [sf 15/ . eulerreplace]  ; 

=  Expand [eulf 11/ . llnname]  ; 

=  Expand [eulf 12/ . llnname]  ; 

=  Expand [eulf 13/ . llnname]  ; 

=  Expand [eulf 14/ . llnname]  ; 

=  Expand [eulf 15/ . llnname]  ; 

=  Expand [ef el/ . llnname] ; 

=  Expand [efe2/ . llnname]  ; 

=  Expand [efe3/ . llnname] ; 

=  Expand [efe4/ . llnname]  ; 

=  Expand [efe5/ . llnname] ; 

Coefficient [renf 11 ,xddotl] / . llncond; 
Coefficient [renf 11 ,xddot2] / . llncond; 
Coefficient [renf 11 ,xddot3] / . llncond; 
Coefficient [renf 11 ,xddot4] / . llncond; 
Coefficient [renf 11 ,xddot5] / . llncond; 

Coefficient [renf 12 ,xddotl] / . llncond; 
Coefficient [renf 12 ,xddot2] / . llncond; 
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mm23  =  Coefficient [renfi2 ,xddot3] /. lincond; 
mm24  =  Coefficient [renfi2 ,xddot4] /. lincond; 
mm25  =  Coefficient [renfi2 ,xddot5] /. lincond; 

mm31  =  Coefficient [renfi3,xddotl] /. lincond; 
mm32  =  Coefficient [renfi3,xddot2] /. lincond; 
mm33  =  Coefficient [renfi3,xddot3] /. lincond; 
mm34  =  Coefficient [renfi3,xddot4] /. lincond; 
mm35  =  Coefficient [renfi3,xddot5] /. lincond; 

mm41  =  Coefficient [renfi4,xddotl] /. lincond; 
mm42  =  Coefficient [renfi4,xddot2] /. lincond; 
mm43  =  Coefficient [renfi4,xddot3] /. lincond; 
mm44  =  Coefficient [renfi4,xddot4] /. lincond; 
mm45  =  Coefficient [renfi4,xddot5] /. lincond; 

mm51  =  Coefficient [renfiB jXddotl] /. lincond; 
mm52  =  Coefficient [renfiB ,xddot2] /. lincond; 
mmB3  =  Coefficient [renfiB ,xddot3] /. lincond; 
mmB4  =  Coefficient [renfiB ,xddot4] /. lincond; 
mmBB  =  Coefficient [renfiB ,xddotB] /. lincond; 

(*  this  gets  the  d.c.  forcing  part  *) 

del  =  renfil/ . lincond; 
dc2  =  renfi2/ . lincond; 
dc3  =  renfi3/ . lincond; 
dc4  =  renfi4/ . lincond; 
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dc5  =  renf 15/ . lincond; 


cl  =  renf 11/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} ; 
c2  =  renf 12/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} ; 
c3  =  renf 13/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} ; 
c4  =  renf 14/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} ; 
c5  =  renf 15/ . {xddotl->0 ,xddot2->0 ,xddot3->0 ,xddot4->0 ,xddot5->0} ; 


dcdxll 

dcdxl2 

dcdxl3 

dcdxl4 

dcdxl5 


D [cl ,xl] / . llncond; 
D [cl ,x2] / . llncond; 
D [cl ,x3] / . llncond; 
D [cl ,x4] / . llncond; 
D [cl ,x5] / . llncond; 


dcdx21 

dcdx22 

dcdx23 

dcdx24 

dcdx25 


D [c2 ,xl] / . llncond; 
D [c2 ,x2] / . llncond; 
D [c2 ,x3] / . llncond; 
D [c2 ,x4] / . llncond; 
D [c2 ,x5] / . llncond; 


dcdx31 

dcdx32 

dcdx33 

dcdx34 

dcdx35 


D [c3 ,xl] / . llncond; 
D [c3 ,x2] / . llncond; 
D [c3 ,x3] / . llncond; 
D [c3 ,x4] / . llncond; 
D [c3 ,x5] / . llncond; 


dcdx41  =  D [c4,xl] /. llncond; 
dcdx42  =  D [c4,x2] /. llncond; 
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dcdx43  =  D [c4,x3] / . lincond; 
dcdx44  =  D [c4,x4] /. lincond; 
dcdx45  =  D [c4,x5] /. lincond; 

dcdxBl  =  D [c5 ,xl] /. lincond; 
dcdx52  =  D [c5 ,x2] /. lincond; 
dcdx53  =  D [c5 ,x3] /. lincond; 
dcdx54  =  D [c5 ,x4] /. lincond; 
dcdx55  =  D [c5 ,x5] /. lincond; 

dcdxdotll  =  D [cl jXdotl] /. lincond; 
dcdxdotl2  =  D [cl ,xdot2] /. lincond; 
dcdxdotl3  =  D [cl ,xdot3] /. lincond; 
dcdxdotl4  =  D [cl ,xdot4] /. lincond; 
dcdxdotlB  =  D [cl jXdotB] /. lincond; 

dcdxdot21  =  D [c2 jXdotl] /. lincond; 
dcdxdot22  =  D [c2 ,xdot2] /. lincond; 
dcdxdot23  =  D [c2 ,xdot3] /. lincond; 
dcdxdot24  =  D [c2 ,xdot4] /. lincond; 
dcdxdot2B  =  D [c2 ,xdotB] /. lincond; 

dcdxdot31  =  D [c3 jXdotl] /. lincond; 
dcdxdot32  =  D [c3 ,xdot2] /. lincond; 
dcdxdot33  =  D [c3 ,xdot3] /. lincond; 
dcdxdot34  =  D [c3 ,xdot4] /. lincond; 
dcdxdot3B  =  D [c3 ,xdotB] /. lincond; 
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dcdxdot41  = 
dcdxdot42  = 
dcdxdot43  = 
dcdxdot44  = 
dcdxdot45  = 

dcdxdotBl  = 
dcdxdot52  = 
dcdxdot53  = 
dcdxdot54  = 
dcdxdotBB  = 


D  [c4,xdotl] / 
D [c4,xdot2] / 
D [c4,xdot3] / 
D [c4,xdot4] / 
D [c4,xdotB] / 

D [cB  jXdotl] / 
D [cB ,xdot2] / 
D [cB ,xdot3] / 
D [cB ,xdot4] / 
D [cB  jXdotB] / 


.  lincond; 
.  lincond; 
.  lincond; 
.  lincond; 
.  lincond; 

.  lincond; 
. lincond; 
.  lincond; 
.  lincond; 
. lincond; 


(*  this  is  the  linearization  of  the  external  forces  *) 


dtaudxdotll 

dtaudxdotl2 

dtaudxdotl3 

dtaudxdotl4 

dtaudxdotlB 


D [renf el ,xdotl] / . lincond; 
D [renf el ,xdot2] / . lincond; 
D [renf el ,xdot3] / . lincond; 
D [renf el ,xdot4] / . lincond; 
D [renf el ,xdotB] / . lincond; 


dtaudxdot21 

dtaudxdot22 

dtaudxdot23 

dtaudxdot24 

dtaudxdot2B 


D [renf e2 , xdot 1] / . lincond ; 
D [renf e2 ,xdot2] / . lincond; 
D [renf e2 , xdot3] / . lincond ; 
D [renf e2 , xdot4] / . lincond ; 
D [renf e2 ,xdotB] / . lincond; 


APPENDIX  A.  APPENDIX 


168 


dtaudxdotSl 

dtaudxdot32 

dtaudxdotSS 

dtaudxdot34 

dtaudxdot35 


D [renf e3 , xdot 1] / . lincond ; 
D [renf e3 ,xdot2] / . lincond ; 
D [renf e3 , xdot3] / . lincond ; 
D [renf e3 , xdot4] / . lincond ; 
D [renf e3 ,xdot5] / . lincond ; 


dtaudxdot41 

dtaudxdot42 

dtaudxdot43 

dtaudxdot44 

dtaudxdot45 


D [renf e4 , xdot 1] / . lincond ; 
D [renf e4,xdot2] / . lincond; 
D [renf e4 , xdot3] / . lincond ; 
D [renf e4 , xdot4] / . lincond ; 
D [renf e4,xdot5] / . lincond; 


dtaudxdotBl 

dtaudxdot52 

dtaudxdot53 

dtaudxdot54 

dtaudxdotBB 


D [renf eB , xdot 1] / . lincond ; 
D [renf eB ,xdot2] / . lincond; 
D [renf eB , xdot3] / . lincond ; 
D [renf eB , xdot4] / . lincond ; 
D [renf eB ,xdotB] / . lincond; 


dtaudxll 

dtaudxl2 

dtaudxl3 

dtaudxl4 

dtaudxlB 


D [renf el ,xl] / . lincond; 
D [renf el ,x2] / . lincond; 
D [renf el ,x3] / . lincond; 
D [renf el ,x4] / . lincond; 
D [renf el ,xB] / . lincond; 


dtaudx21 

dtaudx22 

dtaudx23 


D [renf e2 ,xl] / . lincond; 
D [renf e2 ,x2] / . lincond; 
D [renf e2 , x3] / . lincond ; 
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dtaudx24 

dtaudx25 

dtaudxSl 

dtaudx32 

dtaudx33 

dtaudx34 

dtaudx35 

dtaudx41 

dtaudx42 

dtaudx43 

dtaudx44 

dtaudx45 

dtaudxBl 

dtaudx52 

dtaudx53 

dtaudx54 

dtaudxBB 


D [renf e2 , x4] / . lincond ; 
D [renf e2 ,xB] / . lincond; 

D [renf e3 ,xl] / . lincond; 
D [renf e3 , x2] / . lincond ; 
D [renf e3 , x3] / . lincond ; 
D [renf e3 , x4] / . lincond ; 
D [renf e3 , xB] / . lincond ; 

D [renf e4,xl] / . lincond; 
D [renf e4 , x2] / . lincond ; 
D [renf e4 , x3] / . lincond ; 
D [renf e4 , x4] / . lincond ; 
D [renf e4 , xB] / . lincond ; 

D [renf eB ,xl] / . lincond; 
D [renf eB ,x2] / . lincond; 
D [renf eB , x3] / . lincond ; 
D [renf eB , x4] / . lincond ; 
D [renf eB ,xB] / . lincond; 
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A. 3  Parameters  of  the  Five  D.O.F.  Linear  Model 


In  this  section  I  present  the  algebraic  expressions  for  the  elements  of  the  linear  dy¬ 
namic  model  used  to  study  passive  layout  stability.  These  elements  were  computed 
from  the  symbolic  linearization  of  the  full  non-linear  dynamic  equations  of  the  hve 
d.o.f.  system.  The  equations  of  motion  and  the  associated  matrices  are  repeated  from 
Chapter  4  Equations  4. 1-4.6. 

MI  +  fl(G  +  D)i  +  n\K  +  K')  X  =  0 


where 


M 


G 


D 


mil  0  mi3  0  0 

0  m22  0  0  0 

mm  0  m33  0  0 

0  0  0  mil  0 

0  0  0  0  mss 

0  gi2  0  0  0 

-gi2  0  5-23  0  0 

0  ~5'23  0  0  0 

0  0  0  0  5'4s 

0  0  0  — t/rs  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0 
0  0  0  0  0 

0  0  0  0 
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kii  0  ki3  0  0 

0  7^22  0  0  0 

K  =  k^i  0  k33  0  0 

0  0  0  0  0 

0  0  0  0  7^55 

0  0  0  0  0 

0  0  0  0  0 

K'  =  0  0  2k,hlkl‘^  0  0 

0  0  0  0  0 

0  0  0  0  2kshlkf 

The  state  vector,  x,  is  comprised  of  (in  order):  fl  -  the  tilt  angle,  7  -  the  twist  angle, 
cf)a  -  the  asymmetric  deviation  of  the  two  arms  from  the  nominal  arm  angle,  a  -  the 
deviation  of  the  somersault  angle  from  the  frame  steadily  rotating  at  rate  0,  and 
-  the  symmetric  deviation  of  the  two  arms  from  the  nominal  arm  angle. 

The  individual  elements  of  the  above  matrices  are  given  below. 


^  ^  .  9  I  rb 9^  ml  nil, 

[Ibi  2  III  +2  rb^  m;  +  2  rbg  mi  +  — - ^  ^  — 

Srbs'^mf  ,  ,  ,  16  r/5  mf  cos  ( (Po ) 

- - ^ - Iris  rbs  mi  cos{<Po) - 7- - 72 - 

2  m;  +  mj  (2  m;  +  mj) 

8  rls  rbs  mf  mi  cos{(I)o)  ,  16  r/5  cos(  (Po )  .  ,2  1  ^  i2  . 

(2  mi  A  nil)  2  mi  A  nii 

8  r/5^  mf  cos(<Po)^  4  r/5^  mf  nii  cos(<Po)^ 

(2m/  +  mj)^  (2m/  +  mj)^ 

8  rU^  ml  cos{<I>o  Y  ,  ,  .  ^  ^  n  ,9  .  ^  ^  n9n 

- h  4  rh  rbs  mi  sm(  <I?o  )  A  2  rls^  mi  sm(  <I>o )  ) 


mi3 
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m22 


mss 


77144 


mss 


8  mf  cos{$oY  4  r/5^  mj  cos(<Po)^  8  rlg"^  mf  cos{$oY 
(2  m;  +  TTlfi  Y  (2  m;  +  TTlfi  Y  2  m;  +  mj 

2  rls  rb2  mi  sin(<Po)  +  2  r/5^  m;  sin{$oY) 


[lbs  +  2  m;  +  2  IIs  cos{$oY  +  4  rls  rbs  mi  sin((Po)  + 
2  /4  sin((Po)^  +  2  r/5^  m;  sin((Po)^) 


,  „  „  ,9  /  *  x2  8  r/5^  mf  cos(<Po)^ 

(2//;  +  2  r/5  m;cos(<?o)  4 - - - “2 - h 

(2  m;  +  mt) 

Irls^  mf  mi,cos{<Pof  8  rls^  mf  cos{<Pof  ,  ,  ,2  •  x2x 

- 7- - ^ - 72 - 4 ^ +  m;sin(<?o)  ) 


(2  m;  +  mj)' 


2  mi  +  mi, 


(lbs  +  2  m;  + 


8rb3‘^mf  Irbs^m^mt  8rb3‘^mf 


(2m;  +  mj)^  (2m;  +  mj)^  2  m;  +  mj 


16  rls  rbs  mf  cos{(I)o)  8  rls  rbs  mj  mi  cos{<I>o)  , 
4r/5  rbs  mi  cos(  <Po) - 7:^ - 79^ - 77^ - 79^ - h 


(2  m;  +  miY 


(2  m;  +  miY 


16  rh  rbs  m?  cos(  (In)  /,x2  ^,9  /,x9 

- - - \- 2  II2  cos{  (pg  )  A  2  rls  micos{  (pg  )  + 

2  m;  +  mj 

8  rlg"^  mf  cos{$g)‘^  4  rls^  mf  mi  cos{I>g)^  8  rlg"^  mf  cos{$g)‘^ 

72  I  71  79  77  ^  1“ 


(2  m;  +  mi)' 

2  II3  sin((Po)^) 


(2  m;  +  mi)' 


2  mi  +  mi 


=  (2  III  A  2  rl 3'^  mi  cos{  do  f  +  2  rls ^  mi  sin(  dg)^  A 
4  rlg"^  mf  mi  sin{dg)^  8  rlg"^  mf  sin{dg)^ 


8  rlg^  mf  sin(<Po)^ 
(2  m;  +  mi)^ 


A 


(2  mi  A  mi)' 


2  mi  A  mi 


[Ibi  —  lbs  A  lbs  A  2  III  +4  rbs^  mi  —  2  II2  cos{dg)‘^  A 
2  II3  cos( dg)"^  A  8  rls  rbs  mi  sin( dg)  A  2  II2  sin( dg)^  —  2  II3  sin( dg)^  A 
4  rlg"^  mi  sin{dg)^) 


gi2 
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§23  =  (-2  III  +  2  Ils  cos(I>o  f  -  2  Ils  cos(I>o  f  -  4  rls  rbs  mi  sin(I>o) - 
2  II2  sin{<I>oY  +  2  Ils  sin(<Po)^  —  4  r/5^  m;  sin(<Po)^) 


545  — 


•  \  ,  16  r/5  mf  sin ((Po)  8  r/5  mj  sin( (Po ) 

(4  rls  rbs  mi  sm(  <Po)  -\ - ^ - -2 - h  - 


{2  mi  +  nitY 


{2  mi  +  nitY 


16  rls  rbs  rnf  sin(<Po) 

2  mi  +  mi, 


—  Ills  cos{I>o)  sin(<Po)  +  4  7/5  cos{I>o)  sin(<7o) — 


Irls'^mi  cos(<7o)  sin(<7o)  — 


IQrls^mf  cos{$o)  sin(<7o) 

{2  mi  +  miY 


8  r/5^  mf  mi  cos{$o)  sin(<7o)  ^  16  r/5^  mf  cos{$o)  sin(<7o) , 


(2  m/  +  miY 


2  mi  +  mi 


kii  =  {lbs  -  lbs  -2rb2‘^mi +  2rbs‘^mi  + 


S  rbs'^  mf  ^  Irbs^  mf  mi 


{2mi  +  miY  {2mi  +  miY 

16  rls  rbs  rnf  cos(<7o) 


Srbs^mf 

- - ^ - 4  rls  rbs  mi  cos{  <Po) - - - -  2 

2  mi  A  mi  (2  mi  A  mi) 

8  rls  rbs  mf  mi  cos(<7o)  16  rls  rbs  mf  cos(<7o) 

2  ^ 


{2  mi  A  mi) 


2  mi  A  mi 


A  2  II2  COs(  (pg)  ~ 


m  ,9  /  n9  8  rh^  mf  cos(  (In) 

2  Ils  cos{(I)o)  +2r/5  m/cos((7o)  H - - - 

(2  m/  +  mi) 

4  rls"^  mf  mi  cos{$o)‘^  8  r/5^  mf  cos{$o)‘^ 


A 


(2m/  +  mj)^  2m/  +  mj 

2  II2  sin(<7o)^  +  2  Ils  sin(<7o)^  —  2  r/5^  m/  sin(<7o)^) 


—  4  rls  rbs  mi  sin(<7o)  — 


^13  — 


8  rls  rbs  mf  cos{<Po) 

(-2  rls  rbs  mi  cosfflg) - ; - ^2 - 

(2  m/  +  mi) 

4  rls  rbs  mf  mi  cos(<7o)  ^  8  rls  rbs  mf  cos(<7o) 


(2  m/  +  mj)' 


2  m/  +  mj 


+  2  cos(  )  — 


,9  ,0  .  ^9  8  rln^  mf  COs(  (pn)^ 

2  Ils  cosffPo)  A2rls  micosffPo)  H - ; - ^  .2  + 

(2  m/  +  mi) 

4  r/5^  m^  mi  cos{$o)‘^  8  r/5^  m^  cos(<7o)^ 

(2m/  +  mj)^  2m/  +  mj 

2  //g  sin(<7o)^  +  2  Ils  sin(<7o)^  —  2  r/5^  m/  sin(<7o)^) 


—  2  rls  rbs  mi  sin(<7o)  — 
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k22  =  (  —  Ibi  A  lbs  —  2  Ik  —  2  nil  A  2  Ils  cos(<Po)^  —  I  rig  rbs  nil  sm(<Po)  E 
2  Ils  sin{<I>oY  —  2  rlg"^  nii  sin(<Po)^) 


kss  — 


(-2  rl,  rh,  ,n,  cos{t,]  -  « ’’fa cos(j.»)  _  4  ri,  rt,  m|  m,  cos(j.») 


(2  m/  +  mj) 


(2  m/  +  mj)' 


+ 


8  rh  rb-i  nij  cos ((hn)  ,  .^9  ,  ,9  .^9 

- - - -  A  2  II2  cos{(I)o)  -2II3  cos{<Po)  +  2  r/5^  m/ cos((Po)  + 

2  m/  +  mj 

8  r/5^  mf  cos(<Po)^  ^  4  r/5^  m^  mj  cos(<Po)^  8  r/5^  m^  cos(<Po)^ 


(2  m/  +  mj)' 


(2  m/  +  mj)' 


2  m/  +  mj 


2  //g  sin(<Po)^  +  2  Ils  sin(<Po)^  —  2  r/5^  nii  sin(<Po)^) 


^55  — 


(-2  ri.,  ri.,  cos(4.„)  -  cos(j.»)  _  4  ,4,  rt,  m|  m,  cos(j.») 


(2  m/  +  mj) 


(2  m/  +  mj)' 


+ 


8  rh  rbi  nij  cos  {$o)  ,  /  n9  ,  /  n9  ,9  /n9 

- - - ^  +  2 /4  cos((Po)  -2II3  cos{<Po)  +  2  r/5^  m/ cos((Po)  + 

2  m/  +  mj 

8  rk"^  nif  cos{$oY  4  r/5^  m^  mj  cos(<Po)^  8  r^"^  nif  cos{$oY 
(2m/  +  mj)"  ^  ^  ^ 


{2  nil  +  niiY 


2  nil  A  nii 

8^/  2^3, 


2  /4  sin((Po)^  +  2  7/5  sin((7o)^  -  2  r/5^  m/  sin((7o)^  -  ^  ^^0 ) 

(2  m/  +  nii) 

4  r/5^  mf  nii  An{$oY  8  r/5^  mf  sin(<7o)^ 

(2m/  +  mj)^  2niiA  mb 
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A. 4  Analytic  Solution 

Solving  for  the  torque-free  motion  of  a  rigid  body  tumbling  in  space  is  one  of  the 
classic  problems  of  dynamics.  Given  the  limited  value  of  numerical  integration  as 
a  means  of  studying  dynamics  problems  before  the  advent  of  the  digital  computer, 
much  effort  was  applied  to  hnding  an  analytic  solution  to  the  governing  equations  of 
motion.  The  analytic  solution  provides  the  angular  velocity  of  the  body  in  a  body 
hxed  principle  axes  system  and  the  attitude  of  the  principle  axis  system  with  respect 
to  an  inertial  coordinate  system  [Hughes  86].  The  contribution  that  these  analytic 
solutions  offer  today  is  a  concise  description  of  the  minimal  number  of  states  and 
non-dimensional  parameters  that  govern  rigid  body  rotation.  The  analytic  solutions 
also  give  rise  to  elegant  geometric  interpretations  of  the  motion  that  helps  provide 
some  intuition  of  this  otherwise  complex  movement. 

I  present  the  analytic  solution  for  a  tri-inertial  body.  A  tri-inertial  body  is  one 
in  which  each  of  the  principal  inertias  is  distinct.  Without  loss  of  generality  for 
the  rest  of  this  section  I  will  insist  that,  A  >  I^.  The  analytic  solution  for 

rotation  of  a  tri-inertial  body  uses  the  Jacobian  elliptic  functions.  Elliptic  functions 
are  a  generalization  of  the  circular  functions,  sine  and  cosine,  that  originated  from 
an  effort  to  compute  the  circumference  of  an  ellipse  [Bowman  61].  The  Euler  angles 
describing  the  body  attitude  with  respect  to  inertial  space  can  also  be  written  using 
elliptic  functions. 

Two  integrals  of  motion  are  essential  in  the  analytic  solution  for  torque  free  rigid 
body  motion.  These  represent  conservation  of  angular  momentum  h  and  conservation 
of  kinetic  energy  T. 


12^2  ^  ^  j2^2 

(A.l) 

+  12^2  T  -^3^3 

(A.2) 

Eor  a  hxed  angular  momentum,  the  kinetic  energy  of  rotation  is  bounded  above  and 
below  by  the  kinetic  energy  of  two  pure  spin  solutions.  The  highest  possible  energy 
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state  corresponds  to  rotation  about  the  minor  principal  axis  and  the  lowest  energy 
state  corresonds  to  rotation  about  the  major  principal  axis.  A  continuum  of  solutions 
of  rotational  motion  of  a  rigid  body  exist,  each  corresponding  to  a  different  energy 
level,  between  these  two  bounds. 


<T  < 


The  qualitative  nature  of  the  solution  is  divided  into  two  regions  between  these  bounds 
and  separated  by  the  solution  of  pure  spin  about  the  middle  principal  axis.  The 
following  parameter  is  useful  in  distinquishing  between  the  two  qualitatively  different 
solutions. 


I 


/  has  units  of  inertia.  The  bounds  on  /  corresponding  to  the  bounds  on  energy  are 


h>I>h 


/  =  /i,  /  =  /2,  and  I  =  I3  correspond  to  the  pure  spin  solutions. 

The  maximum  amplitude  of  the  body  hxed  angular  velocities  can  be  obtained 
directly  from  the  two  integrals  of  motion. 


^Im 


^3m 


^2m 


^2m 


^2m 


h 

h 

h 


I-I3 


1/2 


Ihih 

h- 


lh{h 

h- 


Ihih  - 12) 


h 

h 


Ihih  -  h) 


(A.3) 

(A.4) 

{Ii>I>  h)  (A.5) 

(/  =  h)  (A.6) 

{h>I>  h)  (A.7) 
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Next,  the  analytic  solution  for  the  angular  velocity  is  given  for  the  case  of  /i  >  /  >  /2 


CUi 

=  SiUiradn^T]  k) 

(A.8) 

iL>2 

=  S2t02mSn{T-,k) 

(A.9) 

103 

=  S3UJ3^cn{T;  k) 

(A.IO) 

T 


k 


(h-hKi-h) 

1I1I2I3 


1/2 

(i  -  io) 


_(/l_/2)(/_/3)_ 


and  similarly  for  I2  >  I  A  I3 


(A.ll) 

(A.12) 


CUi 

=  SiLOi^cn{T]k) 

(A.  13) 

iL>2 

=  S2t02mSn{T-,k) 

(A.14) 

103 

=  S3UJ3^dn{T;  k) 

(A.  15) 

T 


[(/2-/3)(/l-/) 

II1I2I3 


1/2 

(f  -  to) 


,  ^  \ih  -  hVJ  - 


(A. 16) 

(A.17) 


The  coefficients  Si,S2  and  S3  are  all  ±1  such  that  S1S2S3  =  —1,  allowing  four  com¬ 
binations  of  signs.  The  Jacobian  elliptic  functions,  (sn,cn  and  dn),  depend  on  the 
scaled  time,  r  and  the  modulus  k.  When  k—  >  0  sn,  cn,  dn—  >  sin,  cos,  1. 

To  complete  the  solution  of  the  tumbling  motion  we  require  the  attitude  of  the 
body  as  a  function  of  time.  Since  angular  momentum  is  conserved,  the  inertially 
hxed  angular  momentum  vector  provides  the  reference  with  which  to  measure  body 
attitude.  The  components  of  angular  momentum  as  measured  in  the  body  hxed  axis 
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system  is  given  by 


hi 

IlUi 

h2 

= 

l2^2 

hs 

^3^3 

(A.18) 


Another  expression  for  this  vector  can  be  obtained  by  projecting  the  inertially  hxed 
angular  momentum  vector  onto  the  body  hxed  axes  using  the  matrix  Ca  (2.2).  Since 
the  orientation  of  the  inertial  coordinate  system  is  arbitrary,  we  choose  it  to  oriented 
such  that  the  angular  momentum  vector  is  oriented  intirely  along  the  inertial  y  axis. 
With  this  assumption  in  mind,  equating  the  components  of  these  two  expressions  for 
angular  momentum  gives 


(A. 19) 


h^i 

S^C&h 

l2^2 

= 

C^C&h 

^3^3 

-Seh 

Solutions  for  two  of  the  Euler  angles  are  available  immediately  as 


0 


•  -1 

^3^3 

Sin 

.  h  \ 

tan“^ 

'hup 
.12^2  - 

(A. 20) 

(A.21) 


The  solultion  for  $  requires  signihcantly  more  work.  First  we  use  equation  (2.4)  to 
substitute  for  Ui  in  (A.  19). 


/i(Ai])C*0$  +  C*i])0)  =  S\jiCsh 

l2{CAC&^  —  Ai])0)  =  CAC&h 
Is^  —  Sq^  +  T)  =  —SqU 


Solving  for  $  from  the  hrst  two  of  these  equations  yields 


$  = 


HhSl  +  hCl) 
hh 


(A. 22) 
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To  proceed  we  need  an  expression  for  and  S'!, 


Using  equation  A. 21,  we  obtain 


tan 


(hcoir 


Using  trigonometric  identities  we  obtain 


sin^  T 
cos^  T 


(hcoir 

(/lCUi)2  +  (I2i02)^ 

(-^2^2)^ 

(/lCUi)2  +  (l2C02)^ 


(A. 23) 
(A. 24) 


Using  these  in  equation  A. 22  we  obtain 


$ 


+  -^2^2 

(/lCUi)2  +  (/2CU2)2 


Using  equation  2.15  for  the  denominator  and  equation  2.16  for  the  numerator,  $  can 
be  written  as 


$ 


2T  -  htol 

-  ihiOsy 


At  this  point  we  have  to  specialize  the  solution  for  one  of  the  two  qualitatively 
different  cases.  First  we  solve  for  the  case  I\  >  I  >  l2-  Substituting  for  023  from 
equation  A.  10 

W  -  Ii4‘^lraCn‘^  . 

When  we  use  the  identity  +  cn^  =  1  we  obtain. 


Li  +  L2  stAt 
Ls  +  L4  stAt 


(A. 25) 


Li  —  h{2T  —  htol^)  (A. 26) 

L2  =  hlsLo'^^  (A. 27) 


where 
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Ls  =  h^-lltoL  (A.28) 

iv4  =  (A. 29) 


Integrating  the  equation  above  for  ^  will  provide  the  somersault  angle  as  a  function 
of  time.  This  integral  can  be  written  in  terms  of  the  Legendre  elliptic  integrals. 

In  a  similar  fashion  we  can  solve  for  the  somersault  angle  for  the  second  case 
I2  >  I  >  I3.  In  this  case  substitute  cos  from  equation  A. 15  into  A. 4  and  use  the 
following  identity:  dn^  =  1  —  k^srP .  This  results  in  the  same  form  for  $  with  new 
coefficients. 


Li  +  L2  stAt 
L3  +  L4  stAt 


where 


-  h{2T  -  h^^lm) 

(A. 30) 

L2 

hXuj^^k 

(A.31) 

L3 

7-2  2 

—  h 

(A. 32) 

L4 

r2,  ,2  7  2 

(A. 33) 

To  integrate  equation  A. 25  we  split  the  integrand  into  two  parts  to  obtain 


^dt 


^  Li  r  dt  L2  pI  n  sn 

L3  J  l  +  nsn^  L^n  Jm  l  +  nsn^ 


where  n  =  Lij L3.  We  would  like  to  integrate  this  equation  with  respect  to  the  scaled 
time  T  =  Crt  (A.  11  or  A.  16)  since  sn  =  sni^r).  Therefore,  we  compute  the  differential 
dt  =  jX  and  substitute  it  into  the  integral 

Li  I'  dr  L2  yT  nsn^dr 

LsCr  J  1  +  n  sn^  L^Cr  n  Jto  1  +  n  sn^ 


To  hud  the  change  in  somersault  angle  in  a  set  period  of  time  we  dehne  the  limits  of 
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integration  to  be  the  initial  time,  to,  and  the  hnal  time,  tj  or  equivalently  Tq  and  tj. 
The  change  in  somersault  angle  is  then  given  by  the  following  integral 


Li  rf  dr  L2  Ef  nsn’^dr 

^  °  LsCt  J^o  1  +  n  sn"^  LsCt  n  Jtq  1  +  n  sn"^ 


These  integrals  can  be  expressed  as  standard  elliptic  integrals 


/o  1  +  n  srA{T) 


=  Il{k,  n,  cf)) 


r  n  srA{T)dT 
lo  1  +  n  sn^(T) 


=  T  —  Il{k,  n,  cf)) 


where  sin  ^  =  stit.  Therefore,  breaking  up  the  limits  of  integration  we  obtain. 


Li  Ef  dr  Eo  dr  L2  Ef  nsn’^dr  Eo  nsrAdr 

^  °  LoCr  Jo  1  +  nsrP  Jo  1  +  nsn'^  LoCrU  Jo  1  +  nsn'^  Jo  1  +  nsrP 


Ef  nsrAdr  E°  nsrAdr 


which  is  written  compactly  as. 


«/-«o  =  [n{k,  n,  (ff)  -  U{k,  n,  ^o)]+  .  [(t/  -  U{k,  n,  cff))  -  (tq  -  U{k,  n,  ^0))] 


To  compute  these  integrals  we  must  solve  for  the  initial  scaled  time,  Tq,  for  a  given 
initial  state.  To  this  end,  the  angular  velocity,  ijj2,  is  given  by  the  elliptic  function  sn 


^2(t)  =  S2i02mSn{T-,  k) 


For  a  given  value  of  0J2I {E>2m  ^2)  we  need  to  compute  r.  If  we  assume  that  r  =  0 
when  UJ2  is  zero  and  increasing  then  by  dehnition  S2  =  1.  The  inverse  of  the  Jacobian 
elliptic  function,  sn,  is  the  Jacobian  elliptic  integral  of  the  hrst  kind. 
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Computation  of  this  integral  allows  us  to  compute  Tq  for  a  given  initial  iij2. 
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